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PREFACE 



In this book authors bring out the innovative applications of 
matrices defined, described and developed by them. Here they 
do not include the natural product on matrices newly described 
and defined by them in the book on ‘natural product x n on 
matrices’. 

This book is organized into seven chapters. The first one is 
introductory in nature. In the second chapter authors give the 
unique and new way of analyzing the data which is time 
dependent. We construct three types of matrices called Average 
Time Dependent data matrix (ATD matrix), Refined Time 
Dependent Data matrix (RTD matrix) and Combined Effective 
Time Dependent Data matrix (CETD matrix). The authors 
describe the working of this new type of matrix model by 
working with the real world transportation problem. It is proved 
this new model is effective and elegant. At this juncture the 
authors deeply acknowledge the unostentatious help rendered by 
Dr. Mandalam. 

In chapter three the authors for the first time define the new 
concept of matrices with linguistic variables. These linguistic 
matrices arc used in the place of fuzzy matrices which serve as 
the dynamical system for FCM and FRM models. By this 
method the solution themselves arc linguistic terms so 
mathematical interpretation is not needed. 
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In chapter four authors have used super matrices in the 
following ways. They are used in super fuzzy models like super 
FCMs, super FRMs super FAMs etc. Super linear algebra using 
these super matrices pave way to super eigen values and super 
eigen vectors, infact in almost all applications of matrices in 
linear algebra. The special mention is made to the Leontief 
open production super model and closed super model. By this 
method easy comparison and bulk working is possible. 

We have also used super matrices with entries from finite 
fields in the construction of several types of super codes. 

In chapter five we define interval matrix and matrix interval 
using the natural class of intervals, that is intervals of the form 
[a, b] where a < b or a > b or a and b not comparable. Using 
these interval matrices we have built interval linear algebra and 
these interval matrices when used as stiffness matrices of any 
mechanical problem easily yield a solution. 

Suggestions and literature how we have used and 
constructed DSm matrices and DSm super matrices is given. 
However the authors have constructed DSm vector spaces and 
DSm super vector spaces. The eigen values and eigen vectors of 
these spaces will be refined labels and refined label vectors, 
which comprises chapter six of this book. 

In the final chapter we have defined the new notion of 
bimatrices and n-matrices. Using these new concepts we have 
constructed linear algebra of type I and linear algebra of type II. 
Using these matrices we can build n-eigen values and n-eigen 
vectors. 

Further these n-matrices are used in the construction of 
multi expert fuzzy models like n-FCMs, n-FRMs, n-FAMs etc. 
Finally we have built several types of n-codes using these n- 
matrices whose entries are from finite fields. 

We thank Dr. K.Kandasamy for proof reading and being 
extremely supportive. 



W.B.VASANTHA KANDASAMY 
FLORENTIN SMARANDACHE 
INDRAVENKATBABU 
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Chapter One 



Inttoductiqn 



This book gives the contributions of the authors in the 
innovative applications of matrices to various fields. 

The authors were first to use matrix theory in the analysis 
of raw data. This is described in chapter two of this book. The 
method of applying matrices in this way makes one know the 
results and conclusions by just looking at it. Thus the final data 
is given a graph representation. 

Secondly the authors have constructed matrices where 
entries are just linguistic terms which arc orderable. For more 
about this refer [53]. Using these linguistic values we construct 
matrices. Infact these linguistic matrices are used in studying 
fuzzy linguistic models [53]. Thus we have in the first place 
constructed linguistic matrices and used them in mathematical 
models to study socio economic problems. 

The authors have used super matrices in the first place to 
build super linear algebras. Almost all classical results of linear 
algebras are derived for the super linear algebras [45]. Using 
these super matrices we have built several types of super fuzzy 
models and these super fuzzy models can yield the result in a 
shorter time and comparision of the results are ready made [48]. 




Introduction 



Further these matrices are used in the construction of 
algebraic super codes of different types. For more about 
supercodes please refer [50] . 

We build matrix of refined labels and using these matrices 
we construct DSm vector spaces of refined labels. Several 
properties about DSm vector spaces build using matrix of 
refined labels in elaborately discussed. This is yet another 
means of applying matrix theory to DSm vector spaces of 
refined labels [46]. 

We also construct DSm super vector space of refined labels. 
For these new structures the notion of supcrmatrices is used. In 
the supermatrix if we replace the real entries or any other entries 
by refined labels L we call them as supermatrix of refined labels 
and also they form vector spaces. This is carried out in chapter 
six of this book [52], 

Finally we build a natural class of intervals { [a, b] I a > b 
a < b or a = b, a, b e Q or C or Z or C(Z n ), Z n or R}. Using 
these natural classes of intervals we build interval matrices. 
These interval matrices behave like ordinary matrices as the 
entries are from this special natural class of intervals. Further 
study of stiffness matrix, load vector, mass matrix and damping 
matrix using interval for in this case lower and upper bound 
matrices is based on interval algebra. Also we find for square 
interval matrices, interval eigen values and interval eigen 
vectors [55]. 

Finally we define the notion of bimatrices and n-matrices. 
This is yet another innovative method of applying matrices. 
These n-matrices (n > 2) can be used in the construction of 
fuzzy n-models like n-FCMs, n-FRMs, n-BAMs etc., and their 
mixed structures [54, 56]. Finally we use these n-matrices in the 
construction of different types of n-codes [49]. 




Chapter Two 



Average Teve eepeneent (Aim 
B\taMxrix 



The raw data under investigation is classified under four 
broad heads viz; total number of passengers, total collection, 
number of trips, and hourly occupancy. These four broad heads 
form the columns of the matrices. The time periods of the day 
are represented by the rows of the matrices. Estimating the 
utility rate of a route is a five-stage process. In the first stage, 
we give a matrix representation of the raw data. Entries 
corresponding to the intersection of rows and columns are 
values corresponding to a live network. The initial M x N 
matrix is not uniform i.e., the number of individual hours in 
each time period may not be the same. So, in the second stage, 
we in order to obtain an unbiased uniform effect on each and 
every data so collected, transform this initial matrix into an 
Average Time Dependent Data matrix (ATD matrix). To make 
the calculations easier and simpler, we in the third stage, using 
the simple average techniques, convert the above time 
dependent data matrix into a matrix with entries e,j where, e,, e 
{- 1 , 0 , 1 }. 
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We name this matrix as the Refined Time Dependent Data 
matrix (RTD matrix). The value of the e,j corresponding to each 
entry is determined in a special way. At the fourth stage, using 
the refined time dependent data matrices, we get the Combined 
Effect Time Dependent Data matrix (CETD matrix) which gives 
the cumulative effect of all these entries. In the final stage, we 
obtain the row sums of the combined effect time dependent data 
matrix. We have written a Java program, which estimates all 
these five stages. [41] 

Using the raw data available for any transport corporation, 
we analyze the raw data via matrices and (i) predict the 
maximum utilization time period (peak hours) of a day; and (ii) 
estimate the overall utility rate of the routes. Thus, to be more 
precise our chief problem here is prediction of the peak hours 
and estimation of the most utilized routes using the raw data 
available from any transport organization. We have established 
that the results which we have predicted using the raw data on 
the route 18R which ply from Parry s to Dharmaraja koil 
coincides with the peak hours and the utilization rate of the 
route which can be obtained from the estimation of the observed 
data. Thus, our analysis not only predicts the peak hours and the 
maximum utility routes but also estimates for each route, time 
periods where buses need not be operated as operation of buses 
in those time periods will result in total loss. 

By total loss, we mean that the money collected from the 
passengers at that specified time period may not be even 20% of 
the cost of petrol and service charges spent on plying that 
service. So, even if the concern using the predictions on peak 
hours gets a profit of say 70%, if they do not comply with 
stopping of buses at the non-utilized hours, the transport 
corporation will still result only in loss. Hence, our approach 
via matrices not only suggests the peak hours and the maximum 
utilized routes but also predicts time periods and the routes in 
which the services should be totally avoided or reduced from 
their usual schedule in order to save themselves from loss. 
Thus, by adhering to our predictions the corporation might 
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curtail the loss incurred by the operation of services in the non- 
utilized time periods. 

The predictions can be explicitly seen even by a layman just 
by the observation of the combined effect time dependent data 
matrix mentioned earlier, which will be described and discussed 
in detail in the following. 

Hence, using these results the Transport Corporation can 
operate more number of buses at the peak hours and on the 
routes that are highly utilized and also stop operating buses at 
odd hours or reduce the frequency at the odd hours of the day. 
Thus, the Transport Corporation can achieve a meager gain or 
curtail the loss by using these predictions. 

We have considered the city of Madras, where the Pallavan 
Transport Corporation (PTC) is one of the transport 
organizations that caters to the demands of the community. The 
Pallavan Transport Corporation is a service oriented 
organization run by the state government for the people of the 
city and the neighbouring panchayats in the immediate vicinity 
of the city. It has a fleet of 1343 buses, and covers an overall of 
283 routes and satisfies the needs of more than 75% of the 
urban transit commuters. Apart from the regular services, it has 
the partial services, the night services and the special services. 

The regular services are of four types: normal, limited stop 
service, point-to-point service and the express service. Although 
the source, the destination and the path of travel is common to 
all the regular - services, these services considerably differ in the 
travel time, the travel fare, the comfort and the convenience of a 
passenger. The basic difference observed in the above types of 
services is primarily due to the number of intermediate stops the 
vehicles halt in their course of travel. 

The city transport corporation connects a vast network inter- 
linking various routes. The network has a set of nodes and a set 
of links connecting these nodes. Buses connect these nodes and 
ply over the links and help in passenger and freight 
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transportation. In the city network, we observe that some routes 
have a high patronage almost throughout the day, while other 
routes are busy only at particular hours of a day and there are 
still some routes, which have only very meager passenger 
support throughout the day. Here, we not only predict the high 
utility and peak hours but also suggest the stoppage of meager 
passenger supported services. 

To test the efficiency of the program, we have applied the 
algorithm to a few and established that the observed peak hours 
coincide well with the estimated peak hours for each route 
individually. 

In the first stage of the problem, we split the total hours of a 
day into various time-periods viz. morning, morning peak, 
evening. Each of these time periods consists of individual hours 
of the day viz. {Hi, H 2 , H 3 , ..., H 22 , H 23 , H 24 }, where H ; refers to 
the hour ending of the day. The working hours of the transport 
corporation under consideration is from {H 5 , H 6 , H 7 , .... H 2 i, 
H 22 }. The data obtained from the transport corporation is for 
each hour ending. Hence, the time interval between H, and H i+ i 
is 60 minutes. If we have all the related data viz. the total 
collection, the total number of passengers, the occupancy rate 
and the number of trips etc., for say every 30 minutes or for 
every 15 minutes, then the accuracy in the identification of the 
peak hours will further increase. 

The finer the intervals are divided the better is the 
prediction. This is also established in this chapter. Each time 
period consists of a set of H; individual hours, for example: the 
morning peak time period consists of {H 8 , H 9 , Hio}. These 
time -periods can vary from one analyst to another and from one 
route to another. Even, if we have each individual hour as a 
time period, still the total number of time periods must not 
exceed to total working hours of a day. 

Let M represent the total number of time periods under 
consideration. The utility rate of a time period is influenced by 
the various attributes acting on them. Let N represent the 
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number of attributes considered for the analysis. The various 
time periods are treated as rows and the various attributes are 
treated as columns of a matrix. The raw data corresponding to 
the intersection of each time period and attributes are treated as 
the entries of a matrix. Thus, from the raw data, we obtain the 
initial M x N matrix. 

The number of individual hours in each time period might 
not be the same. Hence, in the second stage, in order to obtain 
consistency in the initial matrix we use the “hourly concept”. 
The entries corresponding to the intersection of each of the time 
periods and the attributes of the initial matrix are transformed so 
that each new entry corresponds to the hourly rate viz. total 
passengers per hour, total collection per hour, total number of 
trips per hour, total occupancy per hour etc. 

Thus, we convert the above initial matrix into the Average 
Time Dependent Data matrix (ATD matrix) i.e., [a^MxN. 

In the third stage, we use the average and the standard 
deviation to convert to above average time dependent data 
matrix into a matrix with entries ey, c,, e {-1, 0, 1}, where i 
represents the i th row and j represents the j lh column. We call 
this newly formed matrix as the Refined Time Dependent Data 
matrix (RTD matrix) i.e., [eij] MxN . 

The value of the entry e,j corresponding to each intersection 
is determined from an interval. This interval is obtained strictly 
by using the average and the standard deviation calculated from 
the raw data. The choice of the interval made by us might not 
be able to identify the accurate peak hour of the route; hence, 
we introduce and define a parameter a, which enables us to get 
a better solution. 

We calculate the mean p, and the standard deviation a, for 

each attribute j, j = 1, 2 N using the data given in the 

average time dependent data matrix. For varying values of the 
parameter a, where a e [0, 1]. We follow a rule and determine 
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the value of the entry e,j in the refined time dependent data 
matrix. 

For each of the attributes j (j = 1, 2, . N) we have the rule: 

If ay < (pj - a * Gj) then eij = -1 ; 

else, 

if a;j g (pj - a * Gj, pj + a * Gj) then e^ = 0; 

else 

if ay > (pj + a * Gj) then e^ = 1 ; 

here denotes the usual multiplication. 

Thus, for different values of a, we obtain different refined 
time dependent data matrices. The main putpose of introducing 
the refined time dependent data matrix is only to minimize the 
time involved in performing the simple arithmetic calculations 
and operations on the matrix. 

In the fourth stage, we bring in the notion of the Combined 
Effect Time Dependent Data matrix (CETD matrix) i.e., [Cij] Mx N, 
which gives the combined effect of all the refined time 
dependent data matrices obtained by varying the parameter a. 
In the final stage, we add up the rows of the combined effect 
time-dependent data matrix. The overall time period utility of a 
route is obtained by inferring the row sums of a combined effect 
time-dependent data matrix. The highest positive value is taken 
as the highly utilized time period of a route and next lower 
value is taken as the next peak hour for the same route. Thus, 
for a particular route, we grade the utilization rate of the 
different time periods. 

The computation starts by computing the ATD matrix from 
the initial M x N matrix. Then mean and standard deviations 
are computed for the ATD vector. Based on mean, standard 
deviation and ATD vector, the RTD and CETD matrices are 
computed either physically or using a programme in Java or 
C++. 
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It is pertinent to mention here this technique of analysing 
raw data can be used in case of any problem which is time 
dependent. This is a very new way of applying matrix theory to 
find solution for a collected raw data. 

Parameters play a vital role in estimating the peak hours of 
a route. This is clearly illustrated for the route 18JJ and the 
route 17D by varying the parameter a. First, consider the route 
18JJ. We have the following: 

The initial M x N matrix obtained from the raw data is as 
below: 



81.00 63.00 

925.00 157.00 

1040.00 181.00 

1159.00 271.00 

1596.00 400.00 

1673.00 451.00 

322.00 160.00 

The corresponding average I 



3.00 


147.20 


16.00 


1394.55 


15.00 


1585.50 


13.00 


1916.05 


18.00 


2615.30 


16.00 


2889.60 


6.00 


522.50 



dependent data matrix is 



40.50 


31.50 


1.50 


73.60 


308.33 


52.33 


5.33 


464.85 


346.67 


60.33 


5.00 


528.50 


386.33 


90.33 


4.33 


638.68 


532.00 


133.33 


6.00 


871.77 


557.67 


150.33 


5.33 


963.20 


161.00 


80.00 


3.00 


261.25 
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The refined time dependent data matrix corresponding to a 
single value of a = 0. 1 is 



-1 -1 -1 -1 
-1 -1 1 -1 
0-110 
110 1 
1111 
1111 
-1 -1 -1 -1 



The corresponding row sum of the above refined time 
dependent data matrix is given by 

Row sum for the first row is -4 
Row sum for the second row is -2 
Row sum for the third row is 0 
Row sum for the fourth row is 3 
Row sum for the fifth row is 4 
Row sum for the sixth row is 4 
Row sum for the seventh row is -4. 



From the row sums of the above refined time dependent 
data matrix, we observe that time period corresponding to the 
fifth row i.e., time period {Fl 16 , Fin, H 18 } and the sixth row i.e., 
time period {Fl^, FFo, FFi} are the peak hours of the route 
followed by the fourth row viz. time period {Fin, H ]4 , FF 5 }. 
Thus, graphically, we have this as below: 




Average Time Dependent (ATD) Data Matrix 
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Graph depicting the peak periods of the 
Route 18JJ for a = 0.1 




The refined time dependent data matrix corresponding to 
the value of a = 0.4 follows: 

“-1 -1 -1 -f 

0-110 
0-110 
0 0 0 0 
1111 
1111 
-1 0 -1 -1 

From the above, we have the row sums as given below: 

Row sum for the first row is -4 
Row sum for the second row is 0 
Row sum for the third row is 0 
Row sum for the fourth row is 0 
Row sum for the fifth row is 4 
Row sum for the sixth row is 4 
Row sum for the seventh row is -3. 

From the row sums of the above refined time dependent 
data matrix, we observe that time period corresponding to the 
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fifth row i.e., time period {Hi 6 , Hn, H i8 } and the sixth row i.e., 
time period {H 19 , H 20 , H 2 i} are the peak hours of a day. Time 
periods {H 7 , H 8 , H 9 }, {H 10 , H n , H 12 } and {H 13 , H 14 , H 15 } have a 
row sum of zero indicating that these time periods are neither 
the peak hours of a day nor the non-peak hours of a day. Hence, 
the number of trips made in these time periods should be kept 
unaltered for the concern to maintain the profit. 

Thus, graphically, we have this as below: 



Graph depicting the peak hour for the 
Route 18JJ for a = 0.4 




The refined time dependent data matrix corresponding to 
the value of a = 0.7 follows: 



-1 


-1 


-1 


-1 


0 


-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


1 


1 


1 


1 


1 


0 


1 


-1 


0 


-1 


-1 
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From the above, we have the row sums as given below: 

Row sum for the first row is -4 
Row sum for the second row is -1 
Row sum for the third row is 0 
Row sum for the fourth row is 0 
Row sum for the fifth row is 4 
Row sum for the sixth row is 3 
Row sum for the seventh row is -3. 

From the row sums of the above refined time dependent 
data matrix, we observe that time period corresponding to the 
fifth row i.e., time period {FF 6 , Hn, Hi 8 } is the first peak hour 
of the day followed by the sixth row i.e., time period {Hi 9 , H 20 , 
H 2 i}. Time periods {H 7 , H 8 , H 9 }, {H 10 , H n , H i2 } and {H u , H 14 , 
H i5 } have a row sum of zero indicating that these time periods 
are neither the peak hours of a day nor the non-peak hours of a 
day. Thus, graphically, we have this as below: 



Graph depicting the peak hour for the route 
18 JJ for a = 0.7 




Hours of the day 
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The refined time dependent data matrix corresponding to 
the value of a = 0.9 follows: 

"-1 -1 -1 -f 
0 0 0 0 
0 0 0 0 
0 0 0 0 
1111 
110 1 
-1 0 -1 -1 



From the above, we have the row sums as given below: 

Row sum for the first row is -4 
Row sum for the second row is 0 
Row sum for the third row is 0 
Row sum for the fourth row is 0 

Row sum for the fifth row is 4 
Row sum for the sixth row is 3 
Row sum for the seventh row is -3. 



The row sums of the above matrix indicates time periods 
corresponding to the fifth row i.e., { Hi 6 , Fin, H i H } is the first 
peak hour of the day followed by the sixth row i.e., { H 19 , FF 0 , 
H 2i }. Time periods {H 7 , H 8 , H 9 }, {H 10 , H n , H 12 } and {H 13 , H 14 , 
H 15 } have a row sum of zero indicating that these time periods 
are neither the peak hours of a day nor the non-peak hours of a 
day. Thus, graphically, we have this as below: 
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Graph depicting the peak hour for the 
route 18 JJ for a = 0.9 



From the above analysis, we observe that the peak hours of 
a route vary from one time period to another with the change in 
the value of the parameter from 0 to 1. The row sums matrix 
obtained in the above cases were specific to only one value of 
the parameter. But the combined effect time dependent matrix 
for all the values of a s [0, 1] is given below: 

The combined time dependent data matrix for all a e [0, 1] 
is 



-10 -10 -10 -10 
-1 -8 6 -2 
0-640 
3 10 3 

10 10 10 10 
10 10 6 10 
-9 -1 -9 -9 



The row sum of the above combined time dependent data 
matrix is 




Hours of the day 
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Row = 1 , row sum - -40 
Row = 2, row sum - -5 
Row = 3, row sum - -2 
Row = 4, row sum = 7 
Row = 5, row sum = 40 
Row = 6, row sum = 36 
Row = 7, row sum = -28. 

From the row sums of the above matrix, we observe that the 
time period corresponding to the fifth row i.e., time period { H i 6 , 
Hn, Hi 8 } is the first peak hours of the day followed by the sixth 
row i.e., time period {H 19 , H 2 o, FFi} and the fourth row 
corresponding to the time period { H i3 , Hu, Hi 5 }. All the other 
time periods have a row sum of negative value indicating that 
these time periods are the non-peak hours of a day. The first 
time period has the least negative value indicating that services 
can be curtailed to the greatest extent in this time period. 

Thus, graphically, we have the peak hour as depicted below: 
From the graph below, we observe that the peak hours as 
obtained from the combined time dependent data matrix, given 
by the cumulative effect of all the values of a s [0, 1], gives the 
true peak hours of a day. These results coincide well with the 
observed hours of the route 18JJ. 



Graph depicting the peak hour for the 
Route 18 JJ for a e [0, 1] 




50 

40 

30 



40 

50 



Hours of the day 
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We give a grading of the peak hours since they are essential 
and helpful to the transportation sector in a number of ways: 

1. It is efficient in identifying the highly utilized time 
periods along each route and enables the transport 
sector to operate more number of buses during these 
hours. 

2. The transport sector identifies the poorly utilized time 
periods and makes necessary arrangements to overcome 
the loss incurred due to the services scheduled in these 
time periods. They might either discard all the trips 
made during these hours or operate minimum number 
of services. 

3. The row sum of the combined time dependent data 
matrix of each route helps in identifying the overall 
utility rate of each route. In other words, if all the row 
sums of the combined time dependent data matrix are 
positive, then we conclude that the particular' route is 
highly utilized throughout the day. 

4. A positive value of the row sum of the combined time 
dependent data matrix indicates the maximum 
utilization of that specific time period. 

5. A negative value of the row sums of the combined time 
dependent data matrix indicates that the passenger 
patronage and the collection in those time periods are 
very less. Hence, the transport sector can curtail the 
number of services totally, partially or operate mini 
buses to overcome the loss. 

6. A zero value of the row sum of the combined time 
dependent data matrix indicates that the corresponding 
time periods are neither the peak hours of a day nor the 
non-peak hours of a day. 
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It is to be noted here that as the total number of time periods 
(M) increases, the accuracy in identifying the utilization rate of 
the different time -periods also increases. We substantiate our 
claim by taking two different values for the total number of time 
periods and this is carried out in the illustration I. 

Illustration I (Route 1 8R) 

The number of time periods is five (M = 5). 

Here, we take five time periods viz. early morning Ci:{H 6 , 
H 7 }; morning C 2 :{H 8 , H 9 , H 10 , H n }; noon C 3 : {Hi 2 , H 13 , H 14 , 
H 15 }; evening C 4 :{H 16 , H 17 , H, 8 , Hi 9 }; and night C 5 :{H 20 , H 2 i, 
H 22 }. We obtain the row sum of the combined time dependent 
data matrix as follows: 

The combined time dependent data matrix corresponding to 
M = 5 is given below: 



-10 -8 -10 -10 
6-265 
2-16 1 . 

10 10 6 10 
-4 -6 0 -3 

The row sums of the corresponding combined time 
dependent data matrices is given by: 

Row sum for the time period early morning = -38 
Row sum for the time period morning =15 
Row sum for the time period noon = 8 
Row sum for the time period evening = 36 
Row sum for the time period night = -13 

The number of time periods is seventeen (M = 17). 

Here, we take 17 time periods viz. early morning Cj : { H 6 } ; 
morning C 2 : {H 7 }; C 3 : {H 8 }, C 4 : {H 9 }; early noon C 5 : {H 10 }; 
Cg: { Hi i } , C 7 : { H 42 } ; noon C 8 : {Hi 3 }; C 9 : {Hi 4 }, Cio: { Hig } ; 




Average Time Dependent (ATD) Data Matrix 



25 



evening C n ;{H 16 }; C l2 : {H 17 }, C n : {H 18 }; late evening C, 4 : 
{H 19 }, Ci 5 :{H 20 }, Ci 6 :{H 2 i}; and night C n : { H 22 } . The row 
sums of the combined time dependent data matrices are 
obtained as follows: 

The combined time dependent data matrix corresponding to 
M = 17 is as follows: 



-10 


-2 


-10 


-10 


-10 


-7 


0 


-10 


0 


-5 


0 


-5 


4 


-2 


0 


2 


7 


-1 


0 


5 


0 


-2 


0 


0 


1 


-3 


10 


2 


3 


0 


0 


2 


2 


-1 


0 


0 


-10 


-2 


-10 


-10 


10 


1 


10 


10 


6 


-1 


0 


6 


10 


3 


0 


10 


-4 


10 


-10 


-4 


10 


0 


10 


10 


-4 


-2 


0 


-3 


-10 


-7 


0 


-10 



The corresponding row sum of the above combined time 
dependent data matrix is given as below: 

Row sum for the time period C i : { H 6 } = -32 
Row sum for the time period C 2 : {H 7 } = -27 
Row sum for the time period C 3 : { H 8 } =-10 
Row sum for the time period C 4 : { H 9 } =4 
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Row sum for the time period C5 :{Hio} = 11 
Row sum for the time period C 6 : {Hu } = -2 
Row sum for the time period C 7 : { H 12 } = 10 
Row sum for the time period C x : { H , ? } = 5 

Row sum for the time period Cg: { 11 !4 } = 1 
Row sum for the time period Cio: {H 15 } =-32 
Row sum for the time period Cn: { Hi 6 } = 31 
Row sum for the time period Ci 2 : {H 17 } = n 

Row sum for the time period C13: {Hi S } = 23 
Row sum for the time period Ci 4 : {H 19} - -8 
Row sum for the time period C15: {H 20 } = 30 
Row sum for the time period Ci6 : {H 2 i } = -9 
Row sum for the time period C17: { H 22 } = -27. 

From the above Illustration I. on the analysis of the row 
sums of the combined time dependent data matrices of M = 5, 
we notice that the time period corresponding to row = 4 namely 
the evening hours {Hi 6 , Hn, Hig, H19} are the busy hours 
followed by row = 2 namely the morning hours {H 8 , H y . Hi 0 , 
Hn}. 



On analyzing the row sums of the combined time dependent 
data matrices of M = 17, we observe that row =11 viz. evening 
{ H 16 } corresponds to the busiest hour of the route followed by 
row = 15 viz. late evening { H 20 } , row = 13 viz. evening { Hi S } - 

Thus, we observe that as the number of time periods 
increases, the accuracy in identifying the busier hours of a route 
also increases. Also, we can observe the added effects of one or 
more attributes on the various time periods. 

Thus, for a e [0, 1] with a step function of a = 0.1, the 
peak hours of the route are Hi 6 , H 20 and H !X . Time period 
corresponding to H 6 is poorly utilized and hence can be 
discarded. 
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The same problem is analyzed and the combined time 
dependent data matrix for varying values of a e [0, 1] with a 
step function of a = 0.05 is given below: 

The combined time dependent data matrix, the row sum and 
the corresponding graph is given as below: 

The combined time dependent data matrix 



-10 


-4 


-10 


-10 


-10 


-10 


0 


-10 


-1 


-10 


0 


-10 


8 


-4 


0 


5 


10 


-3 


0 


10 


-1 


-5 


0 


1 


2 


-7 


10 


5 


7 


-1 


0 


4 


4 


-3 


0 


1 


-10 


-4 


-10 


-10 


10 


2 


10 


10 


10 


-2 


0 


10 


10 


7 


0 


10 


-8 


10 


-10 


-9 


10 


-1 


10 


10 


-9 


-5 


0 


-7 


-10 


-10 


0 


-10 
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The row sums of the above combined time dependent data 
matrix is 

Row = 1 , row sum = -34 
Row = 2, row sum - -30 
Row = 3, row sum = -21 
Row = 4, row sum = 9 

Row = 5, row sum = 17 
Row = 6, row sum = -5 
Row = 7, row sum = 10 
Row = 8, row sum = 10 

Row = 9, row sum = 2 
Row =10, row sum = -34 
Row =11, row sum = 32 
Row = 12, row sum =18 

Row =13, row sum = 27 
Row =14, row sum = -17 
Row =15, row sum = 29 
Row =16, row sum = -21 
Row = 17, row sum = -30. 

We illustrate and describe the problem using 18R service. 

We have applied our algorithm to a single route and 
obtained a comparison of the observed and the calculated data. 
We have taken the route 18R plying from Parrys to Dharmaraja 
Koil and the data was obtained from the Pallavan Transport 
Corporation Limited, Madras, India. 

Here, we take 17 time periods viz. early morning Ci:{H 6 }; 
morning C 2 : {H 7 }; C 3 :{H 8 }, C 4 : {H 9 }; early noon C 5 : {Hi 0 };C 6 
: { Hi i } . C 7 : { H12 } ; noon C 8 : {Hn}; C 9 :{Hj 4 }, Cio : { H15 } ; 
evening C n : {H 16 }; C i2 : {H 17 }, C i3 : {H 18 }; late evening C u : 
{H 19 }, C15 : {H 20 }, Cie : {H 2i }; and night C n : {H 22 } and also 
four attributes viz. A] : total population; A 2 : hourly occupancy; 
A 3 : number of hips and A 4 : total collection. 
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We view this problem from the aspect of both the operators 
as well as the users. 

The initial M x N Matrix is as follows: 



96.00 


64.00 


1.00 


156.95 


71.00 


21.00 


2.00 


103.55 


222.00 


39.00 


2.00 


244.90 


269.00 


63.00 


2.00 


353.50 


300.00 


70.00 


2.00 


392.80 


220.00 


60.00 


2.00 


328.30 


241.00 


53.00 


3.00 


353.65 


265.00 


78.00 


2.00 


348.65 


249.00 


70.00 


2.00 


327.85 


114.00 


63.00 


1.00 


151.60 


381.00 


91.00 


3.00 


526.40 


288.00 


73.00 


2.00 


407.05 


356.00 


112.00 


2.00 


515.15 


189.00 


389.00 


1.00 


252.75 


376.00 


78.00 


3.00 


569.40 


182.00 


59.00 


2.00 


261.90 


67.00 


20.00 


2.00 


81.35 
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The average time dependent data matrix is as follows: 



96.00 


64.00 


1.00 


156.95 


71.00 


21.00 


2.00 


103.55 


222.00 


39.00 


2.00 


244.90 


269.00 


63.00 


2.00 


353.50 


300.00 


70.00 


2.00 


392.80 


220.00 


60.00 


2.00 


328.30 


241.00 


53.00 


3.00 


353.65 


265.00 


78.00 


2.00 


348.65 


249.00 


70.00 


2.00 


327.85 


114.00 


63.00 


1.00 


151.60 


381.00 


91.00 


3.00 


526.40 


288.00 


73.00 


2.00 


407.05 


356.00 


112.00 


2.00 


515.15 


189.00 


389.00 


1.00 


252.75 


376.00 


78.00 


3.00 


569.40 


182.00 


59.00 


2.00 


261.90 


67.00 


20.00 


2.00 


81.35 



The refined time dependent data matrices corresponding to 
different values of alpha are calculated and from them we obtain 
the combined time dependent data matrix. 



Since, the refined time dependent data matrices are obtained 
from the average time dependent data matrices they have a size 
equivalent to that of the average time dependent data matrices 
but, they have entries only in {-1, 0, 1}. For different values of 
the parameter we have different refined time dependent data 
matrices. 
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The refined time dependent data matrix corresponding to 
the value of a = 0. 1 is 



-1 


-1 


-1 


-1 


-1 


-1 


0 


-1 


0 


-1 


0 


-1 


1 


-1 


0 


1 


1 


-1 


0 


1 


0 


-1 


0 


0 


1 


-1 


1 


1 


1 


0 


0 


1 


1 


-1 


0 


0 


-1 


-1 


-1 


-1 


1 


1 


1 


1 


1 


-1 


0 


1 


1 


1 


0 


1 


-1 


1 


-1 


-1 


1 


0 


1 


1 


-1 


-1 


0 


-1 


-1 


-1 


0 


-1 



The refined time dependent data matrix corresponding to the 
value of a = 0.2 is 



-1 


-1 


-1 


-1 


-1 


-1 


0 


-1 


0 


-1 


0 


-1 


1 


-1 


0 


1 


1 


0 


0 


1 


0 


-1 


0 


0 


0 


-1 


1 


1 


1 


0 


0 


1 


1 


0 


0 


0 


-1 


-1 


-1 


-1 
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1 


0 


1 


1 


1 


0 


0 


1 


1 


1 


0 


1 


-1 


1 


-1 


-1 


1 


0 


1 


1 


-1 


-1 


0 


-1 


-1 


-1 


0 


-1 



The refined time dependent data matrix corresponding to the 
value of a = 0.3 is 

-1 0 -1 -1 
-1 -1 0 -1 
0 - 10-1 
10 0 0 
10 0 1 
0 0 0 0 
0-110 
10 0 0 
0 0 0 0 
-1 0 -1 -1 
10 11 
10 0 1 
110 1 
-1 1 -1 -1 
10 11 
-10 0-1 
-1 -1 0 -1 
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The refined time dependent data matrix corresponding to the 
value a = 0.4 is 



-1 0 
-1 -1 
0 -1 
1 0 
1 0 
0 0 
0 0 
0 0 
0 0 
-1 0 
1 0 
1 0 
1 0 
-1 1 
1 0 
-1 0 
-1 -1 

The refined time dependent d 
value a = 0.5 is 



-1 -1 
0 -1 
0 -1 
0 0 
0 1 
0 0 
1 0 
0 0 
0 0 
-1 -1 
1 1 
0 1 
0 1 
-1 -1 
1 1 
0 0 
0 -1 

a matrix corresponding to the 



-1 0 -1 -1 

-1 -1 0 -1 

0-10-1 
0 0 0 0 

10 0 1 
0 0 0 0 

0 0 10 
0 0 0 0 

0 0 0 0 

-1 0 -1 -1 
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10 11 

10 0 1 

10 0 1 

0 1-10 
10 11 

0 0 0 0 

-1 -1 0 -1 



The refined time dependent data matrix corresponding to the 
value of a = 0.6 is 

-1 0 -1 -1 

-1 -1 0 -1 

0 0 0 0 

0 0 0 0 

10 0 0 

0 0 0 0 

0 0 10 

0 0 0 0 

0 0 0 0 

-1 0 -1 -1 

10 11 

10 0 1 

10 0 1 

0 1-10 
10 11 

0 0 0 0 

-1 -1 0 -1 
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The refined time dependent data matrix corresponding to the 
value of a = 0.7 is 



-1 0 
-1 -1 
0 0 
0 0 
1 0 
0 0 
0 0 
0 0 
0 0 
-1 0 
1 0 
0 0 
1 0 
0 1 
1 0 
0 0 
-1 -1 

The refined time dependent d 
value of a = 0.8 is 



-1 -1 
0 -1 
0 0 
0 0 
0 0 
0 0 
1 0 
0 0 
0 0 
-1 -1 
1 1 
0 0 
0 1 
-1 0 
1 1 
0 0 
0 -1 

a matrix corresponding to the 



-1 0 -1 -1 

-10 0-1 
0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 10 
0 0 0 0 

0 0 0 0 

-1 0 -1 -1 




36 



Innovative Uses of Matrices 



10 11 

0 0 0 0 

10 0 1 
0 1-10 
10 11 

0 0 0 0 

-10 0-1 



The refined time dependent data matrix corresponding to the 
value of a = 0.9 is 

-1 0 -1 -1 

-10 0-1 
0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 10 
0 0 0 0 

0 0 0 0 

-1 0 -1 -1 

10 11 
0 0 0 0 

10 0 1 
0 1-10 
10 11 
0 0 0 0 

-10 0-1 
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The refined time dependent data matrix corresponding to the 
value of a = 1 .00 is 



-1 0 -1 -1 

-10 0-1 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 10 
0 0 0 0 

0 0 0 0 

-1 0 -1 -1 

10 11 

0 0 0 0 

10 0 1 
0 1-10 
10 11 

0 0 0 0 

-10 0-1 




38 Innovative Uses of Matrices 



The combined time dependent data matrix 



-10 


-2 


-10 


-10 


-10 


-7 


0 


-10 


0 


-5 


0 


-5 


4 


-2 


0 


2 


7 


-1 


0 


5 


0 


-2 


0 


0 


1 


-3 


10 


2 


3 


0 


0 


2 


2 


-1 


0 


0 


-10 


-2 


-10 


-10 


10 


1 


10 


10 


6 


-1 


0 


6 


10 


3 


0 


10 


-4 


10 


-10 


-4 


10 


0 


10 


10 


-4 


-2 


0 


-3 


-10 


-7 


0 


-10 



The row sums corresponding to the above combined time 
dependent data matrix are: 

Row sum for the time period Ci:{H 6 } = -32 
Row sum for the time period C 2 :{H 7 } = -27 
Row sum for the time period C3:{H 8 } = -10 
Row sum for the time period C 4 :{H<, } = 4 
Row sum for the time period C5: { Hi 0 } = 11 

Row sum for the time period Cg:{Hn } = -2 
Row sum for the time period C 7 :{Hi 2 } = 10 
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Row sum for the time period C 8 : { H j = 5 
Row sum for the time period C^jH^} = 1 



Row sum for the time period Cio:{Hi 5 } = -32 
Row sum for the time period Cn:{Hi 6 } = 31 

Row sum for the time period Ci2:{Hi 7 } = 11 
Row sum for the time period C13: {Hi S } = 23 
Row sum for the time period QpjHig} = -8 
Row sum for the time period Ci5:{H 2 o} = 30 
Row sum for the time period Ci6:{H 2 i} = -9 
Row sum for the time period Ci 7 :{H 22 } = -27 



We observe that peak hours of a route differ for varying 
values of alpha. Hence, we form the combined effect time 
dependent data matrix which gives the combined effect of all 
the values of alpha in the interval [0, 1]. Thus, the highest value 
of the row for the corresponding time period identified from the 
row sums of the combined time dependent data matrix is 
observed as the best peak period of the particular route. The 
next value of the row sum in the small matrix is observed to be 
the next peak period of the same route. In this manner, we rank 
the different time periods of a route. 

In the above estimated method, we notice that the time- 
period corresponding to row = 11 viz. evening { Hj 6 } 
corresponds to the busiest hour of the route followed by row = 
15 viz. late evening {H 20 }, and row = 13 viz. evening {H 18 }. 
Thus, we rank the peak periods for the route 18R to be as 
evening Cn:{Hi 6 }; late evening C15: {H 20 } and evening C13 : 
{H 18 } in the order of first, second, third respectively. 

We now compare the predictions with the ridership method. 

To establish the efficiency of our method, we compare our 
results with the ridership, approach followed by Dr. J. 
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Vasudevan in his research on the rationalization of city bus 
routes. 

He had taken Madras as a case study where, he had 
discussed about the rationalization of the routes and had studied 
for a total of 90 routes. A small section of it deals with the 
average hourly bus passenger load for the city. For each hour 
ending, say 5 a.m., 6 a.m., ..., 21 p.m., 22 p.m., he had 
calculated the total number of passengers who had traveled 
during that particular hour. 

He had also calculated the percentage of the passenger load 
at any hour against the total strength of passengers who had 
traveled that day. This is clearly depicted by the following table. 



Determination of the peak hour based on the 
Passenger load per hour 



Hour ending 


Passengers / hour 


Percent to total 


6 


96 


2.47 


7 


71 


1.83 


8 


222 


5.71 


9 


269 


6.92 


10 


300 


7.72 


11 


220 


5.66 


12 


241 


6.20 


13 


265 


6.82 


14 


249 


6.41 


15 


114 


2.93 


16 


381 




17 


288 


7.41 


18 


356 


9.16 


19 


189 


4.86 


20 


376 


9.68 


21 


182 


4.68 


22 


67 


1.72 


Total 


3886 


100 






Average Time Dependent (ATD) Data Matrix | 41 



We have applied the matrix method to this data in order to 
obtain a comparison and establish that our method is better than 
the ridership method. From the last column of the above table, 
we identify the first peak period as the row corresponding to the 
greatest percent to the total, i.e., 9.80 corresponding to the hour 
ending 16, i.e., evening 4 p.m., followed by 9.68 corresponding 
to the hour ending 20 i.e., night 8 p.m. 

The third peak period is the percent 9.16 corresponding to 
the hour ending 18 i.e., evening 6 p.m. In this fashion, we rank 
the peak periods of the route based on the passenger load. 

Thus, we have established that our matrix method gives us a 
better result that the ridership method of Vasudevan. 



Comparison Table depicting the peak hours 
as estimated by the three methods 



Data as 
obtained 
from the 
transport 
corporation 


Identification of the peak and non-peak 
periods from the ridership method, the 
revenueship method and comparison of the 
two methods with the matrix method 


Hour 

ending 


Observed 
Method I: 

Passengers 

/hour 


Observed 
Method II: 
Revenue/hour 


Matrix 
Method: 
row sum of 
the CETD 
Matrix 


6 


96 


156.95 


-32 


7 


71 


103.55 


-27 


8 


222 


244.9 


-10 


9 


269 


353.5 


4 


10 


300 


392.8 


11 


11 


220 


328.3 


-2 


12 


241 


353.65 


10 
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13 


265 


348.65 


5 


14 


249 


327.85 


1 


15 


114 


151.6 


-32 


16 


381(1)* 


526.4 (II)* 


31 (I)* 


17 


288 


407.05 


11 


18 


356 (III)* 


515.15 (III)* 


23 (III)* 


19 


189 


252.75 


-8 


20 


376 (II)* 


569.4 (I)* 


30 (II)* 


21 


182 


261.9 


-9 


22 


67 


81.35 


-27 



*denotes the ranking of the time periods as given within parenthesis 



From the above table, we obtain a comparative study of the 
observed methods and the estimated method. From the last 
column of the above table, we observe that there are rows that 
have a negative value. 

These negative values indicate that the corresponding time 
periods are comparatively poor in performance i.e., in the 
passenger usage and in the collection. Flence, the services in 
these time periods can be discarded or the frequency of the 
services can be reduced. 

This is one of the innovative means of applying matrices in 
analysis problems which has only raw data. Thus we have not 
given the graphs of them. 





Chapter Three 



Fuzzy Jjngusiic Mvt 
Applications 



N iK 



AND THEIR 



In this chapter we introduce the new notion of fuzzy 
linguistic matrices and use them in fuzzy linguistic models and 
apply them to social problems. We need for this the operations 
on the fuzzy linguistic matrices. We briefly describe the 
operations on them. Let M be any matrix if the entries of M are 
taken from the fuzzy linguistic set we call such matrix M to be a 
fuzzy linguistic matrix for more refer [ 53 ]. We describe 
operations on them. Let L denote the collection of fuzzy 
linguistic terms. 



Let x = (ai, a2, a3, a.4, as) be a fuzzy linguistic row matrix 
(vector) with entries from L. 



a i 



x' — (ai, a2, a 3 , a.4, as ) 1 — 



a 2 

a 3 

a 4 

a 5 



is the fuzzy linguistic column matrix (vector). 
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That is if x = (0, bad, good, fair, very bad) be the fuzzy 
linguistic row matrix, then transpose of x is denoted by 

0 

bad 
good 
fail- 

very bad 

which is a fuzzy linguistic column matrix. Suppose we have 
product of a fuzzy linguistic row matrix x; it can be operated 
with min or max operation with a fuzzy linguistic column 
matrix y if and only if x is a 1 x t row matrix then y must be a 
t x 1 column matrix. 

We show how the ‘min’ operation looks like. 

Suppose 



x‘ = (0, bad, good, fair, very bad) 1 = 



y = 



and x = (hi, b 2 , ...,b t ) 



then min {min {x, y } } = min {min {a b t^}, min {a 2 , b 2 }, ..., 
min {a„ b,}}. 

Likewise we can have {max; max} operation; {max; min} 
operation and {min; max} operation. All these will be 
illustrated by examples. 



Let x = (x b x 2 , ..., x 20 ) and y = (y lf y 2 , ..., y 20 ) with x : , y ; e 
L; 1 < i < 20. L be a fuzzy linguistic matrix. 

max {max (x, y)} = max {max (xi, y i ), max {x 2 , y 2 }, ..., 
max { x 2 o, y 20 }}. 
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Let x = (x ls x 2 , ... x 15 ) and y = (y l5 y 2 , ..., y i5 ) 
max (min {x, y}) = max {min {xi, y i } , min {x 2 , y 2 }, .... min 

{xis, yi 5 } }. 

Likewise min (max {x, y}) = min {max {xi, y i } , 
max { x 2 , y 2 }, ..., max {xi 5 , yi 5 } }- 



We will see how for a typical problem the four operations 
give different sets of answers. 

Let 



good 

bad 

fair 



x = 



0 

best 

bad 



and 



0 



good 



y = (bad. 0, good, very bad, best, fair, best, better) 
be two fuzzy linguistic matrices, 
min {min (y, x)} 

= min {bad, 0, fair, 0, best, bad, 0, better} = 0 
max { min (y, x) } 

= max {bad, 0, fair, 0, best, bad, 0, better} = best. 
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min { max (y, z) } 

= min { good, bad, good, very bad, best, fair, 
best, good} 

= very bad and 

max (max {y, x}) 

= max {good, bad, good, very bad, best, fair, best, good} 

= best. 

Suppose we are interested in finding min {x, y} then we 
have 

min {x, y} = 



bad 


0 


good 


very bad 


good 


fair 


good 


better 


bad 


0 


bad 


very bad 


bad 


bad 


bad 


bad 


bad 


0 


fair 


very bad 


fair 


fair 


fair 


fair 


0 


0 


0 


0 


0 


0 


0 


0 


bad 


0 


good 


very bad 


best 


fail- 


best 


better 


bad 


0 


bad 


very bad 


bad 


bad 


bad 


bad 


0 


0 


0 


0 


0 


0 


0 


0 


bad 


0 


good 


very bad 


good 


fair 


good 


better 



This is the way a fuzzy linguistic column matrix of order 
8 x 1 is multipliced with a 1 x 8 fuzzy linguistic row matrix 
with min operation gives a 8 x 8 fuzzy linguistic matrix. 
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We can also find max { x, y } = 



good 


good 


good 


good 


best 


good 


good 


good 


bad 


bad 


fair 


bad 


best 


bad 


bad 


better 


fail- 


fair 


fair 


fair 


best 


fail- 


fair 


better 


bad 


0 


fair 


0 


best 


bad 


0 


better 


best 


best 


best 


best 


best 


best 


best 


best 


bad 


bad 


fair 


bad 


best 


bad 


bad 


better 


bad 


0 


fair 


0 


best 


bad 


0 


better 


good 


good 


good 


good 


best 


good 


good 


good 



We see clearly max {x, y} ^ min {x, y}. 

Further 0 dominates in min and best dominates in max. 
When an expert wants to boost the results on the positive side 
he can use the max operation. If min operation is used it gives 
the worst state of affairs. 

Thus we see we get different results for these four types of 
operations. According to need one can use any one of the 
operations. 

Now we find the transpose of a 6 x 5 fuzzy linguistic matrix 
M. 



good 


best 


bad 


fair 


best 


0 


bad 


worst 


0 


very bad 


bad 


fair 


very fair 


bad 


better 


worst 


0 


good 


better 


0 


0 


fair 


0 


worst 


best 


good 


fair 


best 


worst 


0 



Now the transpose of this fuzzy linguistic 6x5 matrix M 
denoted 
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good 


0 


bad 


worst 


0 


good 


best 


bad 


fair- 


0 


fair 


fair 


bad 


worst 


very fair 


good 


0 


best 


fair 


0 


bad 


better 


worst 


worst 


best 


very bad 


better 


0 


best 


0 



Likewise if 



best 


0 


bad 


fair 


bad 


0 


0 


very fair 


good 


good 


0 


best 


best 


worst 


0 


0 


best 


good 


good 


better 


best 



is any fuzzy linguistic 7x3 matrix, to find N T . 



N t = 



best fair 0 good best 0 good 

0 bad very fair 0 worst best better . 

bad 0 good best 0 good best 



We see N T is a 3 x 7 fuzzy linguistic matrix. 

Now we can find the product of two rectangular linguistic 
matrices M and N if M is a n x t matrix then N must be a t x m 
matrix then only MN is defined how ever NM is not defined in 
this case and MN is a n x nr fuzzy linguistic matrix. 

However product of any n x n matrix with itself is always 
defined. 



We will illustrate these two situations by some examples. 
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Let 



M = 



good 


bad 


0 


worst 


bad 


0 


good 


best 


0 


good 


bad 


0 


good 


bad 


good 


bad 




ro 


0 0 


ol 



= min (min (M, M)) = 



0 0 0 0 
0 0 0 0 
0 0 0 0 



We see if ‘O' occurs atleast once in every row and atleast 
once in every column then min (min (M, M)) = (0). 

Now we find 



max (min (M, M)) = 



We now find the value of 



good 


bad 


bad 


bad 


good 


good 


good 


bad 


bad 


bad 


good 


good 


good 


good 


good 


good 



good 


good 


good 


best 


best 


best 


best 


best 


good 


good 


good 


best 


good 


good 


good 


best 



max (max (M, M)) = 
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min (max (M, M) = 



0 


bad 


bad 


0 


bad 


0 


bad 


bad 


bad 


bad 


0 


bad 


bad 


bad 


bad 


bad 



We see max {max (M, M)} gives an extreme or better 
values and min {min {M, M}} gives an extreme low values, 
where as min max non negative values and max min more 
positive values. 

Such four types of operations can be used as per need of the 
problem. 

Now if M is a n x m matrix and N is a m x t matrix we can 
find min (min {M, N}), min (max {M, N}), max {max (M, N)} 
and max (min {M, N}). 

We will illustrate these four types of operations. 

Let 



low 


0 


high 


0 


very low 


0 


high 


0 


low 


high 


high 


0 


medium 


0 


low 


low 


low 


0 


medium 


0 





low 


0 


medium 


0 




high 


low 


0 


high 


= 


0 


very high 


low 


0 




medium 


0 


high 


low 




very low 


medium 


0 


high 
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be two fuzzy linguistic matrices associated with temperature of 
an experiment. 



Clearly min (min {M, N}) = 



0 

0 

0 

0 



0 

0 

0 

0 



0 

0 

0 

0 



0 

0 

0 

0 



max {min (M. N)} = 



low 


high 


low 


v.low 


high 


medium 


low 


high 


low 


medium medium 


low 


medium 


low 


medium 


low 



maxjmax (M. N)}= 



Now min {max (M. N)} 



high 


very high 


high 


high 


high 


very high 


high 


high 


high 


very high 


high 


high 


high 


very high 


high 


high 




v.low 


0 


0 


low 


)} = 


0 


0 


low 


0 


low 


0 


0 


low 




0 


low 


0 


0 



Wo see this four types of operations gives four types of 
4x4 fuzzy linguistic matrices. 

Now we can also find linguistic row matrix with a square or 
a rectangular fuzzy linguistic matrix which is compatible. 

Consider X = (fast, slow, very slow, just fast, very fast) to 
be the fuzzy linguistic row matrix. 
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Take 



slow 


medium 


fast 


very slow 


slow 


medium 


fast 


slow 


fast 


fast 


medium 


slow 


just fast 


slow 


very slow 



be a 5 x 3 fuzzy linguistic matrix. Consider min {min {X, M} } 
= (very slow, very slow, very slow). 

Consider max (min {X, M}) = (just fast, medium, fast). 

Consider min (max {X, M}) = (slow, slow, medium). 

To find max (max {X, M}) = (very fast, very fast, very fast). 



Here we describe the notion of Fuzzy Linguistic Cognitive 
Models (FLCM). 

For Fuzzy Linguistic Cognitive maps please refer [53]. 

Example 3.1: Suppose we are interested in studying the child 
labour problem. Let (Ci, C 2 , ..., C f ,j be six attributes / concepts 
associated with it. 

Ci - Child Labour 

C 2 - Good Teacher 

C 3 - School Drop out 

C4 - Poverty 

C5 - Public encouraging child labour 

Let L = { 0, often, + often, -often, very much, much, not that 
much, little, very little, more etc. }. 
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So the state vectors as well as the related fuzzy linguistic 
matrices take their values from the set L. Also the vertices of 
the fuzzy linguistic graph take their values from L. 

Now using the experts opinion we have the following fuzzy 
linguistic graph. 




Now using this linguistic graph we have the following fuzzy 
linguistic matrix M. 

C, c, c 3 c 4 c 5 

C, 0 0 0 + often + very much 

C 9 0 0 - often 0 0 

M = 2 

C 3 + often - much 0 0 0 

C 4 0 0 + often 0 0 

C 5 + very much 0 0 0 0 

Now our state vector takes values from L. The min and 
max operations on L are defined as follows: 
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Let min {0, a ; } = 0 and max { 0. a, } = a, for all a, e L. 

a; g L min {a ; , a ; } = a ; and max {a;, a;} = a ; . 
min {aj - a ; } = - ai and max { ai, -a, } = a ; for all a, e L. 

For instance min {often, very often} = often 
max {often, very often} = very often 
min {much, often} = often, 
max {much, often} = much 
min {-much, often} = - much 
and max {-much, often} = often. 

Like this operations on L are performed. 

Now we find XjM using as before min {min (ai, my)} where 
x, = (ai, a 2 , ..., a 6 ) and M = (my); my, a; e L, 1 < i, j < 6. 

Thus XjM = (+ often, 0, + often, + often, 0) x 



0 0 0 + often + very much 

0 0 - often 0 0 

+ often - much 0 0 0 

0 0 + often 0 0 

+ very much 0 0 0 0 



= (+ often, 0, + often, + often, + often) 
leading to a fixed point. 

Example 3 . 2 : Consider the problem of finding the seven transit 
system which includes the level of service and the convenience 
factors. We have the following eight attributes. 

Ci - Frequency of the service along a route 
C 2 - In-vehicle travel time along the route 
C3 - Travel fare along the route 
C 4 - Speed of the vehicles along the route 
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C5 - Number of intermediate points in the route. 

C 6 - Waiting time 

C 7 - Number of transfers in the route 

The fuzzy linguistic terms associated with Q, C2, C 7 as 
well as the problem are L = { 0 , often, always, a little, much, 
very much, usually, some times } . 

We now give the fuzzy linguistic graph whose vertices are 
Ci, C2, . . ., C 7 and the edges take values from L are as follows. 






C, 


c 2 


C 3 


c 4 


c 5 


c 6 


c 7 


C, 


0 often 
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0 


0 


often 
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c 2 


0 
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very much often 


0 


0 


u 

11 
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0 


0 


very much much 


0 


some 


C 4 
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much 
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C 5 
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often 
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0 


0 


much 


C 6 
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0 


0 


0 
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c 7 
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0 


0 


0 


much 


0 


0 
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Let X = (often, 0, some, 0, often, some, 0) is the fuzzy 
linguistic state vector. 

To study the effect of X on M using max {max {X, M}} 
after updating max {max {X, M}} is (often, often, some, very 
much, often, some, much) = X t . 

Now max {max {X], M} } after updating is say X 2 ; 

X 2 = (often, very much, some, very much, often, some, very 
much). 

Now we find max { max (X 2 , M) } after updating we get say 
X 3 ; 



X 3 = (often, very much, some, very much, often, some, very 
much). 

Now suppose another expert wants to use for the same 
fuzzy linguistic state vector X and the same dynamical system, 
to study the effect using the ‘max min’ operation. 

Now max {min (X, M)} after updating we get say X h 

Xi =(often, often, some, 0, often, some, often). 

We find max {min (X, M)} after updating we get say X 2 . 
This is the way the operations are performed on the fuzzy 
linguistic matrix. 

X 2 = (often, often, some, 0, often, some often). 

For about Fuzzy linguistic relational map-model refer [53]. 

We can use any of the four types of operations {max, max] 
or {min, max] or {max, min] or {min, min] in these fuzzy 
linguistic matrices [53]. 

We will illustrate this situation by an example. 
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Example 3.3: Let us study the employee - employer fuzzy 
linguistic relational model. Suppose we have the following 
fuzzy linguistic concepts / attributes associated with the 
employee taken as the domain space 

D| - Pay with allowances and bonus 

D 2 - Only pay to employee 

D 3 - Pay with allowance to employee 

D 4 - Best performance by the employee 

D 5 - Average performance by the employee 

D 6 - Employee works for more number for hours. 

Suppose the following nodes / concepts are taken as the 
range space of the employer. 



R, - maximum profit to the employer 

R 2 - Only profit to the employer 

R 3 - Neither profit nor less to the employer. 



The fuzzy linguistic terms associated with the fuzzy 
linguistic domain and range spaces be taken as L. 

L = {0, gain, loss, no loss no gain, just gain, just loss, gain, 
heavy loss, good gain}. 
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We give the associated fuzzy linguistic matrix N of the 
fuzzy linguistic graph. 





R 1 


r 2 


D, 


good gain 


0 


d 2 


gain 


0 


N = D 3 


no loss no gain 


0 


d 4 


gain 


0 


d 5 


0 


no loss no gain 


d 6 


just loss 


0 



R 3 

just loss 
0 
0 
0 
0 
0 
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Now we find the resultant of any fuzzy linguistic vector on 
N, the dynamical system associated with the problem. 

Let X = (gain, 0, loss, 0, gain, loss) be the given fuzzy 
linguistic state vector. The effect of X on the fuzzy linguistic 
dynamical system N is as follows: 

Max {max (X, N)} = ((good gain, gain, gain)} = Y 

We find max (max {Y, N T }) = (gain, gain, loss, gain, gain, 
loss) = X! (say) after updating. 

We find max (max {Xi, N}) 

= (good gain, gain, gain) = Y, = Y. 

This is the way the operations using linguistic matrices are 
performed on these models. 

Finally for the first time we introduce the notion of fuzzy 
linguistic relation equations and describe some of the properties 
related with them. 

Let Li and L 2 be any two fuzzy linguistic sets that is both Li 
and L 2 contain fuzzy linguistic terms or L, = L 2 otherwise. Let 
R be a fuzzy linguistic relation that is to each fuzzy linguistic 
term of Li two or more fuzzy linguistic terms in L 2 are 
assigned, that is to each Xi e Li to the domain space (Li is the 
domain linguistic space) we associate a y 2 e L 2 and here the 
degree of membership is not a value between [0, 1] but a fuzzy 
linguistic value from L. 

This we will first illustrate by some example. 

Let L] = {xi, x 2 , x 3 , x 4 , x 5 , x 6 , x 7 } be seven fuzzy linguistic 
terms and L 2 = {y i, y 2 , y 3 , y 4 , y 5 } be some five fuzzy linguistic 
terms with the following relation. 
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We have fuzzy linguistic membership matrix R associated 
with the above map. 





yi 


y 2 


y 3 


y 4 


y 5 


X 1 


good 


fair 
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0 
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X 2 
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worst 
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0 


X 3 


0 


fair 


0 


0 


0 


X 4 


0 


0 


0 


bad 


good 


X 5 


best 


good 


0 


0 


0 


X 6 


0 


0 


v.bad 


0 


best 


X 7 


0 


0 


0 


fair 


0 



Thus for fuzzy linguistic sets Li and L 2 we can have a fuzzy 
linguistic membership function R. The representation by the 
diagram will be known as the fuzzy linguistic sagittal diagram 
and the fuzzy linguistic matrix will be known as the fuzzy 
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linguistic membership matrix where the membership values are 
also fuzzy linguistic terms. 

We can compose two fuzzy linguistic binary relation 
provided both of them take their fuzzy linguistic membership 
values from the same fuzzy linguistic set L. 

We will just illustrate how such compositions are made. 

Suppose 

good bad fair best 

bad fair good good 

0 good fair good 

good bad good fair 

and 

good bad good fair 0 best 

bad good best good bad 0 

best fair best 0 good bad 

0 fair good bad good good 

be two fuzzy linguistic membership matrices of a fuzzy 
linguistic binary relation with the values taken from the same 
fuzzy linguistic space L. 

Thus we have used the innovative methods of applying 
matrices by constructing in the first place fuzzy linguistic 
matrices and secondly using them in fuzzy models like Fuzzy 
Cognitive Maps, Fuzzy Relational Maps and Fuzzy Relational 
Equations. 
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Now 
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bad 


fair 


best 




oQ = 
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good 


good 
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0 


good 


fail- 


good 
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bad 


good 


fair 





good 


bad 


good 


fair 
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good 
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good 


bad 
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fair 


good 


bad 
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good 



where R is found in the following way. 





r i 


U 


r 3 


U 


r 5 


f 6 


If R = 
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f 8 


f 9 


r io 


hi 


r i2 




hi 


f 14 


r i5 


b6 


b7 


f 18 




_ r i9 


Uo 


r 21 


r 22 


r 23 


r 24 



ly = max [min (good, good), min (bad. bad), 
min (fair, best), min (best, 0)] 

= max [good, bad, fair, 0] = good 

r 2 = max [min (good, bad), min (bad, good), 
min (fair, fair), min (best, fair)} 

= max [bad, bad, fair, fair} = fair 

r 3 = max [min {good, good}, min [bad, best}, 
min {fair, best}, min {best, good}} 

= max {good, bad, fair, good} = good 
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r 4 = max {min {good, fair}, min {bad, good}, 
min {fair, 0}, min {best, bad}} 

= max { fair, bad, 0, bad } = fair. 

and so on. 



We have the following fuzzy linguistic binary relation 
which is described by the fuzzy linguistic sagittal diagram and 
the related fuzzy linguistic membership matrix. 
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The fuzzy linguistic membership matrix associated with X 





X 1 
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X 3 
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X 3 
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slow 
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fast 
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slow 
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fast 



Operations can be performed on these matrices. 
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The concept of super matrices was first introduced by [17]. 
Later the authors developed a new algebra using these matrices 
called super linear algebras [45]. Thus for the first time we use 
only super matrices to build super linear algebras. This 
application of super matrices is not only innovative and has lot 
of applications. Further one can find bulk eigen values using 
super square diagonal matrices. This method helps in easy 
comparison of eigen values and also bulk performances takes 
lesser time [45]. 

Apart from this super fuzzy matrices are used in the 
construction of super fuzzy models like New Super Fuzzy 
Relational model, New Super Fuzzy Cognitive Maps model, 
New Super Fuzzy Associative Maps model, we make use of 
them in the construction of super codes [50]. 

In this section for the first time we introduce the new notion 
of Super Fuzzy Relational Maps (SFRMs) models and they are 
applied to real world problems, which is suited for multi expert 
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problems. When we in the place of fuzzy relational matrix of 
the FRM (using single expert) use fuzzy supermatrix with multi 
experts we call the model as Super Fuzzy Relational Maps 
(SFRMs) model. 

We just recall the definition of Domain Super Fuzzy 
Relational Maps (DSFRMs) model. 

DEFINITION 4.1: Suppose we have some n experts working on a 
real world model and give their opinion. They all agree upon to 
work with the same domain space elements / attributes / 
concepts; using FRM model but do not concur on the attributes 
from the range space then we can use the special super fuzzy 
row vector to model the problem using Domain Super Fuzzy 
Relational Maps (DSFRMs) Model. 

The DSFRM matrix associated with this model will be given 
by S M 




= |>ir I s 2 m |... \s” m ] 

where each S' M is a m x t' r matrix associated with a FRM 

given by the i th expert having D /, ..., D m to be the domain 
attributes and (t\ t\ ... t‘ r ) to be the range attributes of the i th 

expert, i = 1, 2, n and S M the DSFRM matrix will be a 
special super row vector / matrix (1 <i <n). 
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However if n is even a very large value using the mode of 
programming one can easily obtain the resultant vector or the 
super hidden pattern for any input supervector which is under 
investigation. 

These DSFRMs will be known as Domain constant 
DSFRMs for all the experts choose to work with the same 
domain space attributes only the range space attributes are 
varying denoted by DSFRM models. 

We see these super row vectors (super row matrices) find 
their applications in the DSFRMs model. 

Next we proceed on to recall the notion of super FRMs with 
constant range space attributes and varying domain space 
attributes. 



DEFINITION 4.2: Let some m experts give opinion on a real 
world problem who agree upon to make use of the same space 
of attributes / concepts from the range space using FRMs but 
want to use different concepts for the domain space then we 
make use of the newly constructed special super column vector 
as the matrix to construct this new model. 

The column attributes i.e. the range space remain constant 
as Ri, ..., R s for all m experts only the row attributes for any i‘ h 
expert is If , Df ..., If ; i = 1, 2, ..., m; vary from expert to 

expert. This system will be known as the Range constant fuzzy 
super FRM or shortly denoted as RSFRM model. 

Thus the associated special super column matrix S M is 
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These super row vectors (super row matrices) find 
applications in the RSFRM models. 

DEFINITION 4.3: Suppose we have m experts who wish to work 
with different sets of both row and column attributes i.e. domain 
and range space using FRMs, then to accommodate or form a 
integrated matrix model to cater to this need. We make use of 
the super diagonal fuzzy matrix, to model such a problem. 
Suppose the first expert works with the domain attributes 
D\, ..., D' t and range attributes R\, . R\ , The second expert 

works with domain attributes D~ , and with range 

attributes Rf , ..., Rf t and so on. Thus the m th expert works with 

D'f , . D'f domain attributes and R'f , ..., R'f range 

attributes. We have the following diagonal fuzzy supermatrix to 
model the situation. We are under the assumption that all the 
attributes both from the domain space as well as the range 
space of the m experts are different. 
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The super fuzzy matrix S associated with this new model is 
given by 



r>l pi pi p2 p2 p2 pm pm pm 

K x K 2 ... K n K x K 2 . . . K n z ... K x K 2 ... 



Mi 


(0) 




(0) 


(0) 


m 2 




(0) 


(0) 


(0) 




(0) 


(0) 


(0) 







where each M, is a t , x n, matrix associated with the FRM, we 
see except, the diagonal strip all other entries are zero. 

We call this matrix as a special diagonal super fuzzy matrix 
and this model will be known as the Special Diagonal Super 
FRM Model which will be denoted briefly as (SDSFRM). 

The super diagonal matrices are used in the SDSFRM 
models. 

DEFINITION 4.4: Suppose one is interested in finding a model 
where some mn number of experts work on the problem and 
some have both domain and range attributes to be not 
coinciding with any other expert and a set of experts have only 
the domain attributes to be in common and all the range 
attributes are different. 
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Another set of experts are such that only the range 
attributes to be in common and all the domain attributes are 
different, and all of them wish to work with the FRM model 
only; then we model this problem using a super fuzzy matrix. 
We have mn experts working with the problem. 

Let the tj expert wish to work with domain attributes P ' , 
Pf ..., P' Hh) and range attributes q\, q\, q) l(h) . 

The t 2 expert works with P' , Pf ..., P' n{h) as domain 
attributes and the range attributes q\, qf qf l{t ^ and so on. 
Thus for the ti expert works with Pf Pf ..., Pf t ) as domain 
space attributes and q[, q\ , < 7 ' as range attributes (1 <i 

<m(tj) and i <n (tf). 

So with these mn experts we have an associated super FRM 
matrix. Thus the supermatrix associated with the Super FRM 
(SFRM) model is a supermatrix of the form 



S(m) = 



A n 

^hnit^nft 1) 


A 12 

yi w(d)n(/2) 




An 


A lx 

l m(t2)n(tl ) 


A 22 

^hn(t2)n(t2) 




A 2 " 

J ^ L m(t2)n(tn) 










A’" 1 

/rL m(tm)n(tl ) 


A"' 2 

rL mftni)n{t2') 




jmn 

' rL m{tm)n{tn) 



where 
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A'J 



q[ qi qL) 



p; 

pi 



m(tj ) 



>»(»,)) 



1 < i <m and 1 <j <n. S(m) is called the super dynamical FRM 
or a super dynamical system. 

This matrix Afi t )n(l , corresponds to the FRM matrix of the 
(ij) th expert with domain space attributes Pfi Pfi ..., Pfi t) and 
range space attributes q(, qf q J „ (t j , 1 < i < m and 
1 <j < n. 

Super matrices are used as the dynamical system of SFRM 
models. 

Next we proceed onto show how super row vectors (super 
row matrices) are used in the construction of SDBAM models. 
We recall the definition of SDBAM model. 

DEFINITION 4.5: Suppose a set of n experts choose to work with 
a problem using a BAM model in which they all agree upon the 
same number of attributes from the space F x which will form the 
rows of the dynamical system fanned by this multi expert BAM. 
Now n distinct sets of attributes are given from the space F v 
which forms a super row vector and they form the columns of 
the BAM model. 

Suppose all the n experts agree to work with the same set of 
t-attributes say (xj xi ... x,) which forms the rows of the synaptic 
connection matrix M. Suppose the first expert works with the pi 
set of attributes given by (y\y\ the second expert with 
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P 2 set of attributes given by ( y\ y\ . . . y' Pi ) and so on. Let the i' h 
expert with p t set of attributes given by ( y\ y\ • ■■y' p ) for i = 1, 

2, ..., n. Thus the new BAM model will have its elements from F y 
where any element in F y will be a super row vector, T = 
(y\y\ ■■■y\ I yfy 2 2 ■■■y 2 Pl I - I y"y" ...y n Pn ). Now the synaptic 

projection matrix associated with this new BAM model is a 
special row supervector M r given by 



y! y\ 



y Pi yiy 2 



y 



Pi 



M r = 



y t y i —y 



n 

Pn 



Here the elements /attributes from F x is a simple row rector 
where as the elements from F y is a super row vector. 

We call this model to be a multi expert Special Domain 
Super-vector BAM (SDBAM) model and the associated matrix is 
a special row vector matrix denoted by M r . Let X = (x/ x 2 ... xj 
gF x . Y= [y x y 2 ...y pi \ y t y 2 ...y f2 \ ... \ y t y 2 ...yj € F y . 

IfX = (xi x 2 ... x t ) £ F x is the state vector given by the expert we 
find 



XM r = Y, e F y 

YMr = X,eF x ... 

and so on. This procedure is continued until an equilibrium is 
arrived. Similarly if the expert chooses to work with Y = 

[y\y\ ■■■y l Pl I yf y 2 ■■■y 2 Pl \ - I y" y " 2 ■■■y n P J e F y then we 

find the resultant by finding 
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YM, t 1 — > then find XM r and proceed on till the system 

arrives at an equilibrium state. This model will serwe the 
purpose when row vectors from F x are a simple row vectors and 
row vectors from F y are super row vectors. 

In the following we describe how the super column vectors 
(super column matrices) we used in the construction of synaptic 
matrix of the SRBAM models. 

DEFINITION 4.6: Suppose we have a problem in which all m 
experts want to work using a BAM model. If they agree to work 
having the simple vectors from F y i.e., for the columns of the 
synaptic connection matrix i.e. there is no perpendicular 
partition of their related models matrix. 

The rows are partitioned horizontally in this synaptic 
connection matrix i.e., the m experts have distinct sets of 
attributes taken from the space F x i.e. elements of F x are super 
row vectors. The resulting synaptic connection matrix M c is a 
special super column matrix. 

Let the 1 st expert have the set of row attributes to be 
(x\ x\ ... x l ), the 2 nd expert have the set of row attributes given 

by (xf xl ... x 2 ^ ) and so on. Let the i th expert have the related 
row attributes as ( x\ x 2 ...x' q ); i = 1, 2, ..., m. 

Let the column vector given by all them is [yi ... y,J. The 
related super synaptic connection matrix 
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M c is a special super column vector / matrix. 

Suppose an expert wishes to work with a super row vector X 
from F x then X = [x\ x\...x l qi \ xf x\ ...x 2 qi | ...| x"' x"‘ ...x”' ] 

we find X o M c { — > Y e F y , YM C T = X; ^ F x , we repeat the 
same procedure till the system attains its equilibrium i.e., a 
fixed point or a limit cycle. 

This model which performs using the dynamical system M c 
is defined as the Special Super Range BAM (SRBAM) model. 

We now proceed onto describe the SDSBAM model which 
makes use of the super diagonal matrices. We describe mainly 
these models to show how these models makes use of the super 
matrices. This SDSBAM model synaptic connection matrix is a 
diagonal super matrix. 
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DEFINITION 4.7: Suppose we have n experts to work on a 
specific problem and each expert wishes to work with a set of 
row and column attributes distinct from others using a BAM 
model. Then how to obtain a suitable integrated dynamical 
system using them. Let the first expert work with (x\ x\ ... x 1 ) 
attributes along the row of the related synaptic connection 
matrix of the related BAM and (y\ y\ ... y l pi ) the attributes 

related to the column , let the second expert give the row 
attributes of the synaptic connection matrix of the BAM to be 
( x 7 x 7 ; ... x 7 u ) and that of the column be (y 7 y\ ... y ' 2 ) and so 
on. Let the i expert give the row attributes of the synaptic 
connection matrix of the BAM to be (x\ x 2 ... x‘ n ) and that of 

the column to be (y[ y\ ...y' p ) for i = 1, 2, ..., n, the 
supermatrix described by 
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where A. is the synaptic connection matrix using the BAM 
model of the i th expert, i = 1, 2, ..., n where 



y i y 2 






A*,= 



i - X, 



(0) denotes the zero matrix. 

Thus this model has only non zero BAM synaptic connection 
matrices along the main diagonal described by M D . The rest are 
zero. 

The dynamical system associated with this matrix M D is 
defined to be the Special Diagonal Super BAM (SDSBAM) 
model. 

Y et another application (or use) of super matrices is in the 
super BAM models which is described below. 

DEFINITION 4.8: Suppose we have mn number of experts who 
are interested in working with a specific problem using a BAM 
model; a multi expert model which will work as a single 
dynamical system is given by the Super BAM (SBAM) 
model. Here a few experts have both the row and columns of 
the synaptic connection matrix of the BAM to be distinct. Some 
accept for same row attributes or vectors of the synaptic 
connection matrix but with different column attributes. Some 
accept for same column attributes of the synaptic connection 
matrix of the BAM model but with different row attributes to 
find the related supermatrix associated with the super BAM 
model. The supermatrix related with this new model will be 
denoted by M s which is described in the following: 
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where A\ is the synaptic connection matrix of an expert who 
chooses to work with (x\ x 2 ... x p ) along the row of the BAM 
model and with (y{ y J 2 ... y j q ) along the column of the BAM 
model i.e., 

y[ y{ ••• y’ q , 

i 

x x 

A'j = x 2 



1 <i < m and 1 <j <n. Thus for this model both the attributes 
from the spaces F x and F y are super row vectors given by 

X— [x x x 2 ... x pi | x l x 2 ... | ••• | 




y.m m m i 

a 2 * * * A 'p m J 
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in F x and 



\r r l l l 

Y= fyi y 2 y q 



2 2 
j't y 2 



••t; 



n n n 7 

x t 2 ■■■yj 



from the space or the neuronal field F v . 

The supermatrix M s is called the synaptic connection 
supermatrix associated with the multi expert super BAM model 
(SBAM model). Now having defined the multi expert super BAM 
model we proceed on to describe the functioning of the super 
dynamical system. 

LetX= [x\ x\ ... x l pi j x\ x\ ... x 2 p2 \ ... \ x” x” ... x p ] € F x 

be the super row vector given by the expert, its effect on the 
multi super dynamical system M s . 

X o M s 1 — > Y 

= [y\ y* 2 ■■■y' q , I y, vl -y 2 J -I y" y" - y"J 

Yo M T s ^ X, € F x . 

X l0 M s ^ Y, € F y ; 

and so on and this procedure is repeated until the system attains 
a equilibrium. 

Finally we just indicate how the super matrices are used in 
super fuzzy associative memories. For more about this refer 
[48]. 

A fuzzy set is a map p : X — » [0, 1] where X is any set 
called the domain and [0, 1] the range i.e., p is thought of as a 
membership function i.e., to every element x e X, p assigns a 
membership value in the interval [0, 1]. But very few try to 
visualize the geometry of fuzzy sets. It is not only of interest but 
is meaningful to see the geometry of fuzzy sets when we discuss 
fuzziness. Till date researchers over looked such visualization 
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[24], instead they have interpreted fuzzy sets as generalized 
indicator or membership functions; i.e., mappings p from 
domain X to range [0, 1], But functions are hard to visualize. 
Fuzzy theorist often picture membership functions as two- 
dimensional graphs with the domain X represented as a one- 
dimensional axis. 

The geometry of fuzzy sets involves both domain X = 
(xi,..., x n ) and the range [0, 1] of mappings p : X — > [0, 1]. The 
geometry of fuzzy sets aids us when we describe fuzziness, 
define fuzzy concepts and prove fuzzy theorems. Visualizing 
this geometry may by itself provide the most powerful argument 
for fuzziness. 

An odd question reveals the geometry of fuzzy sets. What 
does the fuzzy power set F(2 X ), the set of all fuzzy subsets of X, 
look like? It looks like a cube. What does a fuzzy set look like? 
A fuzzy subsets equals the unit hyper cube I n = [0, l] n . The 
fuzzy set is a point in the cube I n . Vertices of the cube F define 
a non- fuzzy set. Now with in the unit hyper cube I n = [0, l] n we 
are interested in a distance between points, which led to 
measures of size and fuzziness of a fuzzy set and more 
fundamentally to a measure. Thus within cube theory directly 
extends to the continuous case when the space X is a subset of 
R n . 



The next step is to consider mappings between fuzzy cubes. 
This level of abstraction provides a surprising and fruitful 
alternative to the prepositional and predicate calculus reasoning 
techniques used in artificial intelligence (AI) expert systems. It 
allows us to reason with sets instead of propositions. The fuzzy 
set framework is numerical and multidimensional. The AI 
framework is symbolic and is one dimensional with usually 
only bivalent expert rules or propositions allowed. Both 
frameworks can encode structured knowledge in linguistic form. 
But the fuzzy approach translates the structured knowledge into 
a flexible numerical framework and processes it in a manner 
that resembles neural network processing. The numerical 
framework also allows us to adaptively infer and modify fuzzy 
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systems perhaps with neural or statistical techniques directly 
from problem domain sample data. 

Between cube theory is fuzzy-systems theory. A fuzzy set 
defines a point in a cube. A fuzzy system defines a mapping 
between cubes. A fuzzy system S maps fuzzy sets to fuzzy sets. 
Thus a fuzzy system S is a transformation S: F — > I p . The n- 
dimensional unit hyper cube l n houses all the fuzzy subsets of 
the domain space or input universe of discourse X = {xi, ..., 
x n } . I p houses all the fuzzy subsets of the range space or output 
universe of discourse, Y = {y b ..., y p }. X and Y can also denote 
subsets of R n and R p . Then the fuzzy power sets F (2 X ) and F 
(2 y ) replace I n and I p . 

In general a fuzzy system S maps families of fuzzy sets to 
families of fuzzy sets thus S: I" 1 x...xl" r — »I Pl x...xP ! Flere 
too we can extend the definition of a fuzzy system to allow 
arbitrary products or arbitrary mathematical spaces to serve as 
the domain or range spaces of the fuzzy sets. We shall focus on 
fuzzy systems S: I n — > I P that map balls of fuzzy sets in I n to 
balls of fuzzy set in I p . These continuous fuzzy systems behave 
as associative memories. The map close inputs to close outputs. 
We shall refer to them as Fuzzy Associative Maps or FAMs. 

The simplest FAM encodes the FAM rule or association (A;, 
Bi), which associates the p-dimensional fuzzy set B! with the n- 
dimensional fuzzy set A;. These minimal FAMs essentially map 
one ball in F to one ball in I p . They are comparable to simple 
neural networks. But we need not adaptively train the minimal 
FAMs. As discussed below, we can directly encode structured 
knowledge of the form, “If traffic is heavy in this direction then 
keep the stop light green longer” is a Hebbian-style FAM 
correlation matrix. In practice we sidestep this large numerical 
matrix with a virtual representation scheme. In the place of the 
matrix the user encodes the fuzzy set association (Heavy, 
longer) as a single linguistic entry in a FAM ha nk linguistic 
matrix. In general a FAM system F: I n — > I b encodes the 
processes in parallel a FAM bank of m FAM rules (Ai, Bi), ..., 
(A m B m ). Each input A to the FAM system activates each stored 
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FAM rule to different degree. The minimal FAM that stores (A;, 
B ; ) maps input A to Bf a partly activated version of B,. The 
more A resembles A ; , the more B;’ resembles Bj. The 
corresponding output fuzzy set B combines these partially 

activated fuzzy sets Bj , B \ , . . . , B[ n . B equals a weighted 
average of the partially activated sets B = w, Bj + ... + w m B 1 m 
where w ; reflects the credibility frequency or strength of fuzzy 
association (A;, 13 ). In practice we usually defuzzify the output 
waveform B to a single numerical value yj in Y by computing 
the fuzzy centroid of B with respect to the output universe of 
discourse Y. 

More generally a FAM system encodes a ha nk of compound 
FAM rules that associate multiple output or consequent fuzzy 
sets B'i, ..., B ; s with multiple input or antecedent fuzzy sets A, 1 , 
..., A ; r . We can treat compound FAM rules as compound 
linguistic conditionals. This allows us to naturally and in many 
cases easily to obtain structural knowledge. We combine 
antecedent and consequent sets with logical conjunction, 
disjunction or negation. For instance, we could interpret the 
compound association (A 1 , A 2 , B), linguistically as the 
compound conditional “IF X 1 is A 1 AND X 2 is A 2 , THEN Y is 
B” if the comma is the fuzzy association (A 1 , A 2 , B) denotes 
conjunction instead of say disjunction. 

We specify in advance the numerical universe of discourse 
for fuzzy variables X 1 , X 2 and Y. For each universe of discourse 
or fuzzy variable X, we specify an appropriate library of fuzzy 
set values A/, ..., A k 2 Contiguous fuzzy sets in a library 
overlap. In principle a neural network can estimate these 
libraries of fuzzy sets. In practice this is usually unnecessary. 
The library sets represent a weighted though overlapping 
quantization of the input space X. They represent the fuzzy set 
values assumed by a fuzzy variable. A different library of fuzzy 
sets similarly quantizes the output space Y. Once we define the 
library of fuzzy sets we construct the FAM by choosing 
appropriate combinations of input and output fuzzy sets 
Adaptive techniques can make, assist or modify these choices. 
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An Adaptive FAM (AFAM) is a time varying FAM system. 
System parameters gradually change as the FAM system 
samples and processes data. Flere we discuss how natural 
network algorithms can adaptively infer FAM rules from 
training data. In principle, learning can modify other FAM 
system components, such as the libraries of fuzzy sets or the 
FAM-rule weights w,. 

We now define the notion of super fuzzy associative 
memories of different types. 

DEFINITION 4.9: We have a problem P on which n experts 
wishes to work using a FAM model which can work as a single 
unit multi expert system. Suppose all the n-experts agree to 
work with the same set of attributes from the domain space and 
they want to work with different and distinct sets of attributes 
from the range space. Suppose all the n experts wish to work 
with the domain attributes (xj X 2 ... x t ) from the cube 
/' = \0, 7] x ... x [ 0 , 1] . Let the first expert work with the range 

t —times 

attributes (y\ y\ ... y 1 ) and the second expert works with the 
range attributes (yf y\ . . . y 2 ) and so on. Thus the i ,h expert 
works with the range attributes (y\ y‘ 2 ... y’ p ), i = 1, 2, n. 

Thus the range attributes 

xr /II 1 I 2 2 2 1 I n n n \ 

Y = (y x y 2 ...y pi \ y x y 2 ... y^ | ... | y r y 2 ...y p J 
are taken from the cube I Pl+P2+ ''' +Pn = [0, 1] x \0, 1\ x . . . x [0, 1] . 

Pi + P2 + ■■ - + Pn ti mes 

This we see the range attributes are super row fuzzy vectors. 

Now the matrix which can serve as the dynamical systems 
for this FAM model is given by F R . 
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yly'i — y* 


2 2 2 
yiy 2 ---y P2 ••• 


y"y n 2 --y" Pn 


X, 












4 


4 




4 


x t 











Clearly F R is a special super row fuzzy vector. Thus F: V 15 

j Pl +p 2 +...+ Pn Suppose using an experts opinion we have a fit 
vector, A = (an a 2 , a,) ; a,- £{0,1}, then A o F R = max min 

(a„ fij); a, e A and fij e F R . Let A o F R = B = (bj), then F R o B = 
max min (fij, b'fi and so on, till we arrive at a fixed point or a 
limit cycle. The resultant fit vectors give the solution. This F R 
gives the dynamical system of the new model which we call as 
the Fuzzy Special Super Row vector FAM model (SRFAM 
model). 

DEFINITION 4.10: Suppose we have n experts working on a 
problem and they agree upon to work with the same range 
attributes and wish to work with distinct domain attributes 
using a FAM model. We built a new FAM model called the 
special super column fuzzy vector FRM model (SCFAM) and its 
related matrix is denoted by F c . The fit vectors of the domain 
space are simple fit vectors where as the fit vectors of the range 
space are super row fit vectors. 

Now we describe the special column fuzzy vector FAM, F c 
by the following matrix. The column attributes of the super fuzzy 
dynamical system F c are given by 

(y,y 2 ..- y s ) e V =[0, 7]x...x[0, 7], 

v v J 

s -times 



The row attributes of the first expert is given by 
(x\, x\, ...,x l p J, the row attributes of the second expert is given 

by (xf, x\, ...,x 2 p fi. Thus the row attributes of the i ,h expert is 

given by (x\ x\ ... x‘ p ), i = 1, 2, n. 
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I'Ve have 



y i y s 

A 

A 



to be a special super column fuzzy vector / matrix, where 

; ; i y 2 ••• y, 

a r 

A, = X 2 



i = 1, 2, ..., n is a fuzzy pi x s matrix. Suppose the expert wishes 
to work with a fit vector (say) [ x\ x\ ... x ^ j xf x\ ... x 2 \ ... \ 

x" x" ... x p ]. Then X o F c = B where B is a simple row vector 
we find F o B = \0, 1] x ... x [0, 1] = I +p ; we proceed on 

Pi+ p 2 +. . .+ p n times 

to work till we arrive at an equilibrium state of the system. 




F r = x 



pi 
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Next we proceed on to define FAM when n expert give 
opinion having distinct set of domain attributes and distinct set 
of range attributes. 

DEFINITION 4.11: Let n experts give opinion on a problem P 
and wish to use a FAM model, to put this data as an integrated 
multi expert system. 

Let the first expert give his/her attributes along the column 
as (y\ y\ ■■■ y' qi ) an d those attributes along the row as 

(x\ x\... x\). 

Let (yf y\ ■ ■■ y 2 qi ) and ( x 2 x; ... x 2 j ) be the column and 

row attributes respectively given by the second expert and so 
on. 



Thus any i th expert gives the row and column attributes as 
(x\ x\ ... x‘ p ) and (y[ y\ ... y‘ q ) respectively, i =1, 2, 3, ..., n. 

So for any i th expert the associated matrix of the FAM would 
be denoted by A' where 



Fi y 2 ••• y q 



A' = 



Now form the multi expert FAM model using these n FAM 
matrices A 1 , A 2 , ..., A" and get the multi expert system which is 
denoted by 
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This fuzzy > supermatrix F D will be known as the diagonal 
fuzzy supermatrix of the F AM and the multi expert system which 
makes use of this diagonal fuzzy supermatrix F D will be known 
as the Fuzzy Super Diagonal FAM (SDFAM) model. Now the 
related fit fuzzy supervectors of this model F x and F y are fuzzy 
super row vectors given by X = (x\ x\ ... jc^ | xf fi 2 ... x 1 ^ \ ... \ 

x" x" ... x p J er F x and Y= (y\ y\ ... j y\ j 2 2 ... yl_ \ ... \ 

y" y" ■■■ y'i ) ^F y . 

Now this new FAM model functions in the following way. 



Suppose the expert wishes to work with the fuzzy super state 
fit vector X = (x\ x\ ... x x \ x\ x 2 ... x 2 pi \ ... \ x" x" ... x" p ) 

then Y= (y\ y\...y\ x \ y[ y 2 2 ... y 2 b \ ... \ y" f 2 ... y"J eF r 
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Now F D o Y = X] g F x and X,o F D = Yj€ F y and so on. 

This procedure is repeated until the system equilibrium is 
reached. 

DEFINITION 4.12: Let us suppose we have a problem for which 
mn experts want to give their opinion. Here some experts give 
distinct opinion both for the row attributes and column 
attributes. Some experts concur on the row attributes but give 
different column attributes and a few others have the same set 
of row attributes but have a different set of column attributes. 
All of them concur to work using the FAM model. To find a 
multi expert FAM model which can tackle and give solution to 
the problem simultaneously. 

To this end we make use of a super fuzzy matrix F s which is 
described in the following. Let the mn experts give their domain 
and column attributes as follows. 

The first expert works with the domain attributes as 
(x\ x\ ... x l p fi and column attributes as (y\ y\ ... y l q J. The 
second expert works with the same domain attributes viz 
(x\ x\ ... x l ) and but column attributes as ( yf y\ ... yfi ). The 

i th expert, 1< i < n works with ( x\ x\ ... x' n ) as the domain 
attributes and (y[ y' 2 ... y' q ) as the column attributes. 

The (n + l) th experts works with the new set of domain 
attributes (x\ x 2 ... x p [ ) but with the same set of column 

attributes viz. (y\ y\ ... y l q J. 



Now the n + j th expert works with using ( x\ x\ ... x 2 pi ) as 
the domain attribute and (y\ y\ ... y‘ q ) as the column attribute 
1 <j <n. The (2n + l) th expert works with (x\ x 2 ... x p ) as the 
row attribute and (y\ y\ ... y' 9i ) as the column attribute. 
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Thus any (2n + k) th expert uses (x\ f ... x ^ ) to be the row 
attribute and ( y l y 2 ... y q k ) to be the column attribute 1< k < 
77. Thus any (t n + r) th expert works with (x\ x 2 ... x' ) as the row 
attribute (1 <t <m) and (y[ y[ ... y r q ) as the column attribute 

1 <r < 77 . 

Now as 1 <t <777 and 1 <r <n we get the FAM matrices of 
all the mn experts which is given by the supermatrix F s . 




where A- is a fuzzy matrix associated with if' expert 
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y( y J 2 



y„ 



x \ 

A-i = x\ 



1 < i <m and l<j <n. This model is known as the multi expert 
fuzzy Super FAM (SFAM) model. The fit vectors associated with 
them are super row vectors from F x and F y . 

The fit super row vector X from F x is 



X = [x\ x\ ... x l | x 2 x 2 ... x 2 



pi 



X, X n 



■■■<] 



and the 



Xe F x = / /, ' +ft+ '" +p " = [0, 1] x \0, 1] x ... x [0, 1 ] . 

Pi+P2+---+P m times 

The fit super row vector Y from F v is 

Y=(y l y 2 ...y qi \ y x ...y b \ ... I y x y 2 ... yj ; 

Y g F y = I qi+q2+ - +q " =[0, 7] x \Q, 7] x ... x [Q, 7] 

qj+qi+.-.+im times 

Thus if 



X = [X x X 2 ...X^ | ^ X, ...X pi | ... | X; x 2 ...x p J 

is the fit vector given by an expert; its effect on F s is given 
by X o F s = Y e F v ; now F s o Y = X 1 £ F x then find the effect of 
X 1 on F s ; X 1 o F s = Y 1 £ F y and so on. 
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We repeat this procedure until we arrive at a equilibrium 
state of the system. 

For exact use of the super matrices please refer [17]. We 
have shown how the abstract model is associated with a super 
matrix. 

We have used super matrices to construct super linear 
algebra [45]. We have developed super linear algebras using 
super matrices is an analogous way like linear algebra and 
derived several classical theorems. Flere we recall a few of 
them for more about the innovative use of super matrices refer 
[45]. 

DEFINITION 4.13: Let V = (V/\ ... I V n ) be a super vector space 
and T s = (Tj \ ... \ T„) be a linear operator in V. IfW= (W i \ ... \ 
WJ be a super subspace ofV; we say that W = (Wj \ ... \ WJ is 
super invariant under T if for each super vector a = (a/ \ ... 

| of) in W= (Wi | ... | W„) the super vector T fa) is in W= (Wj | 
... | WJ i.e. if TfW) is contained in W. When the super subspace 
W = (W, | ... | WJ is super invariant under the operator T s = 
(Ti \ ... j TJ then T s induces a linear operator (T s ) w on the super 
subspace W= (W/ | ... | WJ. 

The linear operator (TJ W is defined by (T s ) w (a) = Tfa) for 
a in W = (Wj | ... | WJ but (TJ W is a different object from T s = 
(Tj | ... | TJ since its domain is W not V. 

When V = (Vi | ... | V n ) is finite (n b . . ., n n ) dimensional, the 
invariance of W = (Wi ... W n ) under T s = (Ti | ... | T n ) has a 
simple super matrix inteipretation and perhaps we should 
mention it at this point. Suppose we choose an ordered basis B 
= (Bi | ... | Bn) =(«;...< I ...|a[ ...< n ) for V = (V, | ... 

V n ) such that B' = (aj ... a* | ... | a” ... a° ) is an ordered basis 

for W = (Wj | ... | W n ); super dim W = (r b ..., r n ). Let 
A = [TJ b so that 




Supermatrices and their Applications | 91 



n l 


n n 


YA\.a\ 1 . 

^ !lJl *1 1 


■ i 


_i { =l 


i n =1 



Since W = (Wi ... W n ) is super invariant under T s = (T t | 
... | T n ) and the vector T s aj = (^a’ | . . . | T n a" ) belongs to W = 

(Wi | ... | W n ) for j t < r t . This means that 






j, < r t ; t = 1, 2, . .., n. In other words Ah = (Ah | ... | Ah ) = (0 
... | 0) if j, < r, and i t > r t . 




where B t is an r, x r t matrix, C t is a r t x (n t - r t ) matrix and D t is 
an (n t - r t ) x (n t - r t ) matrix t = 1,2, n. 



In view of this we prove the following interesting lemma. 




92 



Innovative Uses of Matrices 



LEMMA 4.1: Let W = (Wi \ ... \ W n ) be an invariant super 
subspace for T s = (Tj | ... | T„). The characteristic super 

polynomial for the restriction operator 

(T s ) w - \ . . ,\(T n ) w ) divides the characteristic super 

polynomial for T s . The minimal super polynomial for 
(T) w = (( f) Wi | ...\(T n ) w ) divides the minimal super polynomial 

for T s . 

Proof: We have [T S ] B = A where B = {Bi ... B n } is a super basis 
for V = (Vi | ... | V n ); with B ; = |aj . . . a„. j a basis for Vi, this 

is true for each i, i = 1, 2, ..., n. A is a super diagonal square 

matrix of the form 




where each 





for i = 1,2, ..., n; i.e. 
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and B = [(T S ) W ] where B' is a basis for the super vector 

subspace W = (Wi | ... | W n ) and B is a super diagonal square 
matrix; i.e. 




Now using the block form of the super diagonal square matrix 
we have super det (xl - A) = super det (xl - B) x super 
det (xl - D) 

i.e. (det (xIi-Ai) ...| det (xI n -A n ) ) 

= (det (xl[ -Bj)det(xl" — Dj) | | 

dct(xl' -B n )dct(xT" -DJ). 

This proves the restriction operator (T s ) w super divides the 
characteristic super polynomial for T s . The minimal super 
polynomial for (T s ) w super divides the minimal super 
polynomial for T s . 

It is pertinent to observe that I[ , I ", 1 , represents 

different identities i.e. of different order. 

The K th row of A has the form 
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where Cf' is some r t x (n t - r t ) matrix; true for t = 1, 2, n. 
Thus any super polynomial which super annihilates A also 
super annihilates D. 

Thus our claim made earlier that, the minimal super 
polynomial for B super divides the minimal super polynomial 
for A is established. 



Thus we say a super subspace W = (Wi ... | W n ) of the 
super vector space V = (Vi | ... | V n ) is super invariant under 
T s = (T, | ... | T n ) if T S (W) c W i.e. each T,(W,) c W { ; for i = 1, 
2, ..., n i.e., if a = (ai | ... | a n ) e W 



then T s a = (Tioti | ... | T n a n ) where 

a! = x!a! + ... +x‘ a‘ ; 

11 r i r l 

2 2 2 1 

a 2 = Xj a, + ... + x‘a r 



and so on 



T s a = (tj xjaj + ...+t) x‘ a) | ...|t”x”a” + ...+t" x” a” ) . 



Now B described in the above theorem is a super diagonal 
matrix given by 
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Thus the characteristic super polynomial of B i.e. (T s ) w ) is 

g = (gl | ... | g n ) = 

((x- c ;r i ...(x-c^) eL 'i...i(x-c;r r ...(x-c” n ) e -) 
where e* = dim W, 1 for i = 1, 2, ..., K, and t = 1,2, ..., n. 

Now we proceed onto define T s super conductor of any a 
into W = (Wi| ... | W n ). 

DEFINITION 4.14: Let V= (Vj \ ... I V n ) be a super vector space 
over the field F. W = (W i | ... | W n ) be an invariant super 
subspace ofV for the linear operator T s = (T t \ ... \T n ) of V. Let 
a= (a] | ... | a,j be a super vector in V. The T-super conductor 
of a into W is the set S T (a;W) = (S T {a^W^) |...|*S , r (cc^Wf)) 

which consist of all super polynomials g = (gi \ ... \ g,J (over 
the scalar field F) such that g(TJa is in W, i.e. (gfiTdai \ ... \ 
gn(T,Ja,J gW= (W, | ... | WJ. i.e. gilTJad € WJor every i. Or 
we can equivalently define the T s - super conductor of a in W is 
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a Ti conductor of a, in Wjfor every > i = 1, 2, .... n. Without loss 
in meaning we can for convenience drop T, and write the super- 
conductor of a into Was S(a;W ) = (S(a l ;W l )\...\S{a„;W„)) . 

The collection of polynomials will be defined as super 
stuffer this implies that the super conductor, the simple super 
operator g(Tj = ( gi(T )) \ ... \ g n (TJ) leads the super vector a 
into W. In the special case W = (0 \ ... | 0), the super conductor 
is called the T s super annihilator of a. 

The following important and interesting theorem is proved. 

THEOREM 4.1: Let V = (Vj \ ... I V,J be a finite dimensional 
super vector space over the field F and let T s be a linear 
operator on V. Then T s is super diagonalizable if and only if the 
minimal super polynomial for T s has the form 

p = (pi\ ... | p n ) = ((x - c,' )...(x-c' K )\...\(x- C ")... (x - c n Kn )] 

where (cj ...c^ |...|c" ...c n K ) are such that each set 
c[,...,c K are distinct elements of F for t = 1, 2, ..., n. 

Proof: We have noted that if T s is super diagonalizable, its 
minimal super polynomial is a product of distinct linear factors. 
To prove the converse let W = (Wi | ... | W n ) be the super 
subspace spanned by all of the characteristic super vectors of T s 
and suppose W = (Wi | ... | W n ) ^ (V| | ... | V n ) i.e. each W; ^ 
Vi. By the earlier results proved there is a super vector a not in 
W = (Wj | ... | W n ) and a characteristic super value 
e j =(Cj ,...Cj ) of T s such that the super vector |3 = (T-Cjl)a 

i.e. (Pi |...|p n ) = ((T 1 -cJ i I 1 )a 1 |...|(T n -c”I n )a n ) lies in W = 
(Wi | . . . | W n ). Since (P: | ..i'\ p n ) is in W, 



p=(P! + ... + p^ipf + ... + p^i...ipr + ...+P" n ) 
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where P t = P[ +...+P^ t for t = 1,2, ..., n with T s Pi = CiPi; 1 < i < 
K i.e. (Tjpj, | . . . | T n p" ) = (c^ p^ | . . . | c” p” ) ; (1 < i t < K t ) and 
therefore the super vector 

h(T s )P = (hj (cj ) P{ + . . . +hj (c^ ) P^ | . . . 

h n (c”)P” +---+h n ( C K n )Pk„) 

= (h 1 (T I )P 1 |...|h n (T n )P n ) 

is in W = (Wi | ... | W n ) for every super polynomial h = (Iq | . . . 
h„). 

Now (x - Cj) q for some super polynomial q, where p = (pi 
. . . | p„) and q = (qi | . . . | q n ). 

Thus p = (x - Cj) q implies 
P = (Pi I - - - I Pn) 

= (( x-c j i )q 1 |...|(x-c” )q n ) 

i.e. (q! -q 1 (c’ i )|...|q n -q n (c” ))= ((x-cj i )h 1 |...|(x-c? )h n ) . 
We have 

q(T s )a-q(Cj)a = (q, (T, )ot, - q, (c^ )a t |... | 

q n ( T „K-q n (c" n K) 

= h(T s )(T s - C jI)a=h(T s )P 
= (h l (T 1 )(T 1 -c‘l 1 )a 1 |...|h n (T n )(T n -c;;i n )a n ) 

~(h 1 (T 1 )P 1 |...|h n (T n )P n ) . 

But h(T s )P is in W = (W 1 1 . . . | W n ) and since 

0 = P(T s )a = (p,(T l )a, |...|p n (T n )a n ) 

= (T s -c j I)q (T s )a 

- ((T, - c‘ I, ) q, (T, ) a i: ). . . | (T n - c” I n ) q n (T n )a n ) ; 
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the vector q(T s )a is in W. Therefore q(cj)a is in W. Since a is 
not in W we have q(c j ) = (q 1 (c] i )|...|q n (c? ))=(0|...|0) . Thus 
contradicts the fact that p = (pi | ... | p n ) has distinct roots. 

If T s is represented by a super diagonal square matrix A in 
some super basis and we wish to know if T s is super 
diagonalizable. We compute the characteristic super polynomial 
f = (fi | ... | f n ). If we can factor 

f = (ft I I fn) 

= ((x-c!) d| ...(x-c^) 4 ‘|...|(x-cr) dr ...(x-c^f») 

we have two different methods for determining whether or not T 
is super diagonalizable. One method is to see whether for each i 

= 1, 2, ..., n we can find independent characteristic super 

vectors associated with the characteristic super values c ; . The 
other method is to check whether or not 

(T, -c,I)...(T s — c k I) i.e. ((T, - c'f) . . .(T, -c^) 

I — I (T n -c”)I n ...(T n -c^I n )) 

is the super zero operator. 

Several other interesting results in this direction can be 
derived. Now we proceed onto define the new notion of super 
independent subsuper spaces of a super vector space V. 

DEFINITION 4.15: Let V = (Vj \ ... I V,J be a super vector space 
over F. Let W x = (W x x \ ...\W"), W 2 =(W 2 n |.. .\W 2 n )... 
W K = (W l K | ...\W”) be K super subspaces of V. We say Wj, ..., 
Wk are super independent if a! + ... + a K = 0; a , ^ W t implies 
each a, = 0. 



a i = (ai\...\a i n )eW i =(W i l \...\W i "); 
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true for i = 1, 2, ..., K. IfWi and W 2 are any two super vector 
subspaces of V = (V t \ ... \ VJ, we say W l = (W* \ ...\ W") and 
W 2 -{Wf...\W") are super independent if and only if 
W x nW 2 = (Wi nW 2 \...\W" nW”) = (0 \ 0 \ ... \ 0). IfW,, W 2 , 
..., W K are K super subspaces of V we say W 2 , W 2 , ..., W K are 
independent if W 1 n W 2 n ... n W K = (W/ nW 2 n...Wl- 1...| 
W" r\W" r\...r\W'f) = (0 | ... | 0) . The importance of super- 
independence in super subspaces is mentioned below. Let 

W’ = Wi + ... + W’ k 

= (W? + ... + Wz\...\W" +...+W£) 

= W 1 ...\w;) 

W/ is a subspace Vi and W’= W[+ ... + W^ true for i = 1, 2, ..., 
n. Each super vector a in W can be expressed as a sum 
a-(a[\...\a' n )-{{a\+... + a' K )\...\a"+... + a n K )) i.e. each 
a' - a\ + ...+ a' K ; ol e W,. If Wi, W 2 , ..., Wk are super- 
independent, then that expression for a is unique; for if 

a = (A + ... + P K ) = (A‘+...+K I ...| A" +-+PD 
e W t ; i = 1, 2, ..., K. /? = + ... + f' n then 

«-« = (0|...|0) = ((a\- p\)+...+(a x K - pf)\...\ 



hence each a‘ - ff = 0 ; 1 <i <K; t = 1, 2, ..., n. Thus W 2 , W 2 , 
..., W K are super independent so we can operate with super- 
vectors in W as K- tuples ((a, 1 , . . . , a l K ); . . . , (a” , . . . , a n K ); a\ e W t ; 1 

< i < K; t = 1, 2, ..., n. in the same way we operate with R K as 
K-tuples of real numbers. 

LEMMA 4.2: Let V = (V 1 | ... | V„) be a finite (nj, ..., n„) 
dimensional super vector space. Let W 2 , ..., Wk be super 
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subspaces of V and let W = (W l 1 +... W' K \ ...\Wj + ... + W") . The 
following are equivalent 

(a) W i, ..., W K are super independent. 

(b) For each j; 2 <j <K, we have Wj n (Wi + ... + Wj-i) = 

{( 0 \...\ 0 )} 

(c) If B; is a super basis ofWi, 1 <i <K, then the sequence 
B = (B i ... B k ) is a super basis for W. 

The proof is left as an exercise for the reader. In any or all of the 
conditions of the above stated lemma is true then the supersum 
W = Wi + ... + W K = (W/ + ... + W^|...|W 1 n +... + W”) where 
W, = ( W, 1 | ... | W" ) is super direct or that W is a super direct 
sum of Wi, ..., Wk i.e. W = Wi © ... © Wk i.e. 
(W, 1 ®...®Wk|...| W, n ©...©W”) - If each of the W ; is (1, ..., 
1) dimensional then W = Wi © ... © W n 

=(w 1 1 ©...©w,;i...|w 1 n ©...© w n n ). 

DEFINITION 4.16: Let V = (V ) 1 ... I V n ) be a super vector space 
over the field F; a super projection of V is a linear operator E s 
on V such that E 2 - E s i.e. E s = ( E ] \ ... \ E„) then 

E 2 - (Ef\...\E 2 ) = (Ei | ... | E n ) i.e. each E,is a projection on 
i = 1, 2, ..., n.. Suppose E s is a projection on V and R = (Rfi 
... | R„) is the super range ofE s and N = (Nj\ ... \ N„) the super 
null space or null super space of E s . The super vector fi = (fifi 
... fi n ) is in the super range R = (If \ ... R,j if and only if E s fi 
= fi i.e. if and only if (Eifii \ ... \ E n fi n ) = (fit \ ... \ fi„) i.e. each 
Efii = fit for i = 1, 2, ..., n.. If fi = E s a i.e. fi= (fifi ... \ fi„) = 
(E/a/ | ... | E„a„) where the super vector a = (afi ... \ a„) then 
E s fi = E 2 a = E s a - fi . Conversely if fi = (fit \ ... \ fi„) = E s fi 
= (E^! | ... | E n fi n ) then fi = (fii \ ... \ fi n ) is in the super 
range ofE s . Thus V = R®N i.e. V = (V, \ ... | V„) = (R, ®N, \ 
... | R„ (DN n ). 

Further the unique expression for = (afi ... \ a„) as a sum 
of super vectors in R and N is a = E s a + (a - E s a) i.e. a, = 
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EjOCi + (dj- EiOCi) for i = 1, 2, ..., n. From what we have stated it 
easily follows that ifR and N are super subspace ofV such that 
V = R@N i.e. V= (V, \ ... \ VJ = (R, ®N, \ ... \ R n ®NJ then 
there is one and only one super projection operator E s which 
has super range R = (Rj \ ... \ R„) and null super space N = (Ni 
| ... | N„). That operator is called the super projection on R 
along N. 

Any super projection E s is super diagonalizable. If 
{(a, 1 ., .a' | ...| a". ..a") is a super basis for R = (Ri \ ... \ R n ) 

and (al +1 ...a' h ...| a" +l ...a" ) is a super basis for N = (Nj \ ... 

| N n ) then the basis B = a) h )\...\a" ...a" ) = (Bj \ ... \ B n ) 

super diagonalizes E s . 




= ([^]|-PJ S „) 



where I, is a r, x r, identity matrix; t = 1, 2, ..., n. Thus super 
projections can be used to describe super direct sum 
decompositions of the super vector space V= (Vj\ ... | V„). 

THEOREM 4.2: Let F be the field of real numbers or the field of 
complex numbers. Let Abe a super diagonal matrix of the form 
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be a (ni xrii, ■■■, n„ xn n ) matrix over F. The super function g = 
(gi | ... | g,J defined by g(X, Y) = Y AX is a positive superform 
on the super space (F"‘ xl | ... | F" ,,xl ) if and only if there exists 
an invertible super diagonal matrix 




Each P t is a n ,• x n, matrix i = 1, 2, n with entries from F 
such that A = P P ; i.e., 
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Definition 4.17: Let 




be a superdiagonal matrix with each A l a n, x 77, matrix over the 
field F; i =1, 2, n. The principal super minor of A or super 
principal minors of A (both mean the same) are scalars 



A* (A) = (A ti (4)| ••• |A k (A„)) 

defined by 




for 1 <k t <n, and t = 1, 2, n. 
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Several other interesting properties can also be derived for 
these superdiagonal matrices. 

We give the following interesting theorem and the proof is 
left for the reader. 

THEOREM 4.3: Let f = (fi I ... I fij be a superform on a finite 
(ni, ..., n,J dimensional supervector space V= (Vj \ ... | VJ and 
let Abe a super diagonal matrix of f in an ordered superbasis B 
= (B i \ ... | B„). Then f is a positive superform if and only if A = 
A * and the principal super minor of A are all positive, 
i.e., 




Note: The principal minor of (Ai | ... f A n ) is called as the 
principal superminors of A or with default of notation the 
principal minors of {A b ..., A n } is called the principal super 
minors of A. 

T s = (Ti | ... | T n ) a linear operator on a finite (ni, ..., n n ) 

dimensional super inner product space V = (V i | | V n ) is said 

to be super non-negative if T s = T* 

i.e., 

(T, I ... I T n )= (T; I ... | T n *) 



i.e. T ; = T* for i = 1, 2, . .., n and 
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(T s a | a) = ((Ti ai | a0 | ... | (T n a n | a n )) > (0 | ... |0) 
for all a = (a i | ... | a n ) in V. 

A super positive linear operator is one such that T s = T s * and 

(T a | a) = ((Ti a, | a,) | ... | (T n a n | a n ))>(0 | ... |0) 

for all a = (a i | ... | a n ) * (0 | ... |0). 

Several properties enjoyed by positive operators and non 
negative operators will also be enjoyed by the super positive 
operators and super non negative operators on super vector 
spaces, with pertinent and appropriate modification. 
Throughout the related matrix for these super operators T s will 
always be a super diagonal matrix A of the form 




where each A! is a n ; x n ; square matrix, 1 < i < n, A = A* and 
the principal minors of each A! are positive; 1 < i < n. 

Now we just mention one more property about the super 
forms. 

THEOREM 4.4: Let f= (fi \ ... f) be a super form on a real or 
complex super vector space V = (Vi \ ... \ V„) and 

ja' . . .or* | . . . | a" ... a" j a super basis for the finite dimensional 

super subvector space W = (Wi | ... | W n ) ofV=(Vj\ ... j V„). 
Let M be the super square diagonal matrix where each M-, in M; 
is a r t x r ; super matrix with entries (1 < i < n). 
- /(a[, a)), i.e. 
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used in this construction. Hence we may refer to E s as the super 
projection of V on W that is determined by the direct sum 
decomposition. 



V = W 0 W'; 

(Vi | ... | v n ) = (w,®w;|... |w n ®w;). 

Note that E s is a super orthogonal projection if and only if 
W' = W 1 = (W, 1 1 ... W n x ) . Now we proceed onto develop the 
analogous of spectral theorem which we call as super spectral 
theorem. 

Finally these super matrices are used to develop the notion 
of Leontief economic super models. 

Matrix theory has been very successful in describing the 
interrelations between prices, outputs and demands in an 
economic model. Here we just discuss some simple models 
based on the ideals of the Nobel-laureate Massily Leontief. Two 
types of models discussed are the closed or input-output model 
and the open or production model each of which assumes some 
economic parameter which describe the inter relations between 
the industries in the economy under considerations. Using 
matrix theory we evaluate certain parameters. 

The basic equations of the input-output model are the 
following: 



" a ll 


a i2 


•' a in ^ 




|V 




|V 


a 21 


a 22 


'• a 2n 




P 2 


= 


p 2 


V a nl 


a „2 ' 


a nn J 




V.P n ) 




vPny 



each column sum of the coefficient matrix is one 

i. pi> 0, i= 1, 2, ..., n. 

ii. ay > 0, i , j = 1, 2, n. 
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iii. ay + a 2 j +...+ a nj = 1 

for j = 1, 2 , n. 




are the price vector. A = (ay) is called the input-output matrix 
Ap = p that is, (1 - A) p = 0. 

Thus A is an exchange matrix, then Ap = p always has a 
nontrivial solution p whose entries are nonnegative. Let A be an 
exchange matrix such that for some positive integer m, all of the 
entries of A m are positive. Then there is exactly only one 
linearly independent solution of (I - A) p = 0 and it may be 
chosen such that all of its entries are positive in Leontief open 
production model. 

In contrast with the closed model in which the outputs of k 
industries are distributed only among themselves, the open 
model attempts to satisfy an outside demand for the outputs. 
Portions of these outputs may still be distributed among the 
industries themselves to keep them operating, but there is to be 
some excess some net production with which to satisfy the 
outside demand. In some closed model, the outputs of the 
industries were fixed and our objective was to determine the 
prices for these outputs so that the equilibrium condition that 
expenditures equal incomes was satisfied. 

X; = monetary value of the total output of the i th industry. 

dj = monetary value of the output of the i th industry needed to 
satisfy the outside demand. 
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a,, = monetary value of the output of the i th industry needed by 
the j th industry to produce one unit of monetary value of its own 
output. 

With these qualities we define the production vector. 




the demand vector 




and the consumption matrix, 




By their nature we have 

x > 0, d > 0 and C > 0. 



From the definition of Gy and Xj it can be seen that the quantity 
G;i Xi + G l 2 X 2 +...+ G ik X k 

is the value of the output of the i th industry needed by all k 
industries to produce a total output specified by the production 
vector x. 
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Since this quantity is simply the i th entry of the column 
vector Cx, we can further say that the i th entry of the column 
vector x - Cx is the value of the excess output of the i th industry 
available to satisfy the outside demand. The value of the outside 
demand for the output of the i th industry is the i th entry of the 
demand vector d; consequently; we are led to the following 
equation: 



x - Cx = d or 
(I - C) x = d 

for the demand to be exactly met without any surpluses or 
shortages. Thus, given C and d, our objective is to find a 
production vector x > 0 which satisfies the equation (I - C)x = 
d. 

A consumption matrix C is said to be productive if (1 - C) _1 
exists and (1 - C) _1 > 0. 

A consumption matrix C is productive if and only if there is 
some production vector x > 0 such that x > Cx. 

A consumption matrix is productive if each of its row sums 
is less than one. A consumption matrix is productive if each of 
its column sums is less than one. 

Now we will formulate the Smarandache analogue for this, 
at the outset we will justify why we need an analogue for those 
two models. 

Clearly, in the Leontief closed Input - Output model, 

Pi = price charged by the i th industry for its total output in reality 
need not be always a positive quantity for due to competition to 
capture the market the price may be fixed at a loss or the 
demand for that product might have fallen down so badly so 
that the industry may try to charge very less than its real value 
just to market it. 
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Similarly ay > 0 may not be always be true. Thus in the 
Smarandache Leontief closed (Input-Output) model (S-Leontief 
closed (Input-Output) model) we do not demand p; > 0, pi can 
be negative; also in the matrix A = (ay), 



aij + a 2 j +. . .+a k j * 1 

so that we permit ay's to be both positive and negative, the only 
adjustment will be we may not have (I - A) p = 0, to have only 
one linearly independent solution, we may have more than one 
and we will have to choose only the best solution. 

As in this complicated real world problems we may not 
have in practicality such nice situation. So we work only for the 
best solution. 

Here we introduce a input-output model which has some p 
number of input-output matrix each of same order say n x n 
functioning simultaneously. 

We shall call such models as input - output super row 
matrix models and describe how it functions. Suppose we have 
p number of n x n input output matrix given by the super row 
matrix A = [A! | ... | A n ] where each A! is a 

n x n input output matrix which are distinct. 



A = [Ai | ... |AJ 
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where a) j + a!,j + ...+ af=l;t = 1, 2, .... p and j = 1, 2, ..., n. 
Suppose 
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P = 



f 1 


2 




P\ 


Pi 


Pi 




Pi 
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2 






Ip„ 


Pn 




Pj 



[Pi 



be the super column price vector then 
A * P = P, the (product) * is defined as A * 




P = P that is 



[A 1 P 1 | • • • | A p P p ] = [Pj | ... | P p ] 
A*P = P 



that is 

(I - A) P = (0 | ... | 0) 

i e., ((I - A t ) P i | | (I — A p ) P p ) = (0 | ... | 0). 



Thus A is an super-row square exchange matrix, then AP = 
P always has a row column vector solution P whose entries are 
non negative. 

Let A = [Ai | ... | A n ] be an exchange super row square 
matrix such that for some positive integer m all the entries of 
A m i.e. entries of each A™ are positive for m; m = 1, 2, ..., p. 
Then there is exactly only one linearly independent solution of 

(I - A) P = (0 | ... | 0) 

i.e., ((I - A0 Pi | ... | (I - A p ) P p ) = (0 | ... | 0) 



and it may be choosen such that all of its entries are positive in 
Leontief open production super model. 

Note this super model yields easy comparison as well as this 
super model can with different set of price super column vectors 
and exchange super row matrix find the best solution from the p 
solutions got from the relation 



(I - A) P = (0 | ... | 0) 

i e., ((I - A0 Pi | ... | (I - A p ) P p ) = (0 | ... | 0) . 
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Thus this is also an added advantage of the model. It can 
study simultaneously p different exchange matrix with p set of 
price vectors for different industries to study the super 
interrelations between prices, outputs and demands 
simultaneously. 

Suppose one wants to make a study of interrelation between 
prices, outputs and demands in an industry with different types 
of products with different exchange matrix and hence different 
set of price vectors or of many different industries with same 
type of products its interrelation between prices, outputs and 
demands in different locations of the country were the economic 
status and the education status vary in different locations, how 
to make a single model to study the situation. In both the cases 
one can make use of the super input-output model the relation 
matrix which is a input-output super diagonal mixed square 
matrix, which will be described presently. 

The exchange matrix with p distinct economic models is 
used to describe the interrelations between prices, outputs and 
demands. Then the related matrix A will be a super diagonal 
mixed square matrix 




Ai, ..., A p are the exchange matrices describing the p-economic 
models. Now A acts as integrated models in which all the p 
entities function simultaneously. Now any price vector P will be 
a super mixed column matrix 
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where each 




fort = 1, 2, p. 

Here each A t is a n t x n t exchange matrix; t = 1,2, . . p. AP 
= P is given by 




AP = 




i.e. A t P t = P t for every t = 1, 2, p. i.e. 
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0 


0 




0 


0~ 


0 




o 










o 


0 




o 



Thus AP = P has a nontrivial solution 



P = 



V 



v p py 



whose entries in each P t are non negative; 1 < t < p. 

Let A be the super exchange diagonal mixed square matrix 
such that for some p-tuple of positive integers m = (mi, . . m p ), 
A”' is positive; 1 < t < p. Then there is exactly only one 
linearly independent solution; 



(I-A)P 



0 


0 




o 


o’ 


0 




o 










0~ 


0 




o 



and it may be choosen such that all of its entries are positive in 
Leontief open production super model. 

Next we proceed on the describe super closed row model 
(or row closed super model) as the super closed model (or 
closed super model). 

Here we have p sets of K industries which are distributed 
among themselves i.e. the first set with K industries distributed 
among themselves, the second set with some K industries 
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distributed among themselves and so on and the p set with some 
K industries distributed among themselves. It may be that some 
industries are found in more than one set and some industries in 
one and only one set and some industries in all the p sets. This 
open super row model which we choose to call as, when p sets 
of K industries get distributed among themselves attempts to 
satisfy an outside demand for outputs. 

Portions of these outputs may still be distributed among the 
industries themselves to keep them operating, but there is to be 
some excess some net production with which they satisfy the 
outside demand. In some super closed row models the outputs 
of the industries in those sets which they belong to were fixed 
and our objective was to determine sets of prices for these 
outputs so that the equilibrium condition that expenditure equal 
income was satisfied for each of the p sets individually. 

Thus we will have 

x* = monetary value of the total output of the i th 
industry in the t th set 1 < i < K and 1 < t < p. 

dj = monetary value of the output of the i th 

industry of the t th set needed to satisfy the 
outside demand, 1 < t < p, I = 1, 2, ..., K. 

ajj = monetary value of the output of the i th 

industry needed by the j th industry of the t th 
set to produce one unit of monetary value of 
its own output, 1 < i < K; 1 < t < p. 

With these qualities we define the production super column 
vector 
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The demand column super vector 




and the consumption super row matrix C = (Cj ... C ) 




X> 



and 0(0 | ... | 0). 
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For the t th set from the definition of al and x* it can be seen 
that the quantity 

a n Xj +a i2 x 2 +... + a iK x K 

is the value of the i th industry needed by all the K industries (of 
the set t) to produce a total output specified by the production 
vector X t . Since this quantity is simply the i th entry of the 
column vector C t X t we can further say that the i th entry of the 
column vector X t - X t C t is the value of the excess output of the 
i th industry (from the t th set) available to satisfy the outside 
demand. 

The value of the outside demand for the output of the i th 
industry (from the t th set) is the i th entry of the demand vector d t ; 
consequently we are lead to the following equation for the t th set 
X t - C t X, = d, or (I - C t )X t = d t for the demand to be exactly met 
without any surpluses or shortages. Thus given C t and d, our 
objective is to find a production vector X t > 0 which satisfies the 
equation 

(I - C t )X t = d t , 

so for the all p sets we have the integrated equation to be 

(I - C)X = d 

i.e., [(I - COXt | ... |(I-C P )X P ] 

= (di | ... | d p ) . 

The consumption super row matrix C = (Q | . . . | C p ) is said to 
be super productive if 

(i-c)- i =[(i-c 1 r i |...|(i-c p r i ] 

exists and 

(I-C) -1 = [(I-Cj) -1 | ... | (I-C p )~'] > [0| ... |0] . 
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A consumption super row matrix is super productive if and 
only if for some production super vector 





f x ri 




"o' 


X = 




> 










V 5 , 



such that X > CX i.e. [Xi | | X p ] > [C1X1 | . . . | CpX p ]. 

A consumption super row matrix is productive if each of its 
row sums is less than one. A consumption super row matrix is 
productive if each of its column super sums is less than one. 
The main advantage of this super model is that one can work 
with p sets of industries simultaneously provided all the p sets 
have same number of industries (here K). This super row model 
will help one to monitor and study the performance of an 
industry which is present in more than one set and see its 
functioning in each of the sets. Such a thing may not be possible 
simultaneously in any other model. 

Suppose we have p sets of industries and each set has 
different number of industries say in the first set output of K4 
industries are distributed among themselves. In the second set 
output of K 2 industries are distributed among themselves so on 
in the p th set output of K p -industries are distributed among 
themselves the super open model is constructed to satisfy an 
outside demand for the outputs. Here one industry may find its 
place in one and only one set or group. Some industries may 
find its place in several groups. Some industries may find its 
place in every group. To construct a closed super model to 
analyze the situation. 

Portions of these outputs may still be distributed among the 
industries themselves to keep them operating, but there is to be 
some excess some net production with which to satisfy the 
outside demand. 
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monetary value of the total output of the i th 
industry in the t th set (or group), 
monetary value of the output of the i th industry 

of the group t needed to satisfy the outside 
demand. 

monetary value of the output of the i th industry 

needed by the j th industry to produce one unit 
monetary value of its own output in the t th set or 
group, 1 < t < p. 

With these qualities we define the production super mixed 
column vector 




and the demand super mixed column vector 



d = 




x; 

d‘ = 
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vector X t (1< t < p). Since this quantity is simply the i th entry of 
the super column vector in 




= [CiXj I ... I CpXp]* 

we can further say that the i th entry of the super column vector 
X t - CX t in 

' x t - C,X p 
X-CX = : : 

v X r- C r X r 

is the value of the excess output of the i th industry available to 
satisfy the output demand. 

The value of the outside demand for the output of the i th 
industry (in t th set / group) is the i th entry of the demand vector 
d t ; consequently we are led to the following equation 



X t - C t X t = d t or (l t - C t ) X t = d t , (1 < t < p), 

for the demand to be exactly met without any surpluses or 
shortages. 

Thus given C t and d t our objective is to find a production 
vector X t > 0 which satisfy the equation (I t - C t )X t = d. 

The integrated super model for all the p-sets (or groups) is 
given by X - CX = d i.e., 
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or 



( Xi -CjX, ~j 




ro 


X 2 -C,X 2 




dj 








X -cx 




d 


VP p v J 




^ pJ 



i.e., 



(Ii-C.) 


0 




0 


0 


I -C 
±2 ^ 2 














0 


0 




Ip-Cp. 



'x,' 



vX 



r A \ 

d, 



v d ry 



(I,-C 1 )X 1 



( I r-C p )X 



f , A 



? J 



d , 



V d py 



where 1 is a K t x K t square identity matrix t = 1, 2, . . p. 



Thus given C and d our objective is to find a production 
super column mixed vector 





M 




(T 


X = 




> 










0 

v / 



which satisfies equation (I — C) X = d 



I— - 

1 

_o 

* 




V 




= 




k-Cp)x P J 




A, 



A consumption super diagonal matrix C is productive if 
(1 - C) 1 exists and i.e., 
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a.-c.r 1 


0 




0 


0 


(i 2 -c 2 r‘ 




0 










0 






dr-Cp)" 



exists and 



(ii-Qr 1 


0 




0 


0 


(i 2 -c 2 r‘ 




0 










0 






dp-V 



0 


0 




0 


1 o 


0 




o 










o 






o 



A consumption super diagonal matrix C is super productive if 
and only if there is some production super vector 





[X] 




V 


X = 




> 






w 




0 

v y 



such that 





[A] 




f 


X > CX i.e. 




> 






w 




v CpXp J 



A consumption super diagonal mixed square matrix is 
productive if each row sum in each of the component matrices is 
less than one. A consumption super diagonal mixed square 




Supermatrices and their Applications | 125 



matrix is productive if each of its component matrices column 
sums is less than one. 

The main advantage of this system is this model can study 
different sets of industries with varying strength simultaneously. 
Further the performance of any industry which is present in one 
or more group can be studied and also analysed. Such 
comprehensive and comparative study can be made using these 
super models. 

Finally we develop the notion of super special codes. 

Flere we define two new classes of super special row codes 
using super row matrix and super mixed row matrix. 



DEFINITION 4.18: Suppose we have to transform some n set of 
k lt ..., k n message symbols a I I a I 2 ...a I ki , a 2 a 2 2 ...a k ,..., a"a"...a k , 

a\ ef f ; 1 <t < n and 1 <i < k, (q a power of a prime) as a set 

of code words simultaneously into n-code words such that each 
code word is of length n„ i = 1, 2, n and nj-ki = n 2 ~ k 2 = 
... = n„ - k„ = m say i.e., the number check symbols of ever\> 
code word is the same i.e., the number of message symbols and 
the length of the code word may not be the same. That is the 
code word consisting of n code words can be represented as a 
super row vector; 



1 1 


1 


2 2 


2 


n n 


n 1 


x 1 x 2 . 


. X n 

n I 


Xj x 2 


..x n 

n 2 


. . X]X 2 . 


J 



h, > k i} 1 <i <n. In this super row vector x' = ah , 1 <j < k i = 
1, 2, ..., n and the remaining n ,• - k t elements x' k +I x k+2 ... x' n 
are check symbols or control symbols; i = 1, 2, ..., n. 

These n code words denoted collectively by x s will be 
known as the super special row code word. 
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As in case of usual code, the check symbols can be obtained 
in such a way that the super special code words x s satisfy a 
super system of linear equations; I I s x J = (0) where H s is a 

super mixed row matrix given by H s = [Hi | H 2 ... | H n ] where 
each Hj is a m x rtj matrix with elements from F q , i = 1,2, . . ., n, 
i.e., 

H s xJ = [Hj | H 2 | ... |H n ] [x[ x^ ... <] T 
- Hi(x:) T | h 2 ( x :) t |... |H 2 (x s n ) T 
= [|(0) I (0) I ... I (0)] 

i.e., each Hj is the partity check matrix of the code words x‘ ; i 

= 1, 2, ..., n. H s = [Hi | H 2 ... | H s ] will be known as the super 
special parity check super special matrix of the super special 
row code C s . C s will also be known as the linear [(nj n 2 . . . n n ), 
(k i k 2 ... k n )] or [(n b ki), (n 2 , k 2 ), ..., (n n , k n )] super special row 
code. 

If each of the parity check matrix Hj is of the form 
( A 1 > I n,-k i );i = h2, ..., n. 



H s = [Hi|H 2 | 


- |H„] 






_( A l’ I n,-k 1 ) ( A 2>!n 2 -k 2 ) ••• 


( A nJ„ n -k n )’ 


.... 


(I) 



C s is then also called a systematic linear ((n, n 2 . . . n n ), (ki k 2 . . . 
k n )) super special code. 

If q = 2 then C s is a super special binary row code; (ki + . . . 
+ k n ) | (nj + n 2 + . . . + n n ) is called the super transmission(or 
super information) rate. 

It is important and interesting to note the set C s of solutions 
x s of H s xJ = (0) i.e., known as the super solution space of the 
super system of equations. Clearly this will form the super 
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special subspace of the super special vector space over F q of 
super special dimension (k! k 2 . . . k n ). 

C s being a super special subspace can be realized to be a 
group under addition known as the super special group code, 
where H s is represented in the form given in equation 1 will be 
known as the standard form. 



Now we will illustrate this super special row codes by some 
examples. 



Example 4.1: Suppose we have a super special binary row code 
given by the super special parity check matrix H s = [Hi | H 2 1 II,] 
where 



H,= 



1110 0 
0 0 0 1 0 , 
1 0 0 0 1 



and 



H, 



0 0 0 1 1 0 0 
0 1 0 0 0 1 0 
1 0 0 0 0 0 1 



h 3 = 



1 0 0 0 0 1 0 0 
0 0 0 1 0 0 1 0 
0 0 1 0 1 0 0 1 



i.e., the super row matrix associated with the super special code 
is given by 



"0 


1 


1 


1 


0 


0 


0 


0 


0 


1 


1 


0 


0 


1 


0 


0 


0 


1 


0 


0 


1 


0 


0 


0 


1 


0 


0 


1 


0 


0 


0 


1 


1 


0 


0 


0 


0 


0 


1 



1 0 0 0 0 1 0 0 
0 0 0 1 0 0 1 0 
0 0 1 0 1 0 0 1 
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= [(A b I 3 ) I (A 2 , 1 3 ) I (A 3 , 13)] ; 

i.e., we have given the super special code which is a binary row 
code. The super special code words are given by 



x s =[aj a l 2 a* x^ xl; 

3 3 3 3 3 3 

a i a 2 a 3 a 4 a 5 X 6 



I a? a; a" a 2 4 x 2 s x 2 x 2 

X s] = [ X s I X s I X s]' 



H s x^ = (0) gives 3 sets of super linear equations i.e., H s x^ = 
(0) is read as 



[H, | H 2 1 H 3 ] [xl|xf Ix^^r^fx:) 1 | H 2 (x s 2 ) T |h 3 (xJ) 



[(0) | (0) | (0)]; 



i.e., 

^( x ;) 1 -( o ) 



results in linear equations solving which we can get the code 
words. 

Thus the number of super special row code words in this 
example of the super special code C s is 8 x 16 x 32. 

The super transmission rate is 12/21. Thus this code has 
several advantages which will be enumerated in the last chapter 
of this book. We give yet another example of super special code 
in which every super special code word is a super row vector 
and not a super mixed row vector. 

Now we proceed on to define the notion of super special 
row repetition code. 

Definition 4.19: Let C s = [c' C; | ... C"] be a super 
special row code in which each of the C‘ s is a repetition code, i 
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= 1, 2, ..., n, then we define C s to be a super special repetition 
row code. Here if H s = [H/\H 2 \ ... \H,J is the super special parity 
check matrix of C s , then each Hi is a t - 1 x t matrix that is we 
have 

~1 1 0 ••• 0~ 

1 0 1 ••• 0 

H 1 = H 2 =...=H n =. . . . 

1 0 0 ... 1 , 

L Jt-i x t 

is the parity check matrix. The super special code words 
associated with C s are just super row vectors only and not super 
mixed row vectors. The number of super special code words in 
C\ is 2". 

We illustrate a super special row repetition code by the 
following example. 

Example 4.2: Let C s = [cj | C s 2 | C 2 | C* J be a super row 

repetition code with associated super special row matrix H s = 

[H| | H 2 1 H 3 1 H 4 ] 

"l 1 0 0 0 1 1 0 0 0 

10 10 0 10 10 0 
1 0 0 1 0 1 0 0 1 0 

1 0 0 0 1 1 0 0 0 1 



1 


1 


0 


0 


0 


1 


1 


0 


0 


o" 


1 


0 


1 


0 


0 


1 


0 


1 


0 


0 


1 


0 


0 


1 


0 


1 


0 


0 


1 


0 


1 


0 


0 


0 


1 


1 


0 


0 


0 


1 



DEFINITION 4.20: Let C, be a super special parity check mixed 
row code i.e., C s =[C S ; | C 2 S \ ... \ where C s is obtained 
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using the super special mixed row matrix H s = [Hj \ H 2 \ ... j H n ] 
where each Hi is a unit row vector having t, number of elements 
i.e., 



~1 1 ■■■ 1 


1 1 1 




1 1 f 


tj times 


t 2 times 




t n times 




where at least one t, & tj for i & j. Any super special code word in 
C s would be of the form 



r ; r 


1 


2 2 


2 


n n n 




1 


2 


n 


i 


..x, t 


Xj x 2 .. 


~x t 

l 2 


... XjX 2 ...x t 




X s 


X s 


... x s 



with H x T = (0); i.e., each x' s would contain only even number 
of ones and the rest are zeros. 



C s = [Ci | C -2 | ... | C n ] is defined to be super special parity 
check row code. C s is obtained from the parity check row matrix 
/ vector H s = [Hi | H 2 | ... H n ] where Hi = H 2 = ... = H n 



= 1 1 ••• 1 



m times 



Here a super special codeword in C s would be 



a super row vector of the form [x^ | x 2 s \ ... | x”] with each 
x' s = [x[ x' 2 ... x* B ] where only even number of x- are ones and 
the rest zero, 1 < j < m and i = 1,2, . . ., n. 



Now we will illustrate the two types of super special parity 
check (mixed) row codes. 



Example 4.3: Let C s = 



c[ 



c; c 



be a super special parity 



check code having H s = [Hi | H 2 1 H 3 ] = [11111|111|1111 
1] to be the super special parity check matrix associated with it. 



Definition 4.21: Let C s = [c' | C 2 S \ ... \ C"~\ be a super 
special row code. Suppose x s = [xj \ x] \ . . . \ x" ] is a 
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transmitted super code word and y s = [ y| | y 2 s \ . . . \ y " ] 
the received supercode word then e s = y s - x„ 



is 



f 1 1 I 2 2 I In nl f 1 \ 2 \ 

|_W -x,. | y, -x, | ... | y s -x s ] = le, | | ... 

the super error word or the super error vector. 



'] is called 



We first illustrate how the super error is determined. 

Example 4.4: Let C s = [c* | C s 2 | C 2 | C s 4 ] be a super special 

code with associated super parity check row matrix H s = [Hi | H 2 
|H 3 |H 4 ] 



"l 


0 


0 


1 


0 


0 


1 


0 


1 


0 


1 


0 


0 


1 


1 


0 


0 


1 


0 


0 


1 


0 


1 


0 


1 


0 


0 


1 


0 


0 


0 


1 


1 


1 


0 


0 


0 


0 


1 



10 10 0 
110 10 
0 10 0 1 



10 10 0 10 0 
110 10 0 10 
0 110 10 0 1 



Clearly y s + e s = x s . 

DEFINITION 4.22: The super Hamming distance d s (x s , y s ) 
between two super row vectors of the super special vector space 
Vs, where x s = [x' s \ x 2 s | ... | x"]and y s = | y ] \ ... \ y"] 

is the number of coordinates in which x' s and y' s differ for i = 
1, 2, ..., n. The super Hamming weight w/xj of the super vector 
x s = [x s ; | x 2 s | ... | jc"] in V s is the number of non zero 

coordinates in each x' s ; i = 1, 2, ..., n. In short w s (x s ) = d(x s , 

( 0 )). 



As in case of usual linear codes we define super minimum 
distance d s mm of a super special linear row code C s as 
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Cin= min r d s (ii s ,v s ), 

u s ,v s eC s 

U s^ v s 

d s (u s , v s ) = d s (u s - v s , (0)) = w s (u s - v s ). 

Thus the super minimum distance of C s is equal to the least 
super weights of all non zero super special code words. 

Now the value of 



d s mm = min d s (u s , v s ) 

u s ,v s e C s 



= mm 



» 1 ^c. d -([ u il u -l"-| u »]*[ vl -l v >|-"l v ”]) 

u s^ v s 

[a(u‘, v;) + d(q, v;) + ... + d(u,\ v;)]. 



Now dU n of the super special row code given in example 
3.1.2 is 7 verified using the fact in C s = [c' C 2 | Cf ] ; d' nin C^ 
= 3, dUC 2 = 2 and dNC 2 = 2 . Hence dU (C s ) = 3 + 2 + 2 = 
7. So we will denote d s mm =min d s (u„v s ) by dN n (C s ), u s , v s 
e C s , u s * v s . 



mm 

hyi^c 

x\*y\ 



[c; 

d(x 



c 



t) 



y + 



Q]-dL[(c;) + (c;) + ...+(c:)] 

d(x s 2 , y 2 ) + ...+ min d(x 

x " 



mm 

i <yl e c; 
* 2 ,*y; ' 



:)• 



Now we proceed on to define the dual of a super special 
row code. 

Definition 4.23: Let C s = [c' I C 2 I ... I C"] be a super- 

special row [(tii, ..., n„), (kj, ..., k n )] binary code. The super- 
special dual row code of C s denoted by 
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where ( C\ ) = {u\ \ u‘ s -v' s -0 for all v' &C' S }, i = 1, 2, ..., n. 

Since in C s we have nj - kj = n 2 - k 2 = ...= n„ - k n i.e., the 
number of check symbols of each and every code in C\ is the 
same for i = 1, 2, ..., n. Thus we see n = 2k, alone can give us a 
dual in all other cases we will have problem with the 
compatibility for the simple reason the dual code of C‘ being 
the orthogonal complement will have n, - k t to be the dimension, 
where as C' will be of dimension k t , i = 1, 2, .... n. Hence we 
can say the super special dual code would be defined if and only 
if «,■ = 2k i and such that m = n 2 = ... = n„. 



We can define the new notion of super special syndrome to 
super code words of a super special row code which is 
analogous to syndrome of the usual codes. 

DEFINITION 4.24: Let C s be a super special row code. Let H s be 
the associated super special parity check matrix of C, the super 
special syndrome of any element y s £ V s where C s is a super 
special subspace of the super special vector space V s is given by 
S(y s ) = H s y T s . S(y s ) = (0) if and only ify s £ C s . 

Thus this gives us a condition to find out whether the 
received super code word is a right message or not. Suppose y s 
is the received super special code word, we find S(y s ) = H s y T s ; 
if S(y s ) = (0) then we accept y s as the correct message if S(y s ) = 
H s y T s ^ (0) then we can declare the received word has error. 

We can find the correct word by the following method. 
Before we give this method we illustrate how the super special 
syndrome is calculated. 

Example 4.5: Let C s =[c* | C s 2 | C 2 C^] be a super special 

row code. Let H s = [Lf | H 2 | H 3 1 H 4 ] be the super special parity 
check matrix of C s . 
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H s = 
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0 


0 


1 





Now we have to shown how to find whether the received 
super special code word is correct or otherwise. It is important 
to note that what ever be the super special code word x s e C s 
(i.e., it may be a super special mixed row vector or not) but the 
syndrome S(x s ) = H s xJ is always a super special row vector 
which is not mixed and each row vector is of length equal to the 
number of rows of H s . 

Now we know that every super special row code C s is a 
subgroup of the super special vector space V s over Z 2 = {0, 1}. 
Now we can for any x s e V s define super special cosets as 

x s + c s = {x s + c s |c s G CJ. 



Thus 



V s = {Z 2 x Z 2 x . . .x Z 2 1 Z 2 x ... x Z 2 1 ... | Z 2 x Z 2 x ... x Z 2 } 
= C s u[x‘ + C s ] u ... u [x* +C S ] 

where C s =[c; \C 2 t | ... | C“] 



and 



x s 



2 


2 


n 


n 1 


X, . 


•X n . 


.. X, . 


,x„ 


1 


n 2 


1 


n " J 



1 2 


n“| 


x c x c .. 


• X 


s s 


s J 



and 
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x s C s 





Now we can find the coset leader of every x‘ +C‘ as in 
case of usual codes described in chapter one of this book. Now 
if 



y*= [y 1 , I y s 2 I ••• I y" ] 

is the received message and |^e‘ + (0) | e 2 + (0) |. . .| e” + (0)] is a 

special super coset leaders then using the relation y s - e s we get 
y s - e s to be super special corrected code word. It is interesting 
to note that each e‘ + (0) has a stipulated number of coset 
leaders depending on n ; , i = 1,2, . . ., n. 

We will illustrate this by the following example. 



Example 4.6: Let C s 




be a super special row code. 



Suppose H s = [Hi | H 2 ] be the super special row matrix 
associated with C s . Let 



H s =[Hi|H 2 ] 



1 0 
1 1 



Now c s =[c.:|c s 2 ] with C\ 



10 10 110 
0 10 110 1 ' 

= {(0 0 0 0), (1 0 1 1), (0 1 0 1), 



(1 1 1 0)}and C 2 = {( 0 0 0 0 0), (1 0 0 1 0), (0 1 0 0 1), (0 0 1 1 

1 ), (1 1 0 1 1 ), (0 1 1 1 0 ), (1 0 1 0 1 ), (1110 0 )}. 



Cg = {[0 0 0 0 0 0 0 0 0], [1 0 1 1 0 0 0 0 0], 



[0 1 0 1 I 0 0 0 0 0 ], [1 1 1 0 
[000 0|0100 1 ], [1 1 10 
[1 0 1 1 I 1 0 00 1 ], [00 00 

[0 1 0 1 I 00 1 1 1 ], [1 0 1 1 

[1 1 1 0 I 1 1 0 1 1 ], [1 0 1 1 

[0 0 0 0 I 0 1 1 1 0 ], [1 0 1 1 

[000 01 1 0 1 0 1 ], [1 0 1 1 
[0 0 0 0 I 1 1 1 0 0 ], [1 1 1 0 



0 0 0 0 0], [0 0 0 0 I 1 0 0 1 0], 

0 1 0 0 1], [0 1 0 1 I 0 1 0 0 1], 

0 0 1 1 1], [1 1 1 0 I 0 0 1 1 1], 

001 1 1], [00 001 1 1 0 1 1], 

1 1 0 1 1], [0 1 0 1 I 1 1 0 1 1], 

0 1 1 1 0], [0 1 0 1 I 0 1 1 1 0], 

1 0 1 0 1], [0 1 0 1 I 1 0 1 0 1], 

1 1 1 00], [1 0 1 1 I 1 1 1 00], 
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[0 1 0 1 | 1 1 1 0 0] and so on}. 

Clearly |C S | = 32 . Now the coset table of C' is given by 



Message 




code 










words 


0 0 


1 0 




0 


1 








1 1 


0 0 0 0 


10 11 




0 


1 


0 


1 




1 1 1 


Other cosets 


















10 0 0 


0 0 11 


1 


1 


0 


1 


0 


1 


1 0 


0 10 0 


1111 


0 


0 


0 


1 


1 


0 


1 0 . 


0 0 10 


10 0 1 


0 


1 


1 


1 


1 


1 


0 0 


coset leaders 



















Now the coset table of C s 2 is given by 



message 


0 


0 


0 






1 


0 


0 






codewords 


0 


0 


0 


0 


0 


1 


0 


0 


1 


0 




1 


0 


0 


0 


0 


0 


0 


0 


1 


0 


other cosets < 


0 


1 


0 


0 


0 


1 


1 


0 


1 


0 




0 


0 


1 


0 


0 


1 


0 


1 


1 


0 


message 


0 


1 


0 






0 


0 


1 






codewords 


0 


1 


0 


0 


1 


0 


0 


1 


1 


1 




1 


1 


0 


0 


1 


1 


0 


1 


1 


1 


other cos ets < 


1 


0 


0 


0 


1 


0 


1 


1 


1 


1 




0 


1 


1 


0 


1 


0 


0 


0 


1 


1 
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message 


1 


1 


0 






0 


1 


1 






codewords 


1 


1 


0 


1 


1 


0 


1 


1 


1 


0 




0 


1 


0 


1 


1 


1 


1 


1 


1 


0 


other cosets < 


1 


0 


0 


1 


1 


0 


0 


1 


1 


0 




1 


1 


1 


1 


1 


0 


1 


0 


1 


0 


message 


1 


0 


1 






1 


1 


1 






codewords 


1 


0 


1 


0 


1 


1 


1 


1 


0 


0 




0 


0 


1 


0 


1 


0 


1 


1 


0 


0 


other cosets < 


1 


1 


1 


0 


1 


1 


0 


1 


0 


0 




1 


0 


0 


0 


1 


1 


1 


0 


0 


0 


Suppose y s 


= 


[1 1 


1 


1 1 1 


1 1 1 1] is 


the 


: received word then 



S(y s ) = U s y T s * [(0) | (0)]. e s = [0 1 0 0 | 0 0 1 0 0] is the super 
set coset leader. Thus x s = y s + e s = [1 0 1 1 | 1 10 1 1] e C s . 

With the advent of computers calculating the super special 
coset leaders is not a very tedious job. Appropriate programs 
will yield the result in no time. 

Now we proceed on to describe/defme the super special row 
cyclic code. 

DEFINITION 4.25: Let C s = [Ci \ C 2 | ... I CJ be a super special 
row code. If every’ C' is a cyclic code in C s we call C s to be a 
super special cyclic row code. H s = [Hj \H 2 \ ... \H„] denotes the 
super special parity check row matrix of the super special cyclic 
code. 

We illustrate this by the following example. 

Example 4.7: Let C s =[c* | | Cf ] be a super special cyclic 

code with an associated super special parity check matrix 



H s = [H! | H 2 1 H 3 ] 
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110 10 0 
0 110 10 
0 0 110 1 

001111001001 
010110010010 
101100100100 

We see each of Cf , and are cyclic codes. 

Now we see in general for any super special mixed row 
code with an associated super special parity check matrix H s 
which happens to be a super mixed row matrix we cannot define 
the super special generator row matrix G s . The simple reason 
being if the code words in each of the C‘ in C s where 
C s =[c‘ | Cj | ... | C“], i = 1, 2, n happens to be of 

different length then it would be impossible to define a super 
generator matrix. So we shall first define the notion of super 
special generator row matrix of a super special row code(mixed 
row code). 

Definition 4.26: Let C = [c' C] | ... | C" ] be a super 
special row code. A super special row matrix which generates 
C s exists if and only if in each C' the codes in C s have the same 
number of message symbols, that is if C, has a super special 
parity check row matrix H s = [Hi \ H 2 | ... | H n ] then we 
demanded each C' must have the same number of check 
symbols. Likewise for the super special generator row matrix to 
exist we must have G s = [G/ \ G 2 \ ... \ G n ] where C' have the 
same number of message symbols which forms the number of 
rows of the super row generator matrix G s . 

We shall first illustrate this by an example. 
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Example 4.8: C s = [c* | C s 2 ] be a super special row code. 
Let G.=[g: |Gf] 



"l 


1 


0 


1 


0 


0 


1 


0 


0 


0 


1 


o" 


0 


1 


1 


0 


1 


0 


0 


1 


0 


0 


1 


1 


0 


0 


1 


1 


0 


1 


0 


0 


1 


0 


1 


1 



One can easily find the super code words using Gs [50]. 

Definition 4.27: Let C s = [c' | C; | ... | C”] be a super 

special mixed row code. If each of the codes C' have the same 
number of message symbols then we have the super special 
generator mixed row matrix G s = [g] | G 2 | ... | G”] 
associated with C s . Number of message symbols in each of the 
G' are equal and is the super special mixed row matrix G s ; 1 < 
i <n. 

We illustrate this by the following example. 

Example 4.9: Let C s = [c, | C s 2 | C 2 C 2 ] be a super special 
mixed row code. Let G s =[g 2 | G s 2 | G 2 ] be the associated 
super special mixed row generator matrix given by 

o.=[g; |g; | a;] 



10 0 110 
0 10 0 11 
0 0 110 1 
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1 


0 


0 


0 


0 


1 


0 


1 


0 


0 


1 


0 


0 


l" 


0 


1 


0 


1 


0 


0 


0 


0 


1 


0 


0 


1 


1 


0 


0 


0 


1 


0 


1 


0 


1 


0 


0 


1 


1 


0 


1 


0 



Theorem 4.4: Let C s = [c‘ | Cf 1 . . . | C” ] be a super special 

row code with H s = [Hj \ H 2 \ ...| H,J , the super special parity 
check matrix. If each H, = (A,-, I n -k), i = L 2, ..., n then G s = [Gi 
| G 2 | ...| G,J with Gi = (Ik, - A T ); 1 < i < n if and only if the 
length of each code word in C' is the same for i = 1, 2, ..., n. 

For proof refer [50]. 

For any super special row code we can FIs the parity check 
matrix and the super special generator matrix G s such that 
G s Fig = (zero super vector). 

Now having defined the new class of super special row 
(mixed row) codes we will now define new classes of mixed 
super classes of mixed super special row codes C s i.e., we may 
have the super special row code to contain classical subcodes as 
Flamming code or cyclic code or code and its orthogonal 
complement and so on [50]. 

DEFINITION 4.28: Let C s = [Cl \ C 2 \ ... I C n ] be a super special 
row code. If some of the C, ’s are Hamming codes, some C, ’s are 
cyclic codes i ^j, some Cfs are repetition codes and some C t ’s 
are codes and C p ’s are dual codes of C, ’s ; 1 < j, k, t, i, p < n 
then we call C s to be a mixed super special row code. 

It is important to mention here that even if two types of 
classical codes are present still we call C s as a mixed super 
special row code. 

We will illustrate them by the following examples. 
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Example 4.10: Let C s = [Ci | C 2 | C 3 | C 4 ] be a mixed super 
special row code. Here Ci is a Hamming code, C 2 the repetition 
code, C 3 a code of no specific type and C 4 a cyclic code. 

Let the mixed super special parity check matrix H s 
associated with C s be given by 





H s = 


[Hi 


1 H 2 1 


H 3 |H 4 ] 






"0 


0 


0 


1 1 1 


1 


1 


1 


0 


0 


0 


1 


1 


0 0 1 


1 


1 


0 


1 


0 


1 


0 


1 


0 1 0 


1 


1 


0 


0 


1 


1 1 


1 


0 


0 


1 


1 


1 


0 


1 


0 


0 1 


0 


1 


0 


0 


1 


1 


1 


0 


1 


1 0 


0 


0 


1 


0 


0 


1 


1 


0 


0 



Example 4.11: Let C s = [Ci | C 2 | C 3 ] be a mixed super special 
row code. Let H s = [Hi | H 2 | H 3 ] be the associated super special 
parity check mixed row matrix. Ci is the Hamming code, C 3 any 
code and C 3 a repetition code. 



11111111 
0 0 0 1 1 1 1 0 
0 110 0 110 
10 10 10 10 



1 0 110 0 0 
1 110 10 0 
1 0 0 0 0 1 0 
0 1 1 0 0 0 1 



110 0 0 
10 10 0 
10 0 10 
1 0 0 0 1 



is the mixed super special parity check mixed row matrix for 
which G s does not exist. 
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We define the new notion of super special Hamming row 
code. 

Definition 4.29: Let Cs = \_c’ s \ C 2 S \ ... | C" ] where each 

C' s is a (2 m - 1, 2 m - 1 - m) Hamming code for i = 1, 2, ..., n. 

Then we call C s to be a super special Hamming row code. If 
H s = [Hi | H 2 | ... j H,J be the super special parity check matrix 
associated with C s we see H s is a super special row matrix 
having m rows and each parity check matrix H, has m rows and 
2"' - 1 columns, i =1, 2, ..., n. 

Further the transmission rate can never be equal to l A. If 
m > 2 then will the transmission rate be always greater than Vi ? 

For more examples refer [50]. 

We define using super column matrices the notion of super 
column codes which is an innovative means of using super 
column matrices. 

DEFINITION 4.30: Suppose we have to describe n codes each of 
same length say m but with varying sets of check symbols by a 
single matrix. Then we define it using super column matrix as 
the super code parity check matrix. Let 



C, 




be a set of m codes, Cj, C 2 , ..., C m where all of them have the 
same length n but have n - kj, n - /o, ..., n - k m to be the number 
of check symbols and k h k 2 , ..., k m are the number of message 
symbols associated with each of the codes Ci, C 2 , ..., C m 
respectively. 



Let us consider 




Supermatrices and their Applications 



143 



H 1 



H 



where each H 1 is the n - k; x n parity check matrix of the code 
Q; i = 1, 2, m. We call H s to be the super special parity check 
mixed column matrix of C s and C s is defined as the super 
special mixed column code. 

The main difference between the super special row code 
and the super special column code is that in super special row 
codes always the number of check symbols in every code in C s 
is the same as the number of message symbols in Q and the 
length of the code C; can vary where as in the super special 
column code, we will always have the same length for every 
code Q in C s but the number of message symbols and the 
number check symbols for each and every code C, in C s need 
not be the same. In case if the number of check symbols in each 
and every code Q is the same. Then we call C s to be a super 
special column code. 

In case when we have varying number of check symbols 
then we call the code C s to be a super special mixed column 
code. 

In the case of super special column code C s = [Ci | C 2 | . . . j 
C m ]' we see every code Q in C s have the same number of 
message symbols. Thus every code is a (n, k) code. It may so 
happen that some of the Q and Cj are identical codes. 

For examples refer [50], 



Definition 4.31: Let C s = [C,\ C 2 | ... C„/' where Ci’s are 
codes of same length m. Suppose each C, is generated by a 
matrix G it i = 1, 2, ..., n, then 
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G s = 




G 



n 



generates the super special column code C s . We call G s the 
super special generator column matrix which generates C s . If in 
each of the codes C, in C s , we have same number of message 
symbols then we call G s to be a super special generator column 
matrix; i = 1, 2, ..., n. If each of the codes C, ’s in C s have 
different number of message symbols then we call G s to be a 
super special generator mixed column matrix. 

We say G s is in the standard form only if each G; is in the 
standard form. Further only when G s is in the standard form and 
G s is a super special column matrix which is not a mixed matrix 
we have H s the super special parity check column matrix of the 
same C s with 



0 




0 



For examples refer [50]. 

We proceed onto define classical super special column 
codes. 

DEFINITION 4.32: Let C s = [Cl \ C2 \ ... I C n ]' be a super special 
column code if each of the code C, is a repetition code of length 
n then C s is a super special repetition column code with C/ = 
C2 =...= C n . 
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The super special column parity check matrix 

H 

H s = T 

H 



1 1 0 0 ... 0 

1 0 1 0 ... 0 

1 0 0 0 ... 1 

1 1 0 0 ... 0 

1 0 1 0 ... 0 

1 0 0 0 ... 1 

where 

1 1 0 0 ... 0 

1 0 1 0 ... 0 

1 0 0 0 ... 1 

It is important to note that unlike the super special 
repetition row code which can have different lengths the super 
special repetition column code can have only a fixed length and 
any super special code word 

I I ••• I <] 

where x J s =(11 ... 1), n-times or (0 0 ... 0) n-times only; 
1 <j < n. 
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DEFINITION 4.33: Let C s = [Ci\ C 2 | ... I C,f be a super special 
parity check column code . Let the super special parity check 
column matrix associated with C s be given by 



Where 



H, 




H 



n 



H, = H 2 = ... =H n =(l 1 ... 1 ). 

m— times 



Thus we see we cannot get different lengths of parity check 
codes using the super special column code. However using 
super special row code we can get super special parity check 
codes of different lengths. 

DEFINITION 4.34: Let C s = [Cl \ C 2 | ... I C,f be a super 
special column code if each of the codes C, in C s is a (2 m - 1, 2 m 
- 1 - m) Hamming code for i = 1, 2, ... , n then we call C s to be 
a super special column Hamming code . It is pertinent to 
mention that each code C , is a Hamming code of same length; i 
= 1 , 2, ... , n. 

DEFINITION 4.35: C s = [Cl \ C 2 \ ...| C„]' is a mixed super 
special column code if some C, ’s are repetition codes of length 
n some Cj’s are Hamming codes of length n, some C k ’s parity 
check codes of length n and others are arbitrary codes, 1 < i, j, 
k <n. 

We illustrate this by the following example. 

DEFINITION 4.36: Let C s = [Cl \ C 2 \ ... \ C m f be a super- 
special column code if each of the codes C, is a cyclic code then 
we call C s to be a super special cyclic column code. However 
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the length of each code C,; i = 1, 2, ..., n will only be a cyclic 
code of length n, but the number of message symbols and check 
symbols can be anything. 



Now we illustrate this by an example. 

Definition 4.37: Let 



C(S) = 



cl 


c 1 




c 1 


Cl 


c 2 




C 2 

n 










c; 


s~im 

^2 




s~im 



where C' are codes 1 <i < m and 1 <j <n. Further all codes 
C'j.C^ , are of same length C' 2 ,C 2 2 , C 2 are of same 

length and C' n ,C 2 n , ..., C 2 are of same length. C' ,C 2 ,...,C] t 
have same number of check symbols, C 2 ,C 2 , ..., C 2 n have same 
number of check symbols and C" 1 ,C"‘ ,...,C'" have same number 
of check symbols. 



We call C(S) to be a super special code. We can have the 
super parity check matrix 



H(S) = 



Hi 


Hi 




Hi, 


Hi 


Hi 




Hi 










H n ; 


h; 




h : 



where H‘ ’s are parity check matrices 1 <i <m and 1 <j < n. 
Further H ^ ,H H l have the same number of rows, 

H] ,H 2 2 , ...,H 2 have same number of rows and so on. 
H'fH'", have the same number of rows. Likewise 

H 2 ,H] have the same number of columns, 
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H[,H l H™ have the same number of columns and so on. 

H ! n ,Hf..., H'" have same number of columns. 

The notion of super matrices are ingeneously applied to 
super fuzzy models, super linear algebra and super special 
codes. Such type of applications is not only innovative but can 
lead to easy comparison and save the working time. 




Chapter Five 



Interval Matrices and Natural 
Gass of Intervals 



In this chapter first we recall properties of intervals and then 
define the notion of matrix interval using matrices with entries 
from R. Also we define the matrix of natural class of intervals 
and illustrate them with examples we however get the relation 
between the interval matrix and matrix interval and the relation 
between the natural class of matrix intervals and the matrix of 
natural class of intervals. We will call [A. B] an interval where 
A > B or A < B or A = B or A and B are not even comparable 
but are bound by some common features. When A and B are 
integers or modulo integers we get intervals which we call them 
as natural class of intervals. For if [15, 20] where 15 and 20 e 
Z 2 5 is the situation when we cannot in any way compare them. 
So if we replace A and B by matrices of same order we call that 
interval as matrix interval. 

With this rough idea now we proceed onto define the row 
matrix interval. 
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By an interval [a, b] we mean all elements x such that a < x 
< b are in [a b]. [a a] is nothing but a. So when in an interval 
a = b it collapses to a single point. 

LetQi= {[a, b] I a, b e Q}. 

Let {[a, b] I a, b e Q + u {0}, Q + denotes the set of 
positive rationals } . 

Zi = { [a, b] I a, b eZj; denotes the set of all intervals on Z, 
i.e., the set of positive and negative integers. 

Z\- {[a. b] I a, b e Z + u {0}; Z + the set of positive 
integers } . 

Ri = {[a, b] I a, b e R} i.e., the collection of all intervals in 
the set of reals. 

R} = {[a, b] I a, b e R + u {0}, R + the set of positive reals}. 

Zj, = { [a, b] I a, b eZJ. 

Thus we have defined the collection of all intervals of 
different types. We first make the following observation which 
is very essential. 

We see Q c: Qi (a = [a, b] e Qj). 

Likewise Z cz Z h Q 1 ci , Z n c Z\ , Rj ^ R, R + c R} . 

Thus these which will be known as interval sets are 
generalized or contain the related sets as subsets. We shall call 
Qi the rational intervals, Z| the integer intervals, R, the real 
intervals and Z l n the modulo integer intervals. 



Further Z c; Z[ c Qi cr R[ 
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ZcQcQicR^ZcQcRcR]. 

[-7, 3] is an integer interval in Z\. 

[3, 5] is a positive integer interval in Z\ . 

We see Z + c Z\ c Z^ 

[-5/2, 3/7] is a rational interval in Qi. 

[9/7,5/2] is a positive rational interval in Q+ ; Q + c Q+ c 

[-a/2, a/5] is a real interval in Rj but [a/3, a/7 ] is a 
positive real interval in Rj\ Thus R + c R, c R[. We see [0, 1] 
[0, 4] and modulo integer intervals of Z\ . 

It is important and interesting to note Z\ has only finite 
number elements in it. 

Infact Zj has only four elements viz. {[0, 0], [1, 1] [0, 1]}. 

Z\ has only 6 elements viz {[0, 0] [1, 1], [2, 2], [0, 1], [0, 

2 ], [ 1 , 2 ]}. 

Z\ has only 10 intervals viz. {[0, 0], [1, 1], [2, 2], [3,3]. [0, 
1], [0, 2], [0, 3], [1,2], [1,3], [2, 3]}. 

Z\ has only 15 modulo intervals. {[0, 0], [1, 1], [2, 2], [3, 
3], [4, 4], [0, 1], [0, 2], [0, 3], [0, 4], [1, 2], [1, 3], [1, 4], [2, 3], 
[2,4], [3,4]}. 



Z\ has only 28 modulo intervals. 




152 



Innovative Uses of Matrices 



Thus Z* has only n+ (n-1) + ... +2+1 = n ^ n ^ + ^ modulo 
intervals. 

Z 15 has only {[0. 0], [1. 1], [2, 2], ..., [14. 14]. [0, 1]. [0, 2], 
.... [0. 14], [1. 2], ..., [1. 14]. [2, 3], ..., [2. 14], .... [13, 14]} - 
Z 1 

+5 • 

Number of intervals in Z [ 5 is 120. 

Now we recall the operations on intervals given by 126] 

Definitions [26] 

[a, b] = [c, d] if and only if a = c and b = d 
[a, b] c [c, d] if only if c < a < b < d. 

[a, b] < [c, d] if and only if b < c. The width of an interval 
[a, b] is denoted by w([a, b]) = b - a and the magnitute of the 
interval l[a, b]l = max (lal, Ibl). 

Degenerate intervals are those intervals which has zero 
width, can be identified with real numbers. Thus a = [a, a]. 

We now proceed onto give the arithmetic operations as 
given by Kuperman. 

If * denotes any one of the symbols +. -, ., I (. Product, / 
division). 



Then the interval [a, b] * [c. d] = 



x * y 



a < x < b 
c< y <d 



except 



we do not define [a.b] I [c, d] if o e [c, d]. (This excludes the 
definition by zero possibility). Alternatively we can define the 
arithmetic operations on intervals by giving the endpoints of the 
intervals resulting from the sum, difference, product or quotient 
of two intervals. 

Equivalent to the definition in I we thus we 
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[a, b] + [c, d] = [a+c, b+d] 

[a, b] - [c, d] = [a-d, b-c] 

[a, b] . [c, d] = [min (ac, ad, be, bd), max (ac, ad, be, bd)] 

[a, b] / [c, d] = [a, b] / [1/d, 1/c], provided 0 g [c, d]. 

For the denominator cannot be zero. Thus in the definition 
of division c and d are either both positive or both negative. By 
way of numerical example we thus have 

(i) [-2, 3] + [-9,-3] = [-11,0] 

[8, 11] + [5,7] = [13, 18]. 

[-3, 8] + [3, 4] = [0, 12] 

(ii) [-3, 8] - [-7, -2] = [-3 - (-2), 8+2] 

= [- 1 , 10 ] 

[8, 17] - [-2, 6] = [8 - 6, 17 - (-2)] 

= [2, 19] 

[2, 7] - [-1, 3] = [2-3, 7— (— 1 )] 

= [- 1 , 8 ] 

(iii) [3, 5] . [7, 9] = [min {21, 27, 35, 45], 

max {21, 27, 35, 45}] 

= [21, 45] 

[-2,1]. [2,5] = [min {-4, 2, 5, -10], 

max {-4, 2, 5, -10}] 

= [-10, 5] 

[-5, 2] . [3, 7] = [min {-15, 6, 14, -35], 

max {-15, 6, 14, -35}] 

= [-35, 14] 

[-2, 7] / [3, 8] = [min {-2/3, -2/8, 7/3, 7/8], 

max {-2/3, -2/8, 7/8, 7/3}] 

= [-2/3, 7/3] 
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[-5, 4] / [-2,6] = [min [5/2, -5/6, -4/2, 4/6], 

max [5/2, -2, 4/6, -5/6}] 

= 1-2, 5/2]} 

[-8, 9] / [3, 7] = [min [-8/3, 8/7, 3, 9/7], 

max [-8/3, 8/7, 3, 9/7} 

= [-8/3, 3] 

Now we recall the properties of interval arithmetic given by 
Kooperman. 

For more about these notions please refer [26] 

We may note that interval arithmetic is a generalization or 
extension of real arithmetic since [a, a] is a real number and the 
definitions hold for intervals of this form as well. 

Thus laws that hold for interval arithmetic must 
automatically hold for real arithmetic. But laws that hold for 
real arithmetic may not always hold for interval arithmetic. 

For example, interval arithmetic is associative and 
commutative with respect to addition and multiplication. Thus, 
given interval numbers I, J, K we have 

(I, J) + K = I + (J + K) 

1+ J= J+I 

But distributive laws does not always hold for interval 
arithmetic. 

[3, 5] ([1,2] + [-7, 2]) 

= [3,5] [1,2] + [3,5] [-7, 2] 

= [min [3, 5, 6, 10], max [3, 5, 6, 10}] + 

[min [-21, 6, -35, 10], max [-21, 6, -35, 10}] 

= [3, 10] + [-35, 10] 

= [-32, 20] = [3, 5] [1, 2] + [3, 5] [-7, 2] 
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Consider [3, 5] ([1, 2] + [-7, 2]) 

= [3, 5] [-6, 4] = [min {-18, -30, 12, 20}, 

max {-18, -30, 12, 20}] 

= [-30, 20] 

Thus [3, 5] ([1, 2] + [-7, 2]) * [3, 5] [1, 2] + [3, 5] [-7, 2], 

It is pertinent to mention here that the distributive laws hold 
good when the interval is a real number say k ; i.e., the 
degenerate interval [k, k] of zero with, then it is easy to verify 

[k, k] ([a, b] + [c, d]) = [k, k] [a, b] + [k, k] [c, d] 

= k([a, b] + [c, d]) = k [a, b] + k [c, d] 

= k [a+c, b+d] = [ka, kb] + [kc, kd] 

= [k(a+c), k(b+c)] = [ka + kc, kb + kd] 

= [ka + kc, kb + kd] 

Thus we have 

b [a - 8a, a + Sa] = b ([a, a] + [-8a, 8a]) 

= b[a, a] + b[-Sa, 8a] 

[ba, ba] + [— Ibl 8 a, Ibl 8a] 

= [ba - IbISa, ba + Ibl 8a] 
b[a - 8a, a + 8a] = [ba - Ibl 8a, ba + IbISa] 

(We may note that 8a is non negative by implication of its 
appearance in [a - 8a, a + 8a] for in any interval [c, d], c <d). 

We now proceed onto recall the concept of rounded interval 
arithmetic. An important use of interval arithmetic is to 
determine upper bounds for the errors due to rounding. 

In rounded interval arithmetic working in either single or 
double precision, we regard each real number as an interval. An 
in the first place we round the lower endpoint down and the 
upper end point up. 
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In this way, the maximum possible effect of rounding is 
taken into account, so that the interval for the final result 
contains the true value. 

In using rounded interval arithmetic it should be clear that 
the final result can be obtained as sharply as desired, i.e., the 
width of the final interval can be as small as desired, by using a 
sufficiently long wordlength. By way of a numerical example. 
Suppose we multiply two number, a = 0.2310 I 581 and b = 
0.8351 I 621 is a decimal coded machine using a precision of 
wordlength of four decimal digits. 

Then, in the first place we write down the intervals 
containing a and b. 

a e [0.2310, 0.2311] and be [0.8351,0.8352] 

Then carrying out the multiplication we have 

ab e [0.2310, 0.2311] [0.8351, 0.8352] 

= [0.1929 118,..., 0.1929192...] 

= [0.1929, 0.1930], 

Thus, rounded interval arithmetic leads to the following 
interval which contains ab : [0.1929, 0.1930], the width of this 
interval being 0.1930 - 0.1929 = 0.0001. 

To obtain the product ab more sharply we clearly have to 
use a longer wordlength for the computation. 

The actual rounding down and up can most easily be 
achieved in machine coding. The lower endpoint should be 
rounded down so that the rounded value is less than or equal to 
the value being rounded. Thus if the value being rounded down 
is positive we need merely truncate while if the value is 
negative we truncate, examine the truncated portion, and, if this 
is not zero, we subtract one unit in the least significant position. 
And the upper end point must be rounded up if a computation is 
being carried out in rounded interval arithmetic, an interval 
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being said to be properly rounded if the endpoints are rounded 
as described above. 

We may mention that the rounding up and down can be 
coded in Fortan necessary although not as efficiently. 

Thus the final interval obtained by rounded interval 
arithmetic takes account of the errors due to rounding and 
contains the true value. 

Interval vectors and matrices and vectors and matrices 
whose elements are interval numbers, the superscript I being 
used to indicate such a vector or matrix. 

Thus given matrices B = (by) and C = (cy) of order n such 

that 

by < cy, i, j = 1,2, ..., n, then the interval matrix A*= [B, C] 
is defined by 

A 1 = [B, C] = {A = (ay) I by < ay < Cy; i, j = 1, 2, ..., n}. 
Kuperman [26] calls B the lower endpoint matrix of the interval 
matrix A 1 = [B, C], C the upper endpoint matrix and M = 
Vi (B+C) the midpoint matrix. 

DEFINITION 5.1: Let A = (aj, ..., a„) and B = (bj, ..., b„) a i} bi 
e R and Z n or Q or Z or C where A and B are of some order, 
viz; 1 x n. We define the row matrix of natural class of interval 
as 

[A, B] - [(a,, a„), (b h .... b n ) ]. 

If each a; > b ; , 1 < i < n we say [A, B] such that A > B if 
each aj < bj, 1 < i < n we say the row matrix natural class of 
interval [A, B] with A < B. For more about properties and 
working with natural class of intervals refer [57-58]. 

If some a ; > b and some aj < b, and some a k = b k we cannot 
say A > B or A < B. We will illustrate the row matrix of natural 
class of interval or in short row matrix interval by some 
examples. 
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Example 5.1: Let X = [(5, 3, -1, 0, 2, 4), (2, 4, 5, 7, -3, 1)] = 
[A, B] is the row matrix interval clearly A is not comparable 
with B. 

Example 5.2: Let Y = [A, B] = [(5, 2, 0, 7, 9, 6, 10), (7, 3, 4, 5, 

9, 10, 8, 12)] be a row matrix interval where A < B. 

Example 5.3: Let M = [(8, 3, 4, 9, -2, 7, 18, -9), (2, 1, 3, 5, - 

10, 2, 9, -14) = [A, B] be the row matrix interval where A > B. 

Example 5.4: Let S = [(3, 4, 2, 1, 5), (3, 4, 2, 1, 5)] = [A, B] be 
the row matrix interval, here A = B. So S = (3, 4, 2, 1, 5). We 
call this row matrix interval as degenerate interval. 

We will now proceed onto define the notion of addition, 
multiplication, subtraction and division of these row interval 
matrices. 

First we define the notion of order in a row interval matrix. 

Let X = [A, B] be a row interval matrix where order of A 
and B is 1 x n. We define the order of X, the row interval 
matrix to be 1 x n. Thus X is a 1 x n row interval matrix. 

We can add two row interval matrices if and only if they are 
of same order otherwise addition is not defined; for we cannot 
add a 1 x 5 row matrix with a 1 x 9 row matrix. 

Definition 5.2: Let X = [A, B] - l (a,, a 2 , ..., a n ), (b h b 2 , ..., 
b„)] and Y ~ [C, D] - [(a, c 2 , ..., c„), (di, d 2 , ..., d„) ] be any two 
1 xn row matrix interval. We define X + Y = [A. B] + [C, D] 

= [A + C,B + D] 

= [( a h ..., a n ) + (c h ..., c„), (bi, ..., b„) + (d h .... d„) ] 

= [(aj + a, ..., a n + c„), (b, + d h .... b„ + d n ) ] 

X + Y is again a 1 xn row matrix interval. 

We will first illustrate this by some examples. 
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Example 5.5: Let X = [A, B] = [(3, 1, 5, 0, -3, 2), (7, 8, 9, 11, 
0, 1)] and Y = [C, D] = [(-5, 2, 3, -7, 8, 10), (-2, 11,10, -4, 5, - 
3)] be any two 1x6 row matrix intervals. 

We see X + Y = [A, B] + [C, D] 

= [(3, 1,5,0, -3,2), (7,8,9, 11,0, 1)] + 

[(-5, 2, 3, -7, 8, 10), (-2, 11,10, -4, 5, -3)] 

= [(3, 1,5,0, -3, 2) + (-5, 2, 3, -7, 8, 10), (7, 8, 9, 1 1, 0, 1) 

+ (-2, 11, 10, -4, 5, -3)] 

= [(3 + (-5), 1+2, 5+3, 0+(-7), -3 + 8, 2+10), 

(7-2, 8+11, 9+10, 1 1+(— 4) 0+5, 1 + (-3)) 

- [(-2, 3, 8, -7, 5, 12), (5, 19, 19, 7, 5, -2)] 
is again a 1 x 6 row matrix interval. 

Example 5.6: Let P = [A, B] = [(3, 7, 1), (5, 1, -3)] and 

R = [C, D] = [(-2, 1, 0), (-7, 5, 2)] be two 1x3 row matrix 
intervals. 

P + R = [A, B] + [C, D] 

= [A + C, B + D] 

= [(3, 7, 1) + (-2, 1, 0), (5, 1, -3) + (-7, 5, 2)]. 

Now as in case of natural class of intervals we may have 
A + B>C + DorA + B<C + DorA + B not comparable with 
C + D or A + B = C + D. 

Example 5. 7: Let 

S = [(0, 0, 0, 0, 0, 0), (3, 1, 2, -5, 7, 2)= [S u S 2 ] and 
P = [(7, 2, 1, 0, 5, 1), (0, 0, 0, 0, 0, 0)] = [P ls P 2 ] be any two 
1x6 row matrix intervals. 

NowS + P = [S 1 ,S 2 ] + [P 1 ,P 2 ] 

= [(Si + Pi), (S 2 + P 2 )] 

= [(0, 0, 0, 0, 0, 0) + (7, 2, 1, 0, 5, 1), (3, 1, 2, -5, 7, 2) 

+ ( 0 , 0 , 0 , 0 , 0 , 0 )] 

= 1(7, 2, 1,0,5, 1), (3, 1,2, -5, 7, 2)] 



is again a 1 x 6 row matrix interval. 
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We observe that addition of 1 x n row matrix intervals is 
commutative for if A = [Ai, A 2 ] and B = [Bi, B 2 ] then A + B = 
B + A. 

Further addition of 1 x n row intervals is associative, that is 
if A = [Ai, A 2 ], B = [B,, B 2 ] and C = [C,, C 2 ] then (A + B) + C 
= A + (B + C). 

We define 0 = [(0, 0, 0), (0, 0, 0)] as the zero 

interval matrix or 1 x n row zero matrix interval or 1 x n row 
matrix zero interval. 0 acts as the additive identity. 

For if A = [Ai, A 2 ] and (0) = [(0). (0)] then A + (0) = 
(0) + A = A for all A; 1 x n row matrix interval. Further if 
A = [A], A 2 ] then -A = [-Ai, -A 2 ] and A + (-A) = [(0). (0)]. 

Let V = { [A. B] = M where A and B are 1 x n row matrices 
with entries from R or Z or Q or Z„ or C}; V is an additive 
abelian group of 1 x n row matrix interval. 

Now we can define product of two 1 x n row matrix 
intervals. Suppose A = [Ai, Bi] and B = [Ci, Di] be two 1 x n 
row matrix intervals; we define product AB as follows: 



Let AB = [A 1 ,B 1 ] [Ci.Dj] 

= [AiCi, BiDi] 

= [(ai, a 2 , ..., a n ), (c l5 c 2 , .... c n ) (bi, b 2 , ..., b n ), 

(di, d 2 , .... d n )] 

= [(ai, Ci. a 2 c 2 , .... a n c n ), (bidi. b 2 d 2 , .... b n d n )]. 

We see AB is again a 1 x n row matrix interval. 

We will illustrate this situation by some simple examples. 
Example 5.8: Let 

V = [A. B] = [(5, 7, -2, 0. 7, 8. 9. 11), (3. 1. +2. 4. 0. 18. 3, 0)] 
and S = [S 1; S 2 ] - [(3. 2, 0. 1. 5, 0, -2, 0). (0, 4. 0. 2, 1. 0. 2. -4)] 
be two 1x8 row matrix intervals. 
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Consider VS = [A, B] [Si, S 2 ] 

= [AS!,BS 2 ] 

= [(5, 7, -2, 0, 7, 8, 9, 1 1), (3, 1, +2, 4, 0, 18, 3, 0), 

(3, 2, 0, 1, 5, 0, -2, 0), (0, 4, 0, 2, 1, 0, 2, -4)] 

= [(15, 14, 0, 0, 35, 0, -18, 0), (0, 4, 0, 8, 0, 0, 6, 0)] 

is again a 1 x 8 row matrix interval. 

Example 5.9: Let V = [(3, 2, 0, 5, 1), (7, 8, 1, 0, 3)] and 
M = [(0, 1, 2, 0, 5), (3, 0, 0, 7, 0)] be two 1x5 row matrix 
intervals. Now VM = [(3, 2, 0, 5, 1), (7, 8, 1, 0, 3), (0, 1, 2, 0, 
5), (3, 0, 0, 7, 0)] 

= [(0, 2, 0, 0, 5), (21,0, 0, 0, 0)] is again a 1 x 5 row matrix 
interval. 

Example 5.10: Let A = [(-3, 2, 5, 1, 0, 2), (7, 0, 0, 8, 0, 1)] and 
B = [(0, 0, 0, 0, 8, 0), (0, 7, 5, 0, -3, 0)] be two 1x6 row matrix 
intervals. 



We see AB = [(-3, 2, 5, 1, 0, 2), (7, 0, 0, 8, 0, 1)] 

[(0, 0, 0, 0, 8, 0), (0, 7, 5, 0, -3, 0)] 

= [(-3, 2, 5, 1, 0, 2), (0, 0, 0, 0, 8, 0), (7, 0, 0, 8, 0, 1) 

(0, 7, 5, 0, -3, 0)] 

= [( 0 , 0 , 0 , 0 , 0 , 0 ), ( 0 , 0 , 0 , 0 , 0 , 0 )]. 

Thus A * [(0), (0)] and B * [(0), (0)], 
but AB = [(0), (0)]. 

Now we see if V = {M = [A, B] I collection of all 1 x n row 
matrix intervals} then V is closed with respect to multiplication. 
The operation multiplication is commutative and associative on 
V; but we see V is only a commutative semigroup under 
multiplication. For we see V has infinitely many zero divisors. 

Example 5.11: Let M = [A, (0)] and [(0), B] = P be two 1 x n 
row matrix intervals. We see MP = [(0), (0)] is a zero divisor. 



Now we define subtraction of two 1 x n row matrix interval. 
Suppose M = [(5, 3, -1, 0, 5), (3, 2, 4, 2, 0)] and P = [(3, 0, 1, 2, 
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4), (2, -1, 5, 4, -3)] be two 1x5 row matrix intervals. We fine 
M - P = [(5, 3, -1. 0. 5), (3. 2. 4, 2, 0)] [(3. 0. 1, 2, 4), (2, -1. 5, 
4, -3)] = [(2, 3, -2, -2, 1), (1, 3, -1. -2, 3)] is again a 1 x 5 row 
matrix intervals. 

The following points are pertinent to be observed. 

(1) The operation subtraction is not commutative. 

(2) The operation of subtraction is non associative. 

Consider A = [Ai, BJ = [(3. 2, -1, 0), (8. 5, 1, 2)] and 
B = [A 2 . B 2 ] = [(6, 3, 4, 8), (-4, 3, 2, 7)] two 1x4 row matrix 
intervals. 

A - B = [Ai, B,] - [A 2 , B 2 ] = [A, - A 2 , B, - B 2 ] 

= [(3, 2, -1. 0) - (6, 3. 4. 8), (8, 5. 1. 2), (-4, 3, 2, 7)] 

= [(-3, -1.-5, -8), (12, 2, -1,-5)]. 



Consider B - A = [A 2 , B 2 ] - [Ai, B i] 

= [96, 3, 4, 8), (-4, 3, 2, 7)] - [(3, 2, -1, 0), (-4, 3, 2, 7) 
-(8,5, 1,2)] 

= [(3, 1,5, 8), (-12, -2, 1,5)]. 

A-B^B-A. 

Thus the operation subtraction is non commutative. 

It is easy to check A - (B - C) = (A - B) - C for any 1 x n 
row matrix interval. 

Consider A = (Ai, Bi), B = (A 2 , B 2 ) and C = (A 3 , B 3 ) be 
three 1x3 row matrix interval. 



A = (A,, B,) = [(3, 2, 0), (-2, 1,5)], 

B - (A 2 , B 2 ) = [(-3, -2, 4), (5, 3, -1)] 
and C = (A 3 , B 3 ) - [(2, 4, 5), (0, 1,2)]; 

(A - B) - C = ([(3, 2, 0), (-2, 1,5)]- [(-3, -2, 4), (5, 3, -1)] 
-[(2, 4, 5), (0, 1,2)] 

= [(6, 4, -4), (-7, -2, 6)] - [(2, 4, 5), (0, 1, 2)] 

= [(4, 0, -9), (-7, -3, 4)]. 
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A - (B - C) = [(3, 2, 0), (-2, 1, 5)] - ([(-3, -2, 4), (5, 3, -1)] 
-[(2, 4, 5), (0, 1,2)]) 

= [(3, 2, 0), (-2, 1,5)]- [(-5, -6, -1), (5, 2, -3)] 

= [(8,8, 1), (-7, -1, 8)]. 

Thus (A - B) - C * A - (B - C). 

We now proceed onto define the notion of division of two 
1 x n row matrix intervals. 

Division of an 1 x n row matrix interval A = [Ai, B | ] can be 
divided by B = [A 2 , B 2 ] if and only if no coordinate in A 2 and B 2 
is zero: 

That is if A = [A,, BJ = [(ai, a 2 , a n ), (b 1? b 2 , b„)] 

and B = [(c,, c 2 , .... c n ), (d lf d 2 , d n )]. 

Then A/B = [A,/A 2 , B,/B 2 ] 

= [(ai/Ci, a 2 /c 2 , a n /c n ), (bi/dr, b 2 /d 2 , b n /d n )] where c, * 
0 and di ^ 0 for 1 < i < n. 

We see clearly A/B ^ B / A even if none of the coordinates 
in Ai and Bi is zero. 

Also the division operation is non associative. 

Now we can give better structure for the 1 x n row matrix 
interval. 

THEOREM 5.1: Let V = { M - [A, B] the collection of all 1 x n 
row matrix intervals with entries from R or Q or Z or Z„ or C or 
C(Z„)}; V is a commutative ring with unit and has zero divisors. 

The proof is direct and hence is left as an exercise to the 
reader. 

Note V in general is an infinite ring. V is a finite ring only 
the entries of M are from Z n or C(Z n ). 
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We will give one or two examples of this ring. 

Example 5.12: Let 

M={P = [A, B] = [(a!, a 2 , a 3 ), (b b b 2 , b 3 )] I a ; , b ; e Z 2 ; 1 < i < 3} 
be the collection of all 1 x 3 row matrix intervals. Clearly M is a 
finite commutative ring with zero divisors. Further IMI = 64. 

Example 5.13: Let 

P = {M = [A. B] = [(ai, a 2 ), (b b b 2 )] where a i; b; e Z, 1 < i < 2} 
be the collection of all 1 x 2 row matrix intervals P is a ring. P 
has no units but zero divisors. P is of infinite order. 

Example 5.14: Let K = {M = [A. B] = [(a b ..., aio), (b b ..., 
bio)] I a„ bi e C, 1 < i < 10} be a collection of all 1 x 10 row 
matrix intervals, K is a ring of infinite order. 

We can define substructures on these structures. This is 
simple and hence is left as an exercise; these rings have subrings 
and ideals. 

Now we bring in the representation / connection between 
1 x n row matrix intervals and 1 x n interval matrices. 

We know N = ([a b bi], ..., [a n , b n ]) is a 1 x n row interval 
matrix and M = [A, B] 

= [(a b ..., a n ), (b b ..., b n )] is a 1 x n row mattix interval. 

We see every M can be made into N and vice verse. 

Consider N = ([a,, bj], . . ., [a n , b n ]) = [(a b . . ., a n ), (b b . . ., b n )] 

(by taking the first components of every interval in N together 
and the second components of every interval in N together. 



Consider M= [(a b ...,a n ), (b b b n )] = ([a b bi], ..., [a„, b n ]). 
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Write M as n intervals by taking the first component of A 
and first component of B to form the first interval [ai, bi] and so 
on we see M = [A, B] becomes N and N becomes M. 

Now we proceed onto define the notion of n x 1 column 
matrix interval. 

Definition 5.3: Let 



M = [A, B] = 



be a pair of column matrices. We define M to be a m x 1 
column matrix interval. 



From now on wards if we say a m x 1 column matrix 
interval we mean a pair [A. B] where A and B are mxl column 
matrices (vectors). 

We will illustrate this situation by an example. 

Example 5.15: Let 
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M is a 7 x 1 column matrix interval. 




166 



Innovative Uses of Matrices 



Example 5.16: Let 




be a 5 x 1 column matrix interval. 



Example 5.17: Let 
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be a 6 x 1 column matrix interval. 

We can add only two column matrix intervals of same 
order; just as in case of usual mxl column matrix interval. 

We will just show how addition of two mxl column 
matrix intervals is obtained. 




Interval Matrices and Natural Class of Intervals | 167 

Example 5.18: Let 




be two 7x1 column matrix intervals. 
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is again a 7 x 1 column matrix interval. 
Example 5.19: Let 
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is a 5 x 1 column matrix interval. 

We see addition of two n x 1 column matrix intervals is 
commutative and associative. 




is the zero column matrix interval. 



Further for every 




Inview of all this we have the following theorem. 
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Theorem 5.2: Let 
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a„ bi eZ or R or C(Z„) or Q or Z„ or Cj 



be the collection of all n x 1 column matrix intervals. V is an 
abelian group under addition. 

The proof is simple and hence is left as an exercise for the 
reader. 

Clearly if M and N be two n x 1 column matrix intervals 
then MN cannot be defined or found. 

Example 5.20: Let 



M=1 



/T„ “I 



P = 






where a ; , b ; e Z 3 , 1 < i < 4 } 



be the collection of all 4 x 1 column matrix interval. M is an 
abelian group under addition of finite order. 

Example 5.21: Let 



W = 



r r_ i 



S = 



r b^ 



where aj, bj e Z, 1 < i < 7 } 



be a 7 x 1 column matrix interval group of infinite order. 
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Now having seen examples of m x 1 column matrix interval 
group. We now proceed onto define m x n matrix intervals 
(m * n). 

We will illustrate this by some examples. 





f 


a n 


• a ,„ 




bn • 


• b ln 


\ 


Let P = 




a 21 


• a 2„ 


9 


b ?.i 


■ b 2n 






V 


_ a ml • 


• a mn_ 




_ b ml ' 


b 

mn _ 


/ 



where a,j, by e Z or Q or R or C or C (Z n ) or Z n . 1 < i < m and 1 
<j<n}. We say P is a m x n matrix interval. 

We will illustrate this by some examples before we define 
some operations of these sets. 



r 


"3 


0 


1 


4" 




"-3 


5 


1 


2" 


\ 




5 


3 


4 


0 




0 


1 


2 


0 






-1 


2 


3 


4 


9 


1 


2 


0 


4 






7 


0 


8 


4 




5 


1 


6 


2 




V 


0 


1 


0 


2 




0 


0 


1 


-1 


y 



P is a 5 x 4 matrix interval. 

We can as in case of row / column matrix intervals add any m x 
n matrix interval. 

Consider 



r 


“3 4 0" 




1 

1 

K> 

oo 


\ 




2 0-1 




0-12 






0 4-3 


9 


1 0 1 




V 


7 1 0 




6 4 0 


y 



M = 



and 
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-2 14 1-2 3 

0 3 1 -4 5 -6 



2 0 4 7 -8 -9 

5 10 0 1-2 



be two 5x5 matrix intervals. We find M + N as follows: 



M + N = 



3 4 0 -2 8 4 

2 0-1 0-12 



0 4-3 1 0 1 

7 1 0 6 4 0 



-2 14 1-2 3 

0 3 1 -4 5 -6 



2 0 4 7 -8 -9 

5 10 0 1-2 



3 4 0 



7 1 0 



-2 1 4l f-2 8 4 



2 0 -1 0 3 1 0 -1 2 -4 5 -6 

0 4 -3 2 0 4 ’ 1 0 1 + 7 -8 -9 



1 0 6 4 0 



1 -2 3 



0 1 -2 



15 4-16 7 

2 3 0 -4 4 -4 



2 4 1 



-8 10 



12 2 0 6 5 -2 



is again a 4 x 3 matrix interval. Thus we can add matrix 
intervals of same order. 



In view of this we have the following theorem. 
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Theorem 3.3: Let V = {[[aij] mxn , [bij] mxn ] where a, p b r , eZ„ or 
C(Z„) or Z or Q or R or C, 1 <i < m, 1 <j < nj be the m x n 
matrix interval. V is an additive abelian group. 

If elements of V are in Z n or C(Z n ) then V is finite otherwise 
V is infinite. 

However multiplication is not defined on V as 
multiplication is not compatible. 

Example 5.22: Let 

f a. l a 2 ^| ( b, b 2 ^j 

V=j a 3 a 4 , b, b 4 a i( b 3 g Z 5 , 1 < i, j < 6} 

v a 5 a 6y v^5 



be the additive 3x2 matrix interval. V is a group of finite order 
and V is commutative. 

Example 5.23: Let 

v — 

a;, bj g Z, 1 <i, j < 8} 

be the additive 2x8 matrix group interval. M is of finite order 
and is commutative. 

Example 5.24: Let 






and T = 




b i 

di 




be two 2x2 matrix intervals. 
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We define 



P + T = 



a b 
c d 



e f 

U h . 



A c i 



\ 

9 j 


" e i 


f, y 


) 


V§1 


h J. 



cl + 3 .^ 


b + bj" 


^e + ej 


f + f, V 


_l C + C l 


d + di/ 


’U + gi 


h + hj_ 



is again a 2 x 2 matrix intervals. 



r 


Y a b^ 




^e 


f V 




r 


"a. 


b,' 




" e i 




“\ 


V - 


,v c d v 


9 | 


vg 


h y. 


-J 


V- 


_V C 1 


d iy 


9 j 


v§l 


h i ) 


Jy 






*i J 



g h 









gi h 



o 

i7 



f aaj +bc. 


ab. 


+ bd," 


' ee i + f gi 


ef, +fhy 


i^ca, +dCj 


cb, 


+ dd i/ 


\g e i +hg! 


gf, +hhj_ 



is again a 2 x 2 matrix interval. 



1 

JX 




f e f V 

c i L i 




^a b^ 






I 

_o 

p- 


9 j 


vSi h i>. 




_ v c dj 


9 | 


vg h y 



^aa; + b,c 

v Cja + djC 



ajb + b,d^) f ae + bg af + bh 'i 
be, +djd j ^ce + dg cf + dhj 



Clearly PT ^ TP. 

Thus the product of two square matrix intervals is in general 
non commutative. 

Further we see the product of square matrix intervals can 
also be zero. 
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Example 5.25: Let 



P = 



T3 0 
2 1 





1 


0 


2" 




9 


0 


2 


1 






1 

K> 


0 


1 


y 



and 





cn 

0 

1 




"l 0 2" 


\ 




1 2 0 


9 


0 3 1 




V 


0 1 2 




6 1 0 


y 



be any two 3x3 matrix interval. 



PR = 



3 0-1 
2 1 0 
0 2 1 

^3 0-1 
2 1 0 
v 0 2 1 



V 



1 0 2 
0 2 1 
2 0 1 

0 3 O 

1 2 0 

0 1 2 j 



0 3 1 

1 2 0 
0 1 2 



V 

is a 3 x 3 matrix interval. 
Consider 



yv 

Y 0 8 O 
1 8 2 
2 5 2 



yv 



^0 3 O 
1 2 0 



0 1 2 



r \ 0 2 a 
0 3 1 



v 



6 1 0 



r 5 

2 

0 



1 0 2 
0 3 1 
6 1 0 



f \ 0 2 a 
0 2 1 

V 



2 0 1 



f 17 2 2 a 
6 7 2 



1 4 



O 0 2 a 
0 3 1 
v 6 1 0 y 



0 -1 



\ 


ri 


0 


2" 




9 


0 


2 


1 




/ 


y 


0 


K 





RP = 
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Y 0 3 A 
1 2 0 



V 



0 1 2 



A 0 -0 
2 1 0 



7 V 



0 2 1 



A 0 2 
0 3 1 



v 



6 1 0 



/v 



1 0 2 V 
0 2 1 
2 0 1 , 





"6 


5 


r 




f 5 


0 


4 " 






7 


2 


-i 


9 


2 


6 


4 






V 2 


5 


2 , 




v 6 


2 







is again a 3 x 3 matrix interval. However RP ^ RP. 
Example 5.26: Let 



M = 





'0 


2 


1 


0 




"4 


0 


1 


o' 






4 


0 


0 


1 




0 


2 


0 


4 






1 


0 


2 


0 


9 


1 


1 


0 


0 






A 


3 


4 






A 


0 


1 


°, 





and 





r-i 


0 


2 


A 




0 


1 


0 


2 




i 


1 


0 


0 




vO 


1 


0 


A 



"0 


1 


0 


A 




2 


0 


1 


0 




1 


0 


2 


0 




A 


2 


0 


A 





be two 4x4 matrix intervals. 



Now N + M= 





r-i 


2 


3 


A 




f 4 


1 


1 


A 






4 


1 


0 


3 




2 


2 


1 


4 






2 


1 


2 


0 


9 


2 


1 


2 


0 






A 


4 


4 


A 




A 


2 


1 


A 





is again a 4 x 4 matrix interval. 
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Consider MN = 

7 o 2 i oY-i 02 1V4 0 1 0Y0 1 0 lY 

4001 0 102 02042010 

1020 1 100 ’ 1 100 1020 
[o 3 4 o|o 1 0 lj Y 0 1 oJ[o 2 0 lj 




3 0 4 

1 8 5 

2 2 1 



14 2 4 

4 8 10 4 

2 111 



4 7 0 6 1 0 2 0 



is a 4 x 4 matrix interval. 

Now NM 

1 02 1Y0 2 1 0V0 1 0 1 Y 4 0 1 oY 

0 1024001 20100204 

1 100 1020 ’ 1020 1 100 
^0 1 0 iJy 3 4 oj {0 2 0 iJy 0 1 oj 



2 1 7 0 (0 2 1 4 
4 6 8 1 9 1 2 0 
4 2 1 1 ’ 6 2 1 0 



4 3 4 1 0 4 1 



is a 4 x 4 matrix interval. 

We see MN ^ NM. We have the following theorem. 

THEOREM 5.4: Let V = [[ciij] mxm , [bij] mxm ] be a collection ofm 
x m interval matrices with ay, by g Z„ or C(Z n ) or Z or Q or R 
or C. V is a non commutative ring with respect to matrix 
interval addition and multiplication. 
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This proof is direct and simple and hence left as an exercise 
for the reader to prove. 

Clearly V is a ring with zero divisors units and idempotents. 

Further we see as in case of row matrix interval we can in 
case of square matrix interval also obtain the square interval 
matrix. Also from a square interval matrix find the square 
matrix interval. 

This we will illustrate by some examples. 

Example 5.27: Let 




be a 4 x 4 matrix interval. We have V = [V), V 2 ] = [(a,j), (by)]. 
Now we can write V = ([ay, by]). 

[5,3] [0,-5] [1,1] [2,2]" 

[3.0] [5,15] [4,-1] [1,5] 

[2.1] [3,2] [0,3] [4,0] 

J0,0] [1,1] [6,0] [0,7], 

'[5,3] [0,-5] 1 2 " 

[3.0] [5,15] [4,-1] [1,5] 

[2.1] [3,2] [0,3] [4,0] 

v 0 1 [6,0] [0,7], 

is clearly a 4 x 4 interval matrix where intervals in V are 
increasing or decreasing or degenerate intervals. 
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Thus given a 4 x 4 matrix interval we can convert it into a 
4x4 interval matrix. 

Now we just show how an interval square matrix M can be 
written as a square interval. Consider 



[0.1] 


[2,3] 


[5,0] 


[U] 


[2,-1] 


[4.2] 


[0,0] 


[7,1] 


[8,0] 


[0,5] 


[U] 


[2,4] 


[3,3] 


[1,-1] 


[8,9] 


[7.1] 


[4,3] 


[1,2] 


[10.0] 


[4,4] 


[8,5] 


[7,7] 


[2,0] 


[5,5] 


[3,7] 



be a 5 x 5 interval matrix. 



Write M = 





'0 


2 


5 


1 


2" 




f 1 


3 


0 


1 


-r 






4 


0 


7 


8 


0 




2 


0 


1 


0 


5 






1 


2 


3 


1 


8 


9 


1 


4 


3 


-1 


9 






7 


4 


1 


10 


4 




1 


3 


2 


0 


4 








7 


2 


5 


3J 




,5 


7 


0 


5 


7 , 





= [Mi, M 2 ]; M is a 5 x 5 matrix interval. Thus we have the 
following theorem the proof of which is direct and simple. 

THEOREM 3.5: Let M = ([my, riy]) be a n xn interval matrix. 
Then M = [(my), (ny)[ is a n xn matrix interval. Conversely if A 
= [(ay), (by)] is a n xn matrix interval then A - ([ay, by]) is a 
n xn interval matrix. 

The proof is left as an exercise to the reader. 

Thus we see we can have in case of square matrix intervals 
the notion of determinant of intervals. If M = [Mi, M 2 ] where 
Mj and M 2 are n x n matrices then the determinant interval M is 
denoted by IMI = l[Mj, M 2 ]l = [IMil, IM 2 I], 
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We will first illustrate this situation by some simple 
examples. 

Example 5.28: Let 



73 


<7 




' 4 


lY 


! 5 


V 


9 | 


v-2 


V. 



be a 2 x 2 matrix interval. 

Determinant interval of P denoted by 





'"3 4' 




" 4 r 






.7 7 


9 


“2 7 





" 3 4 


4 1 " 


5 7 ' 


’ -2 7 



= |1, 30] is again an interval matrix. 
Example 5.29: Let 





(3 1 7 




"7 0 1 " 






0 5 4 


9 


0 1 0 






v 2 0 Ij 




v 5 0 -7 





be a 3 x 3 matrix interval. 









"3 7 r 




"7 0 I s 


IAI = 






0 5 4 


9 


0 1 0 








v 2 0 Ij 




7 0 -1 
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3 7 1 7 0 1 
0 5 4,0 1 0 
2 0 1 5 0 -1 



5 4 


0 


4 


1 0 


5 


- 3 


-7 


+ 1 




0 1 


2 


1 


2 


0 



1 0 


0 


0 1 


1 0 1 1 


7 


-0 


+ 1 


’-H 

1 

O 


5 


-1 \ 


5 0 



= [61, -12] is the interval determinant. 

Now we show how inverse of an interval matrix can be 
found. 

Suppose M = [Mi, M 2 ] is a n x n matrix interval, then 
M _1 = [M[',M“ 1 ]. 

Thus MM 1 = [MjM^M.M' 1 ] 

= [M[ 1 Mj,M“ 1 M 2 ] 

= m~'m 

= [In, In]. 

We can use any known method to determine the inverse of a 
n x n matrix interval. 

We can use Gauss - Jordan method to find inverse for a 
matrix interval. 
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Suppose 




= [A, B] = [(a y ), (by)] then the rows of A are denoted by Rn, 
R12, Rin and the rows of B are denoted by R 2 i, R 22 , •••, R211 
we write M as 




We by row operations on A and B make A and B to I nxn 
identity matrices and all the operations are done on A and B are 
also simultaneously performed on I nxn . 

The adjoined matrices l nxn turns to be A 1 and B 
respectively. That is M 1 = [A -1 , B ']. 

For instance if 



M = [A, B] = [(ay), (by)] 
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'3 i r 




1 5 2' 






0 1 2 


9 


0 0 2 






V 0 0 Ij 




v 0 1 3J 





to find NT 1 . 



M~ 



f ?> 1 11 0 0 A 

0 12 0 10 



r \ 5 2 1 0 0 A 
0 0 2 0 1 0 





Ll° 


0 


1 


0 


0 


ij 


1° 


1 


3 


0 


0 








R 11~ R 13 




- 


"3 


1 


0 


1 


0 


-1> 




fl 


5 


0 


1 


-1 


0" 


-i 


R n 


V 




0 


1 


9 


0 


1 


0 




0 


0 


9 


0 


1 


0 




R 21 ~ R 22 


r — 




z. 


9 


Z, 




R 21 




- 




0 


1 


0 


0 








1 


3 


0 


0 


K 


- 


R 12 — 2 R i3 




r 


"3 


1 


0 


1 


0 


-n 




fl 


5 


0 


1 


-1 


0" 


i 


R 12 


V — 




0 


1 


0 


0 


1 


-2 




0 


0 


1 


0 


1/2 


0 




R 21 /2 


7 — 




’ 












0 


1 


0 


0 


K 




v0 


1 


3 


0 


0 


K 







"3 


0 


0 


1 


-1 


+f 


fl 


5 


0 


1 


-1 


03 






0 


1 


0 


0 


1 


-2 


0 


0 


1 


0 


1/2 


0 






,0 


0 


1 


0 


0 






1 


0 


0 


-3/2 


K 







1 0 


0 


1/3 


-1/3 


1/3' 


fl 


0 


0 


1 13/2 -5' 






0 1 


0 


0 


1 


-2 


0 


1 


0 


0 -3/2 1 






,0 0 


1 


0 


0 


1 , 




0 


1 


0 1/2 0 





Now M 1 





'1/3 


-1/3 


1/3" 




f 1 


13/2 


-5" 






0 


1 


-2 


9 


0 


-3/2 


1 








0 


1 , 




,0 


1/2 


o. 





MM -1 = 



'3 


1 


r 


'1/3 


-1/3 1/33 


0 


1 


2 


0 


1 


-2 


0 


0 


1 , 


0 


0 


1 , 



ri 


5 


2' 


n 


13/2 


-5> 




0 


0 


2 


0 


-3/2 


1 






1 


3, 


,0 


1/2 
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r l 0 0 A 
0 1 0 
0 0 1 



V 



f l 0 0^ 
0 1 0 
0 0 1 



= [I 3 , I 3 ]. 



Also 





A/3 


-1/3 


1/3' 


'3 


1 


r 


f 1 


13/2 


-5' 


fl 


5 


2^ 






0 


1 


-2 


0 


1 


2 . 


, 0 


-3/2 


1 


0 


0 


2 






l 0 


0 


1 , 




0 


K 




1/2 


0 , 


,0 


1 









f 1 


0 


0 " 




0 


1 


0 




,o 


0 


K 



fi 


0 


o' 




0 


1 


0 






0 


K 





- [I3, 13]- 



Thus MM 1 = M 'M. M 1 is the inverse matrix interval of 
M. 



We can using Gauss Jordan method to find the inverse 
matrix interval. 

Now we can use the method of row reduction and bring a 
m x n matrix interval M to row-reduced echelon matrix interval. 

(i) If the m x n matrix interval M = [Mi, Ni] is row 
reduced 

(a) If the first non zero entry in each non zero row of 
Mi and Ni is equal to 1 in M. 

(b) Each column of Mi and Ni in M which contains 
the leading non zero entry of some row has all its 
other entries 0. 

(ii) Every row of Mi and Ni which has all its entries 0 
occurs below every row which has a non zero entry. 
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(iii) If rows R lh R 12 , R lsj and R 2 i, R 22 , R 2s2 are 

non zero rows of Mi and Nj respectively and if the 
leading non zero entry of rows p and q occurs in k p 
and t q respectively of M and Ni; p = 1,2, Si and 
q = 1, 2, s 2 with ki < k 2 < ... < k S| and tj < t 2 < 



We can also describe an m x n row reduced echelon matrix 
interval M = (Mi, Ni) as follows: 



Here M = [(my), (n,)] = [M,, NJ. 



Either every entry in Mi and Ni is 0 or there exists positive 
integers si, s 2 ; 1 < Si, s 2 < m and si and s 2 positive integers, ki, 

k , and ti, t 2 , ..., t with 1 < kj, t; < n and 

S 1 s 2 J 



(a) 


my = 0 and n pv = 0 if i > si and my = 0 if j < kj and 
p > s 2 and n pv = 0 if v < t v . 


(b) 


m iki = Sy; 1 < i < si, 1 < j < si n pt = 8 pv ; 1 < p < s 2 ; 
1 < v < s 2 . 


(c) 


ki < k 2 < . . . < k and ti, t 2 < ... < t . 

t, 2 


We just give examples of 3 x 5 row reduced echelon matrix 
interval R = [R 1? R 2 ] 



0 


1 


-2 


0 


1/3"] 




"0 


1 


-5 


0 


1/7" 


0 


0 


0 


1 


5 




0 


0 


0 


1 


2 


0 


0 


0 


0 


o J 




v° 


0 


0 


0 


o. 



The following theorem is simple and hence is left as an 
exercise to the reader. 

THEOREM 5.6: Every m xn matrix interval M = (Mu Nj) is 
row equivalent to a row reduced echelon m x n matrix internal. 
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P(x) is the polynomial intervals we have introduced the 
notion of interval polynomials in [55] and have discussed their 
related properties. 

Suppose we have two sets of m-linear equations in n 
unknowns. 

say An xi + A 12 x 2 + ... + A ln x n = yi 
A 2 i Xi + A 22 x 2 + . . . + A 2n x n = y 2 



A m iX] + A m2 x 2 + . . . + A mn x n — y m 

where y iy 2 ... y m and Ay; 1 < i < m, 1 < j < n are given elements 
of a field F. 

Any n tuple (xj, ..., x n ) of elements of F which satisfies 
each of the equations is a solution of the systems. 

Consider another set of m-linear equations in n-unknowns. 



say 



Bn 


x; 


+ Bi 2 x 2 


+ . 


. . + B ]n 


= Zi 


B 2 i 


x; 


+ B 22 x ' 2 


+ . 


■ • + B 2n x' 


= Z 2 


B m l 


x; 


B rn2 X2 


+ .. 


• + B mn x n 


= Z m 



where Z\, z 2 ,.... z m and By; 1 < i < m, 1 < j < n are given 
elements of F. 

Any n tuple (xj, ..., x [ ) of elements of F which satisfies 
the equations is called the solution of the systems. 

If yi = y 2 = ■■• = y m = o and Zi = z 2 = ... = z m = 0 is 
homogeneous set of equations or that each of the equations is 
homogeneous. 

We now define the notion of interval system of m linear 
equations in n interval unknowns. 
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Consider 



1 

_> 

7* 


+ . 


■■ + A m x „ 


B n x ! 


+. 








Yi 




z l 




A 21 X 1 


+ . 


■ + A 2„ X n 


b 2 , x ; 


+. 


■ + B 2n X :, 


= 




y 2 


9 


z 2 




_ A ml X l 


+ . 


■+ A „ m X n’ 


B m , x ; 


+. 


• + B mX_ 






i 

. 

3 




_ Z m _ 





(5.1) 



(5.1) is a system of m linear interval equations in interval n 
unknowns. A row matrix interval [(xi, x 2 , ..., x n ), (xj,x 2 , 
x' n )] of elements from the field F which satisfies each of 
equations (5.1) is called the solution interval of the system. 

We will first give some examples. 



Consider 



3xj - x 2 + x 3 
x, + 3 x 3 

-4x, + 2x, -5 x 3 



-Xj + x 2 + 7x 3 
x 2 -4x 3 
7x, -2x 2 



7 

3 

-1 



2 

0 



-8 



is an 3-linear interval equations in 3 unknowns. 

We give yet another example. 

Take 



x, -3x 2 +x 3 -4x 4 + x 5 


-3x^ + 4 x 3 - x 5 






4 




' 1 “ 




8x, + x 2 + x 4 -5 x 5 


2Xj -4 x 2 +x 4 + 2x 5 






2 




2 




7x 2 + 5x 3 - 4x 4 + 2x 5 


Xj +3x 3 -x 4 + 5x, 






3 


9 


3 




x 3 +4 x 5 


X l+ X 2+ X 5 






0 




-2 
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this is again a 4-linear interval equations in 5 unknowns. 
Yet we consider another example. 



5xj - x, + x, + 2 x 4 


x, + 3x 3 -4x 4 






"3" 




f (T 




3Xj + x 3 -2 x 4 


2x, + x 3 + x 4 






-1 




2 




5x, -x, -3 x 4 


3x 3 - x, + x, 


= 




2 


9 


1 




7x, + x 2 + 8x 3 


3xj + x 2 - 4x 3 + x 4 






1 




4 




-2Xj + 5x 3 - 2x 4 


5x 2 -3x 3 +2x 4 






, 0 , 




- 5 , 





is again a 5-linear interval equations in four unknowns. 

Now we can give these linear interval equations a matrix 
interval equations representations which is as follows. 



A n 


A12 


A 

•• A ln 




B n 


B 12 • 


•• B ln ^ 


A 2) 


A 22 


- A 2n 


9 


B 21 


b 22 . 


•• B 2n 


A ml 


A m2 • 


•• A m j 




V B ml 


B m2 • 


•• B^J 



([(Xj. X n ), (x,\ x' )]) 





Yi 




z i 






y 2 


9 


z 2 






_y m _ 




_ Z m _ 





Here [A. B] [X, X'] = [Y. Z] where A = (Ay), B = (By) 

X= (x l5 . . x n ) and X' = ( xj , . . x^ ). 

We call [A. B] the matrix interval of coefficients of the 
system. 
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Thus this m x n matrix interval with entries from the field F 
is a function [A, B] from the set of pairs of integers (i, j); 1 < j < 
m; 1 < j < n, into the field F. The entries of the matrix interval 
[A, B] are scalar intervals [A(i, j), B(i, j)] = [Ay, By] and quite 
often it is convenient to describe this matrix interval by 
displaying it entries in a rectangular interval array having m row 
intervals and n column intervals. [X, X'] is a n x 1 matrix 
interval and [Y, Z] is a m x 1 matrix interval. 

Thus [A, B] [X, X'] = [Y, Z] or [AX. BX'] = [Y, Z] is 
nothing more than a shorthand notation of the system of linear 
interval equation. 

Just like in usual matrices we can define row equivalent 
m x n matrix interval [A. B], [A', B'] such that [AX, BX'] = 
[(0), (0)] and [A'X, BX'] = [(0), (0)] have exactly the same 
interval solutions (solution intervals). 

The proof of this result is direct and can be derived by any 
interested reader. 

We can as in case of usual simple matrices say a n x n 
matrix interval M = [A, B] is an invertible n x n matrix interval 
and N = [C, D] is another n x n invertible matrix interval then 

(1) M' 1 = [A. BF 1 = [A -1 , B '] 

and (M -1 ) -1 = [A~‘ B -1 ] -1 
= [A, B] (since (A -1 ) -1 = A. 
and (B- 1 )- 1 = B ). 



If both M = [A. B] and N = [C, D] are invertible n x n 
matrix interval then so is MN and (MN) 1 = N 'M 

Further it can be easily verified that product of n x n 
invertible matrix intervals is invertible. Further if A = [A l5 BJ 
is invertible matrix interval; that is A row equivalent with 
I = [I n , I n ] the n x n identity matrix interval. 




190 



Innovative Uses of Matrices 



If A = [Ai, Bi] is an n x n invertible matrix interval then 
AX = [(0), (0)]; that is [A^, B^J = [(0), (0)] has only the 
trivial solution X = [Xi, X 2 ] = [(0), (0)]. 

The system of interval equations AX = Y 



[A 1 ,B 1 ][X 1 ,X 2 ] = [Y 1 ,Y 2 ]that is 



[ A i X ] , B,X 2 ] = [Y lf Y 2 ] has a solution 

X = [Xj, X 2 ] for each n x 1 matrix interval Y = [Y h Y 2 ], 

We can define as in case of usual square matrices define 
these notions in case of square matrix intervals. 

Let A = [Ai, Bj] be a n x n matrix interval over the field F. 

The principal interval minors or minor intervals of 
A = [A], B i] = [(Ay), (By)]; 1 < i, j < n are scalars 



A k (A) = [A k (A0, A k (B i)] defined by 
A k (A) = [A k A], A k A 2 ] 





1 

> 


> 

** 




1 

w 


_tD 

£ 




det 


... > 

to 


• A 2k 


, det 


B 2 i • 


• B 2k 






1 

> 

7T 


> 




1 

dd 

** 







where 1 < k < n. 



Let A = [Ai, A 2 ] be a n x n matrix interval we say A is a n x 
n upper triangular matrix interval if both A, and A 2 are upper 
triangular matrix interval. 



Consider M = [Mi, Ni] 
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( M 

iVljl 


m 12 . 


- M, n ^ 




( N 
1 mi 


n 12 


n 13 . 


- N ln ^ 


0 


M 22 


.. M 2n 




0 


n 22 


N 23 . 


- N 2n 


0 


0 . 


- M 3n 




0 


0 


n 33 . 


- N 3n 


0 


0 . 


- M 4n 


9 


0 


0 


0 . 


- N 4n 


V 0 


0 . 


- M nny 




v 0 


0 


0 . 


- N nny 





(5 


7 


2 


o N 




' 2 


1 


0 


4 " 






0 


1 


-5 


-7 




0 


3 


1 


-2 






0 


0 


3 


2 


9 


0 


0 


4 


0 








0 


0 


K 




vO 


0 


0 


7, 





is a 4 x 4 upper triangular matrix interval. 



Take A = [Ai, Bi] be a n x n square matrix interval. If both 
Ai and Bi are n x n lower triangular matrices then we define A 
to be a n x n lower triangular matrix interval. 



We will just give an example before we proceed onto define 
further structures. 



Let 



P = [Pi,P 2 ] = 





f 5 


2 


0 


0 


n 




"7 


6 


4 


0 


2^ 






0 


7 


0 


6 


3 




0 


1 


0 


5 


0 






0 


0 


5 


2 


1 


9 


0 


0 


2 


-2 


2 






0 


0 


0 


6 


-1 




0 


0 


0 


1 


4 






vO 


0 


0 


0 


8, 




vO 


0 


0 


0 


7, 





is a 5 x 5 upper triangular matrix interval. 
Consider P T - |P,. P 2 ] T = [ P, T , P 2 T 1 
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f 5 


0 


0 


0 


0] 




'l 


0 


0 


0 


0^ 






2 


7 


0 


0 


0 




6 


1 


0 


0 


0 






0 


0 


5 


0 


0 


5 


4 


0 


2 


0 


0 






0 


6 


2 


6 


0 




0 


5 


-2 


1 


0 






vl 


3 


1 


-1 


8J 




v2 


0 


2 


4 


V 





we see P 1 is a lower triangular matrix interval. 
Consider A = [A], B|] 



6 0 0 0 

7 10 0 

8 2 4 0 

0 5 12 



2 0 0 0 
0 5 0 0 
12 7 0 
3 0 15 



be a 4 x 4 lower triangular matrix interval. 

% 7 

a' = [ a;,b[] = 



0 12 5 
0 0 4 1 
0 0 0 2 
is a 4 x 4 lower triangular matrix interval. 



0 ^ (2 

0 






3 ^ 

0 



0 0 7 1 
0 0 0 5 



Thus we see by taking transpose we can easily covert a 
n x n lower triangular matrix interval into a n x n upper 
triangular matrix interval and vice versa. 

We can easily verify the following statements are equivalent 
in case of an invertible n x n matrix interval A = [Ai, A 2 ]; 

(a) There is an upper triangular matrix interval 
P = [Pi, P 2 ] = [(^)’( p .i)] with P k ' k = P^ = 1; 1 < k 

< n such that B = AP = [AiPi, A 2 P 2 ] = |B h B 2 ] is a 
lower triangular matrix interval. 
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(b) The principal minors of A = [A], A 2 ]; that is the 
principal minors of A, and A 2 are all different from 0. 

We call A = [Aj, A 2 ] to be a n x n diagonal matrix interval 
if A = [Ai, A 2 ] = [(Ay), (By)] than Ay = 0 and By = 0 if i * j, 
1 < i, j < n. 



Thus [Pi, P 2 ] - 



7 

0 

0 

0 



0 0 
8 0 
0 1 
0 0 





1 

oo 


0 


0 


o 






0 


6 


0 


0 




9 


0 


0 


-1 


0 






1 

o 


0 


0 


1 

-J 





is a 4 x 4 diagonal matrix interval. 



Further it can be easily verified if A = [A], A 2 ] = [(ay), (by)] 
be a n x n square upper triangular matrix interval that is if ay = 
0, by = 0 for i > j, that is every entry below the main diagonal is 
zero. A is invertible if and only if every entry on the main 
diagonal of both Aj and A 2 are different from zero. 



For instance 





'8 


0 


4 N 


/ 




0 


5 


1 


9 






0 




V 



3 1 
0 5 
0 0 



2" 

7 



the upper triangular matrix interval is invertible, where as 





'8 


0 


4 


8^ 




f 1 


2 


0 








0 


0 


5 


9 




0 


1 


2 


3 






0 


0 


7 


4 


9 


0 


0 


0 


5 








0 


0 


2, 




v 0 


0 


0 


-7, 





the upper triangular matrix interval is not invertible. 
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Prove if P = [P,, P 2 ] 



f 1 


1/2 


1/3 


1/4' 




f 1 


1/2 


1/3 


1/4 


1/2 


1/3 


1/4 


1/5 




1/2 


1/3 


1/4 


1/5 


1/3 


1/4 


1/5 


1/6 


9 


1/3 


1/4 


1/5 


1/6 


,1/4 


1/5 


1/6 


1/7. 




,1/4 


1/5 


1/6 


1/7 



is an invertible matrix interval and P 1 = [ P, 1 , P, 1 ] has integer 
entries. 

We can define permutation matrix interval as in case of 
usual matrices. 



We will give an example or two. 



P=[Pi,P 2 ] = 





'0 


1 


0 


o' 






0 


0 


1 


0 






1 


0 


0 


0 


9 






0 


0 


1/ 


V 



1 0 
0 1 
0 0 
1 0 



0 " 

0 

1 



is a permutation matrix interval. 

Now we proceed onto give example of symmetric and skew 
symmetric matrix intervals. 

Consider 



H = [H 1? H 2 ] = 



3-124 
-18 4 0 

2 4 3 9 

4 0 9 2 



8 4 3 1 

4-142 
3 4 0 7 

12 7 2 



we see in this matrix interval a,, = ajj if i ^ j, 1 < i, j < 4. 

Thus if A = [Ai, B i] = [(ay), (by)] be a square matrix 
interval. If in Ai and Bi, ay = ay and by = by if i ^ j, 1 < i, j < n 
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then we define A = [A], Bi] to be a symmetric matrix interval. 
If A = [A], B i] = [(aij), (bjj)] be a square matrix interval. 

If aij - — ajj and by - -b,, for i & j and if i = j then a i; = b„ = 0, 
1 < i < n, 1 < i, j < n, then we define A to be a skew symmetric 
matrix interval. 

Let P = [Pi, Pd 



0 


3 


1 


-4 


5 


8" 


-3 


0 


2 


1 


-3 


0 


-1 


-2 


0 


-7 


2 


9 


4 


-1 


7 


0 


5 


1 


-5 


3 


-2 


-5 


0 


8 


-8 


0 


9 


-1 


-8 


o. 



0 


1 


2 


3 


4 


5 


-1 


0 


-6 


-7 


-8 


-9 


-2 


6 


0 


1 


-2 


3 


-3 


7 


-1 


0 


4 


-5 


-4 


8 


2 


-4 


0 


6 


-5 


9 


-3 


5 


-6 


0 



be a 6 x 6 skew symmetric matrix interval. 

We can as in case of usual square matrices write the square 
matrix interval as a sum of a square symmetric matrix interval 
and a square skew symmetric matrix interval. 



We will illustrate this situation by an example or two. 
Let 



7 -2 
5 0 
0 8 



A = [A,, B ,] = 



f 8 4 3 5 A 

7 12-2 

10 -5 8 0 7 

0 4 0 -7 J -2 



^ 2-4 
0 5 



0 

4 



be a 4 x 4 square matrix interval. 




196 



Innovative Uses of Matrices 



A* = [A,, BJ* = 



>;<b;] = 





"8 


7 


10 


0 ' 




' 2 


0 


7 


- 2 " 




4 


1 


-5 


4 




-4 


5 


0 


4 




3 


2 


8 


0 


5 


3 


7 


5 


0 




v 5 


-2 


0 


-0 




v .1 


-2 


0 


8 . 



be the transpose of the 4x4 square matrix interval. 

Let A = [Pj, P 2 ] + [Q,, Q 2 ] 

= [Pi + Qn ?2 + Q 2 ] 



u „ A i+ A l ^ B l + B i 
where Pi - — 1 L , Qi = — - L 



a.-a; B, -b; 

P 2 = — — — 1 and Q 2 = -4— — 1 



Now Pi - 



' 8 


4 


3 


5 ^ 




"8 


7 


10 


0 s 


7 


1 


2 


-2 


+ 


4 


1 


-5 


4 


10 


-5 


8 


0 




3 


2 


8 


0 




4 


0 


" 7 , 




v 5 


-2 


0 


-0 



06 11 13 5 ^ 

112-3 2 

13 -3 16 0 

5 2 0 -14 



Clearly Pi is a symmetric matrix of order 4x4. 
A - A 1 

Now consider P 2 = — L 
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Y 8 4 3 5 ] (8 7 10 oY 

17 12-2 41-54 

2 10 -5 8 0 3 2 8 0 

^0 4 0 - 1 ) ^5 -2 0 - 1 ) 



"0 - 3-7 5 ' 

1 3 0 7 -6 

_ 2 7 -7 0 0 ‘ 

v -5 6 0 0 , 

Clearly P 2 is a skew symmetric matrix of order four. 



Qi = 



b,+b; 

2 



f f 2-4 3 1 W 2 0 7 2 Y 

1 0 5 7 -2 -4 5 0 4 

27050 3750 

1^-2 4 0 8 J { 1 -2 0 8 J 



' 4 -4 10 -1 

J 4 10 7 2 

” 2 10 7 10 0 

V -1 2 0 16 

It is easily verified that Qi is a symmetric matrix of order 
four. 



Q 2 — 



b,-b; 

2 




198 



Innovative Uses of Matrices 



' 2 


-4 


3 


1 > 




" 2 


0 


7 


- 2 " 


0 


5 


7 


-2 




-4 


5 


0 


4 


7 


0 


5 


0 




3 


7 


5 


0 


.-2 


4 


0 


8 , 




> 1 


-2 


0 


8 , 



2 



0-4-4 3 

4 0 7 -6 

4-700 
-3 6 0 0 



is a skew symmetric matrix of order four. 



Now A = 
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11 


13 


5 'j 




" 0 
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-7 


5 1 
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11 
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-3 
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-6 












+ — 
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13 


-3 


16 
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-7 
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0 




v5 


2 


0 


-14; 




v-5 


6 


0 


0 , 





" 4 


-4 


10 


-f| 




" 0 


-4 


-4 


3 ' 


1 


-4 


10 


7 


2 


1 


4 


0 


7 


-6 


— 










+ — 










2 


10 


7 


10 


0 


2 


4 


-7 


0 


0 




Vi 


2 


0 


16; 




v-3 


6 


0 


0, 



= [A x , BJ 

= (P1+P2), |(Qt + Q 2)]. 

2 2 
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The sign column matrix interval is 



X = [X,, X 2 ] = 



where aj, bj e {-1, 1 } ; 1 < i < n. 





f n 




r-fi 






-i 




i 






i 




i 






i 


•> 


-i 






i 




-i 






-i 




i 















= [Yi, Y 2 ] is a sign column matrix interval or sign column 
vector interval. 

The sign row matrix interval sign row vector interval 
A=[A,.A 2 ] 



= [(ai, a m ), (b 1; b 2 , b m )] where a ; . b ; e {-1. 1}; 1 <i < m. 
Thus 

A = [(1,-1, 1. 1. 1, -1,-1, 1), (-1, 1, 1, -1,-1, 1,-1, 1)] 

is a sign row vector interval or (sign row matrix interval). 

These sort of matrix intervals are useful in factor analysis. 
We will give examples of orthogonal matrix interval. 
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Let 



X = 



1 -7 
1 -2 

1 3 



f 2 O 

-1 2 

V° 1 



7 l 


1 M 


r 2 


-l 


°y 


U 


-2 3/ 


’7 


2 


u. 



be a matrix interval. 

X' = 

is again a matrix interval. 
Consider 

x* x = 1 1 lV2 - 1 



-i 


-2 


3j 


V 


2 


U 










fl 


-r 






i 


1 


o 








-l 








l 


-2 






-i 


(N 

1 








’ 1 


2 






/ 


l 1 


3 , 


V 








r 


"3 


0" 


f 5 


°T 










14, 




sj. 



1 -1 
1 -2 
1 3 



( 



r 

-1 2 

o 1 , 



( 2 O 
-1 2 

0 1 



Thus this is an example of a orthogonal matrix interval. 

If we have X = [A, B] be a matrix interval. 

X 1 = [A 1 , B 1 ] be the transpose of X. 

Suppose X*X = [A‘A. B'B] = [I n , I„] then we call x to be a 
orthogonal interval matrix. 
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Consider 

Y i o Vo oY 
0 0 0 1 
X “ 0 1 ’ 0 0 

.1° °Ju °J_ 

to be the matrix interval. 

f 7i o o oVo o o lV 

X 1 — 

[(o 0 1 oJ\o 1 0 oj 

be the transpose of the matrix interval X. 

Y i o Vo oY 

t Yi o o oVo o o i VI oo oi 

X'X = , X 

^0 0 1 oj ^0 1 0 ojj 0 1 0 0 

Ll° °J U °J. 




= [I 2 , Id. 

Thus X is an orthonormal matrix interval. 

Suppose X = [A. B] be a matrix interval and X 1 = [A 1 . B‘] be 
the transpose of the matrix interval; if X l X = [A l A, B'B ] = [I m , 
I m ] then we say X is a orthonormal matrix interval. 



Interested reader can give examples of them. 




202 Innovative Uses of Matrices 



This result can easily be extend to square matrix interval. If 
X = [A. B] be a n x n square matrix interval and X' = [A*. B‘] be 
the transpose of X be such that X l X = [I n , I n ] then we define X 
to be a orthogonal matrix interval or more specifically X is a 
orthonormal matrix interval. Suppose X l X = [D, C] where D 
and C are just diagonal matrices or symmetric matrices then we 
define X to be only a orthogonal square matrix interval. 

Thus we use matrix theory in a new way to study properties 
of interval matrices constructed using natural class of intervals. 




Chapter Six 



Dfiv[]VkTRix of Refined Labels 



[46, 52] have worked with the new notion of DSm vector 
spaces and DSm super vector spaces. L z = {..., L j, L i, L 0 , 
Li, L 2 , Lj, ...} = {Lj I j e Z} the set of extended labels with 
positive and negative indexes. 



Similarly Lq = L q I q e Q} as set of labels whose indexes 
are fractions. Lq is isomorphic to Q, through the isomorphism 



f<2 (Lq) = 



q 

m + l’ 



for any q e Q. 



On similar lines they define L R = {Lj I j e R}; R the set of 
real numbers. 



L r is isomorphic with R through the isomorphism 

f R (L r ) = — 
m + l 

for any r e R. For more refer [46, 52], 
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Further they have proved {L R , +, x} is a field where + is the 
vector addition of labels and x is the vector multiplication of 
labels defined as DSm field of refined labels. 

If [Li, L 2 ] is a label interval then L 3/2 = Li, 5 is the label in 
the middle of the label interval. 

Also L = -L; that occur in qualitative calculations. We just 
from [ ] recall the operations. 

Let a, b, c e R and the labels 

T _ a T _ b _ c 

L a — , L,k — and L c — . 

m + 1 m + 1 m + 1 

Vector addition of labels 



L a + Lb — L a+b — L a+b . 

m+1 

In this chapter we just indicate how we have used matrix 
theory in the construction of refined label matrices. When the 
refined labels happen to be ordered that is 



Lo < Lj < . . . < L m 

m e N; then we can use these matrices to built vector spaces, 
find eigen values and all properties pertaining to vector spaces 
can be derived. This has been elaborately carried out in [46, 
52], 



So this is an unique and a new way of using matrix theory. 
Not only it is matrix of refined labels whose properties are 
studied, we have also invented to notion of DSm super vector 
spaces, which uses refined labels of super matrices. 

This study is elaborately made in [46, 52], 




Chapter Seven 



n-Mvr 



MK 



and Their Applications 



In this chapter we proceed onto recall the definition, 
properties and applications of n-matrices. This is yet another 
innovative method of using matrices. For more about these 
concepts please refer [47, 51, 54, 56]. 

DEFINITION 7.1: Let A = A] u A 2 where A; and A? are two 

matrices which are distinct and the entries are taken from Q or 
R or Z or Z„ or C(Z„). 

If both A i and A 2 are row matrices we call A to be a birow 
matrix. Similarly, if both Ai and A2 are column matrices we call 
A to be a bicolumn matrix; A is defined as the square bimatrix if 
both Ai and A2 are square matrices. If A / and A2 are 
rectangular matrices we define A to be a rectangular bimatrix. 
If in A - A] u A2 one of A 1 is a row matrix and A2 is not a row 
matrix we define A to be a mixed bimatrix. That is if Aj and A? 
are two different matrices of order m x n and s x t we call A to 
be mixed bimatrix. 



We will illustrate this by some examples. 
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Example 7.1: Let 



A = Ai u A 2 = (3, 1, 5, 6, 0, 2) u (3, 1 0, 0, 5), 



we call A to be a row bimatrix. 



Example 7.2: Let 



A = A! u A 2 = 



be a column bimatrix. 



Example 7.3: Let 



B = Bi u B 2 = (9. 0, 1) u (8, 4. -4), 



B is a row bimatrix. 



Example 7.4: Let 



r n 4 2 7 1 0 

9 0 12 

A = A[UAi= u 0 1 9 -1 8 

" 0 10-4 

L J 5 6 8 1 -4 



A is a rectangular bimatrix. 
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Example 7.5: Let 



A — A[ At — 



3 

-1 



u 



4 

-1 

0 



0 -1 
7 0 

2 0 
1 5 



2 

3 

1 

8 



be a square bimatrix. 



Example 7.6: Let 



M = (8, 0,-1, 4, 3, V3,7) u 



5 

-6 

8 



1 

-4 



be a mixed bimatrix. 



Example 7.7: Let 



T = Tj u T 2 - 




0 2 a 
0 3 



u 



-9 

2 

3 

4 
-5 
6 



be a mixed bimatrix. 
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Example 7.8: Let 



S — S 1 1 — j S 2 — 



9 

2 

3 

0 

-3 

0 




0 

2 



be the mixed bimatrix. 



Example 7.9: Let 



P = P 1 uP 2 = (-3,0, 1, 2) u 



1 

-1 



2 -1 

0 6 



5 7 

-2 0 



be the mixed bimatrix. 

Now we can have bimatrices both mixed or otherwise. 

Suppose A = Ai u A 2 u A 3 u . . . u A n ; oc > n > 2 and A;’s 
are distinct row matrices we define A to be n-row matrix. If 
n = 2 it becomes a row bimatrix. If n = 3 we get the row 
tri matrix and so on. 

If row matrices are replaced by column matrices we call 
them as n-column matrix. Instead of a column matrix we use 
rectangular matrices we call them as n-rectangular matrix or 
rectangular n-matrix. 

Suppose in A = Ai u A 2 u ... uA n the Aj’s are square 
matrices, we define A to be a n-square matrix or a square n- 
matrix. 



We will illustrate all these situations by examples. 
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Example 7.10: Let A = A] u A 2 u A 3 u A 4 u A 5 = (8, -1, 0, 2, 
1) u (0, 5, 8) u (3, 2) u (-7, 1, 2, 3, 4, 5, -7) u (4, 6, 8, 10, 12) 
be a 5-row matrix or row 5-matrix. 



Example 7.11: Let 



M = Mi u M 2 u M 3 u M 4 = 
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be the column 4-matrix. 

Example 7.12: Let V = V) u Y 2 u V 3 u V 4 u V 5 u V 6 u V 7 
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be a 7-square matrix or square 7-matrix. 
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Example 7.13: Let T = Ti u T 2 u T 3 u T 4 u T 5 
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be a rectangular 5-matrix or 5-rectangular matrix. 



Example 7.14 : Let V = V! u V 2 u V 3 u V 4 



= (-9.0. V3.I, 7+i) u 
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-1 

8-i 

5i 
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4 + 8i 



8 -i 7 - i 8 + 4i 
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be a 4-mixed matrix or mixed 4-matrix. 
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Example 7.15: Let A = Ai u A 2 u A 3 u A 4 u A 5 u A 6 



= (-8, 0, 1, 3+2i) u 
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be a 6-mixed matrix or mixed 6-matrix. 



Now we have seen just definition and examples of them. 
We proceed onto give three major applications of them. 

The first application of bimatrices (n-matrices) is in the 
construction of bilinear algebra (n-linear algebra of type 1 and 
type II). For more literature please refer [60-1]. 

One of the major application of these n-matrices (bimatrices 
is that they can be used in the construction of n-eigen values 
(bieigen values) and n-eigen vectors (bieigen vectors). This 
bulk operation helps in saving time as well as easy for 
calculation and comparison. Further bimatrices are used in the 
construction of bifuzzy models like BiFCMs, BiFRMs, n-FCMs 
and n-FRMs; n > 2. 

We just indicate how the bieigen values (n-eigen values) are 
calculated. We only give illustrations for theory please refer 
[60-1], 
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Example 7.16 : Let 
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M = M,uM 2 = 
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be a bimatrix associated with some linear operator. The bieigen 
values associated with M is as follows: 



IM — (A,i vt Ag)I — IMj — ^.[1 IM2 — ^2! 



3-A,j 

8 



0 

-1-^ 



u 



8-A, 2 0 0 

0 -\-X 2 2 

0 6 l-X 2 



= (3 - h) (-1 - h) u ( 8 -X 2 ) [ — ( 1 — A. 2 ) - l 2 ] - 0 u 0 
= - (1-A.t) (3-A.O u -( 8 -A. 2 ) ((1- X\ )+l 2 ) 

= 0u0 



leading to X,i = 3 and -1 and hi = 8 hi = ± v/l3 . 

Thus l = 3u8, 3u± -Jl 3 , 3 u - Vl3 
-1 u 8. -1 u Vl3 , -1 u - yfi3 . 

Thus using these six bieigen values one can calculate the 
related bieigen vectors. By this way we can compare the eigen 
values and also perform the operations simultaneously. 

However the entries of the matrix M are taken from the 
field of reals. If they are taken from the rationals all the 
solution will not exist. We may have partial solution in that 



case. 
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Example 7.1 7: Let P = Pi u P 2 u P 3 
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be the trimatrix with the entries from reals. To find the trieigen 
values and trieigen vectors of P. 

Consider IP - X (=A| u Xi u A, 3 )l 



= IPi-XilulP 2 -X 2 lulP 3 -A,3l 









2-X 2 
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3-X 


1 + 1 
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0 3 


-*2 




|l -X. 
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0 2-X, 0 0 

9 0 3-^ 0 

0 1 0 4 -^ 



- (3-A.O (5-Ai) u (2-a 2 ) (1-Xj) x 

(3-a 2 ) u (I-A 3 ) (2-a 3 ) (3-Xg) (4-A.a) 

Thus the trieigen values are (3, 2, 1), (3, 2, 2), (3, 2, 3), 
(3, 2, 4), (3, 1, 1), (3. 1, 2), (3, 1, 3), (3, 1, 4), (3, 3, 1), (3, 3, 2), 
(3, 3, 3), (3, 3. 4), (5, 2, 1), (5, 2, 2), (5, 2. 3), (5, 2, 4), (5, 1, 2), 
(5, 1, 3), (5, 1. 4), (5. 1. 1), (5, 3, 1), (5, 3, 2), (5, 3, 3), (5, 3. 4). 



We have 24 sets of trieigen values and their corresponding 
24 set of trieigen vectors can be calculated. This is one of the 
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major innovative method of applying the bimatrix (n-matrix) 
theory in bilinear algebra or (n-linear algebra). 

The next set of applications of bimatrices (n-matrices) is in 
the fuzzy models. We just indicate how these bimatrices (n- 
matrices) used. We call a bimatrix A = A] u A 2 to be a fuzzy 
matrix if the entries of both A] and A 2 are from the unit interval 
[0, 1]. We make use of the min max or max min or min min or 
max max operations on them. 

On similar lines we define n-fuzzy matrix. A matrix S = S i 
u S 2 u ... u S n (n > 2) is said to be a n-fuzzy matrix if the 
entries of each Si is from the unit interval [0, 1]; 1 < i < n. We 
use min max or max min or other such operators on Si’s. We 
also at times include -1 as a fuzzy number depending on the 
fuzzy model. The fuzzy models which make use n-fuzzy 
matrices are n- FCMs. n-FRMs and n-FAMs. 

For more about these structures one can refer the books 
[48. 54. 56]. 

Suppose we have a problem related with the school 
students. 

The school student problem is dependent on the teachers 
also dependent on the parents and on the school management. 
So the school children problem is dependent on three groups. 
Thus we can use a fuzzy trimatrix and use the Fuzzy Relational 
Maps (FRMs) model to analyse the problem. 

Let us assume Si, S 2 , ..., S6 are the attributes associated 
with the FRM. 

Let Ti, T 2 , ..., T 7 be the attributes related with the teachers. 

Pi, P 2 , ..., P6 are the attributes associated with the parents 
and Mi, M 2 , ..., M 8 are the attributes associate with the 
management. 
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We get using an experts opinion a trigraph and the matrix 
associated with the trigraph will be the fuzzy trimatrix which 
serves as the dynamical system of the triFRM model. If S 
denotes the triFRM models tridynamical system then S = Si u 

S 2 uS 3 




The entries of the matrices Si, S 2 and S 3 are from the unit 
interval [0, 1]. 

Any state vector X = Xi u X 2 u X 3 = 

— ( x | ... x‘)u(x?...x*)u(xf ... xl) 

where xj e {0, 1}; 1 < i < 6, 1 < j < 3. 

We find X o S = Xi o Si u X 2 o S 2 u X 3 o S 3 
= YiuY 2 uY 3 
= Y. 

where elements of Yj are from {0, 1} Yi will be (y| ... y \ ) Y 2 
will be (yf.y*, yj:) and Y 3 will be ( y 3 , y 3 2 , ..., y 3 s ) yj e 
{0, 1}, 1 <j<3and 1 <i<8. 

Thus these n-fuzzy matrices are used in constructing 
n-fuzzy models, when we have n-experts or these models can 
also be defined as multi expert models. 
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Finally we see these n-matrices when constructed over Z p „ , 

p a prime and co > n > 1 , can be used in the construction of n- 
codes. 

We have defined n-codes [49]. For more about the uses of 
n-matrices in algebraic coding theory refer [49]. 

The construction of n-codes using n-matrices happens to be 
one of the new ways of applying n-matrices to algebraic coding 
theory. 

Thus the concept of n-matrices and their applications are 
innovative contributions in matrix theory by the authors. 
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In this book authors have given the innovative ways in 
which they have used matrices. Starting from ATD, RTD 
& CETD matrices which can be applied to any collected 
data depending on time. A new type of matrices called 
fuzzy linguistic matrices which play a major role in the 
construction of fuzzy linguistic models is defined. 

Then the super matrices are used in the construction of 
super linear algebra and super fuzzy models. Further 
interval matrices using natural class of intervals is 
introduced. The ingenious way of defining DSm matrices f() 
of refined labels is an interesting feature. Finally the 
construction of n-matrices and its use in the building 
of n -codes is described. 
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